% ~~~~~ function x = kepler_U(dt, ro, vro, a)% ~~~~~%{ This function uses Newton's method to solve the universal Kepler equation for the universal anomaly. mu - gravitational parameter (km^3/s^2) x - ...