T=(JD-2451545.0)/36525

theta_0=280.46061837+360.98564736629*(JD-2451545.0)+0.000387933*T^2-(T^3/38710000)