
n = 20;   % plots a 2n by n grid on the surface of the sphere
[xc,yc] = meshgrid(linspace(-3,1,2*n+1), linspace(-1,1,n+1));
[xp,yp,zp] = mapc2m(xc,yc);
surf(xp,yp,zp,'LineWidth',2);   
colormap([1 1 1]);   
daspect([1 1 1])
axis off
