%draw spherical surface
[phi,theta]=meshgrid(0:0.1:pi,0:0.1:2*pi);
size(phi)
size(theta)
X=cos(theta).*sin(phi);
Y=sin(theta).*sin(phi);
Z=cos(phi);
surf(X,Y,Z);