function [x,y,z]=mysphere(n)
teta=0:pi/n:pi;
phi=0:2*pi/n:2*pi;
[T P]=meshgrid(teta,phi);
x=sin(T).*cos(P);
y=sin(T).*sin(P);
z=cos(T);
