clear all
colormap(jet(256));
[a,b,c]=mysphere(100);
surf(a,b,c);
shading('interp');
