add orbitGetTimeOfFlight interface

src/orbits.h

@@ -163,3 +163,5 @@ orbitGetSpecificEnergy(double a, double mu);
 glm::vec2
 polarToRect(double true_anom, double r);
 
+double
+orbitGetTimeOfFlight(system_2body sys, double theta_begin, double theta_end);