magnitudeVenus function
Visual magnitude of Venus. i is phase angle in radians.
Implementation
double magnitudeVenus(double r, double delta, double i) {
final id = toDeg(i);
return -4 + 5 * math.log(r * delta) / math.ln10 + (0.01322 + 0.0000004247 * id * id) * id;
}