pub fn orthonormal_basis<S>(e1: Unit3<S>) -> [Unit3<S>; 3]where S: Real,
Create an orthonormal basis given the first basis component