Für die Berechnung des Sonnen-auf und -untergangs gibt es allgemein Gültige Formeln. Die Richtigkeit wurde aber nicht überprüft. Da diese Formeln nur Näherungswerte ergeben, sind sie nicht als wirklich richtig anzusehen. Die Abweichung ist aber minimal.
Die Formeln sind unter anderem auf der Seite
http://lexikon.astroinfo.org/zeitgleichung/
entsprechend erklärt.
// WOZ - MOZ = -0.1752*sin(0.033430 * T + 0.5474) - 0.1340*sin(0.018234*T - 0.1939)
// Deklination = 0.40954*sin(0.0172*(T-79.35))
// Zeitdifferenz = 12*arccos((sin(h) - sin(B)*sin(Deklination)) / (cos(B)*cos(Deklination)))/pi;
// Aufgang Ortszeit = 12 - Zeitdifferenz - Zeitgleichung
// Untergang Ortszeit = 12 + Zeitdifferenz - Zeitgleichung