Source code:
%reset -f #core import math; #libs import numpy as np; import matplotlib.pyplot as pp; from mpl_toolkits import mplot3d; #function to plot Pi = math.pi; f = lambda x,y: math.sin(x/Pi/2)*math.cos(y/Pi/2)*(x-y); #make data Len = 50; X = np.linspace(-10,10, Len); Y = np.linspace(-10,10, Len); X,Y = np.meshgrid(X,Y); Z = []; for I in range(Len): Z += [[]]; for J in range(Len): Z[I] += [f(X[I][J],Y[I][J])]; #any function of x,y #plot pp.figure(figsize=[10,10]); pp3d = pp.axes(projection="3d", elev=50, azim=70); pp3d.set_title("Function z = f(x,y)"); pp3d.set_xlabel("x values"); pp3d.set_ylabel("y values"); pp3d.set_zlabel("z values"); pp3d.plot_wireframe(X,Y,np.array(Z)); #pp3d.plot_surface(X,Y,np.array(Z), cmap="coolwarm"); #eof
No comments:
Post a Comment