num eclipticLongitude(num M) { var C = equationOfCenter(M); var P = rad * 102.9372; // perihelion of the Earth return M + C + P + pi; }