Function sample_sphere_surface 

pub fn sample_sphere_surface(rng: &mut Lcg) -> [f64; 3]

Sample a point uniformly on the surface of the unit sphere (Marsaglia 1972).

Returns [x, y, z] with x² + y² + z² = 1.