JulianDay JD = new JulianDay(333, 1, 27, 12);
    Debug.Assert(JD.JulianDayNumber == 1842713.0);

    JulianDay JD1 = new JulianDay(1957, 10, 4) + 0.81;
    Debug.Assert(JD1.JulianDayNumber == 2436116.31);

    JulianDay JD = new JulianDay(1987, 04, 10, 0);

    double T = (JD.JulianDayNumber - 2451545.0) / 36525;

    Instant MeanSiderealDate = new Instant(6, 41, 50.54841);
    MeanSiderealDate += 8640184.812866 * T;
    MeanSiderealDate += 0.093104 * Math.Pow(T, 2);
    MeanSiderealDate += 0.0000062 * Math.Pow(T, 3);

    Debug.Assert(MeanSiderealDate.Hour == 13);
    Debug.Assert(MeanSiderealDate.Minute == 10);
    Debug.Assert(Math.Round(MeanSiderealDate.Second, 4) == 46.3668);

    Instant MeanSiderealInstant = new Instant(19, 21, 00);
    MeanSiderealInstant = MeanSiderealInstant * 1.00273790935;
    MeanSiderealInstant += MeanSiderealDate;
    Debug.Assert(MeanSiderealInstant.Hour == 8);
    Debug.Assert(MeanSiderealInstant.Minute == 34);
    Debug.Assert(Math.Round(MeanSiderealInstant.Second, 4) == 57.0896);

    IslamicCalendar MC = new IslamicCalendar(1421, 1, 1);
    JulianDay JD = MC.ToJulianDay();

    Debug.Assert(JD.Year == 2000);
    Debug.Assert(JD.Month == 4);
    Debug.Assert(JD.Day == 6);

    JulianDay JD1 = new JulianDay(1991, 8, 13);
    IslamicCalendar MC1 = IslamicCalendar.FromJulianDay(JD1);

    Debug.Assert(MC1.Year == 1412);
    Debug.Assert(MC1.Month == 2);
    Debug.Assert(MC1.Day == 2);

    double JDE = Equinox.Approximate(1962, Season.Summer);
    Debug.Assert(JDE == 2437837.39245);

    JulianDay JD = Events.Easter(179);
    Debug.Assert(JD.Month == 04);
    Debug.Assert(JD.Day == 12);

    JulianDay JD2 = Events.Easter(1991);
    Debug.Assert(JD2.Month == 03);
    Debug.Assert(JD2.Day == 31);

    JulianDay JD3 = Events.Pesach(2015);
    Debug.Assert(JD3.Month == 4);
    Debug.Assert(JD3.Day == 4);