X=linspace(-pi,pi,100);
Y=X;
[x,y]=meshgrid(X,Y);
z=cos(0.1+x.^2/5-x/3-y/5+y.^2/9);
surf(x,y,z);