double UTCForSolarAngleAux( int rise, double jd,
                            double latitudeRad,
                            double angleRad)
{
    double centuryTime = JulianCenturyFromJulianDay( jd);
    double equationOfTime = EquationOfTime(centuryTime);
    double sunDeclinationRad = SunDeclinationRad(centuryTime);
    double hourAngle = LocalHourAngleSunRad( latitudeRad, sunDeclinationRad,
                                             angleRad);
    if (!rise) { hourAngle = -hourAngle; }

    return 720.0 - (4.0 * RAD2DEG(hourAngle)) - equationOfTime; /* minutes */
}
bool 
TestEquationOfTime( )
{
    bool ok = true;
    cout << "Testing EquationOfTime" << endl;

    int prec = cout.precision( 12 );

    double jd = 2448908.5;
    cout << "jd = " << jd << endl;
    TESTCHECKFE( EquationOfTime( jd ).Days(),
                 TimeIncrement( 0, 13, 42.6 ).Days(), &ok, 1e-4 );

    cout.precision( prec );

    if ( ok )
        cout << "EquationOfTime PASSED." << endl << endl;
    else
        cout << "EquationOfTime FAILED." << endl << endl;
    return ok;
}