Java tutorial
//package com.java2s; public class Main { /** * Calculate the equation of center for the sun. This value is a correction * to add to the geometric mean longitude in order to get the "true" longitude * of the sun. * * @param t number of Julian centuries since J2000. * @return Equation of center in degrees. */ private static double sunEquationOfCenter(final double t) { final double m = Math.toRadians(sunGeometricMeanAnomaly(t)); return Math.sin(1 * m) * (1.914602 - t * (0.004817 + 0.000014 * t)) + Math.sin(2 * m) * (0.019993 - t * (0.000101)) + Math.sin(3 * m) * (0.000289); } /** * Calculate the Geometric Mean Anomaly of the Sun. * * @param t number of Julian centuries since J2000. * @return Geometric Mean Anomaly of the Sun in degrees. */ private static double sunGeometricMeanAnomaly(final double t) { return 357.52911 + t * (35999.05029 - 0.0001537 * t); } }