declination
double declination(double t, double beta, double lambda);
double delta = declination(t, beta, lambda);
double delta = declination(t, 0, lambda);
double delta = declination(ut, 0, lambda);