Main loop optimization
Main loop optimization
'get_ang_momentum_2' subroutine for angular momentum conservation computations
energy, angular momentum, and d_energy & d_ang_momentum
orbit with eccentricity