public static double CalcSunEqofCtr(double GeomMeanAnomSun, double JulianCentury)
{
double SunEqofCtr = Math.Sin(Calc.rad(GeomMeanAnomSun)) * (1.914602 - JulianCentury
* (0.004817 + 0.000014 * JulianCentury)) + Math.Sin(Calc.rad(2 * GeomMeanAnomSun))
* (0.019993 - 0.000101 * JulianCentury) + Math.Sin(Calc.rad(3 + GeomMeanAnomSun)) * 0.000289;
return(SunEqofCtr);
}
public static double CalcObliqCorr(double MeanObliqEcliptic, double JulianCentury)
{
double ObliqCorr = MeanObliqEcliptic + 0.00256 * Math.Cos(Calc.rad(125.04 - 1934.136 * JulianCentury));
return(ObliqCorr);
}
