[th,ph]=meshgrid(0:pi/10:2*pi,0:pi/10:pi);
x=2*cos(th).*sin(ph);
y=2*sin(th).*sin(ph);
z=2*cos(ph);
surf(x,y,z);
axis equal;