mySphere(N)

MATLAB function for generating equidistant points on the surface of a unit sphere.

Input: the target number of points or nodes on the surface of a unit sphere (N).
Output: position column vectors X, Y, and Z of length equal to or close to N.
Points are distributed by choosing circles of constant latitudes separated by d_theta, and on these circles placing points separated by d_phi, such that d_theta ~ d_phi. The area per point (d_theta)*(d_phi) therefore remains constant across the sphere surface.