pub fn velocity( parallax: f64, right_ascension: f64, declination: f64, proper_motion_ra: f64, proper_motion_dec: f64, radial_velocity: f64) -> DVec3