fn gravity_accel(pos: [f64; 3], mu_val: f64) -> [f64; 3]
Compute gravitational acceleration at a position (km/s²).