Apsidal evolution under Poynting–Robertson drag (source code)

= Apsidal evolution under Poynting–Robertson drag