//spherical points distribution
//using golden ratio method
phi=(1+sqrt(5))/2;
nph=100; //half number of points
i=-nph:nph;
lat=asin((2*i)/(2*nph+1));
lon=(2*%pi/phi)*i;
x=cos(lat).*cos(lon);
y=cos(lat).*sin(lon);
z=sin(lat);
scatter3(x,y,z);
isoview
