pub fn velocity(
position: &Vec3,
angular_momentum: Float,
eccentricity: Float,
radius: Float,
parameter_ellipse: Float,
true_anomaly: Float,
longitude_ascending_node: Float,
argument_periapsis: Float,
inclination: Float
) -> Vec3