#!/usr/bin/python
from mpl_toolkits import mplot3d
import numpy as np
import matplotlib.pyplot as plt

func=np.loadtxt("electr.dat")
x=np.linspace(0,1,500)
y=np.linspace(0,1,500)
(X,Y)=np.meshgrid(x,y)

#print np.shape(func)

p3f=plt.axes(projection="3d")
p3f.set_xlabel('X')
p3f.set_ylabel('Y')
p3f.set_zlabel('func')

#p3f.plot_surface(X,Y,func,color="black")
#p3f.plot_surface(X,Y,func,cmap='copper')
p3f.contour3D(X,Y,func,50,cmap='copper')
#p3f.plot_wireframe(X,Y,func,color="black")
#p3f.plot_trisurf(x,y,func,color="black")
#p3f.scatter(X,Y,func,color="black")


plt.grid()
plt.show()