Integrator vyuzivajici metodu Runge-Kutta 4.radu, zahrnujici Poyntinguv-Robertsonuv jev. -- Tento testovaci priklad souhlasi s Bertotti et al. (2003), odhadem tau_d. Casova skala pro cirkularizaci drahy je STEJNA jako pro pad na Slunce?! => e = 0 dosahnu az tesne pred padem! --- POZOR!!! ZDE JE CHYBA v n0 <- musim ho jeste vynasobit AU, aby se jednalo o jednotkovy vektor! Pak se ale objevi ony VELKE perturbace a, e!!! -> asi jsou vporadku, kdyz beta = a_rp/a_g dosahuje =~ 0.25, tak oskulacni elementy nejsou, co byvaly?! ??? Nemel bych je pocitat s JINOU hodnotou GM?! ANO! -- Co kdyz roste R a klesa beta, zmensuji se tyto perturbace?!