To help me with calculating the visual magnitude of the International Space Station I need to be able to calculate the phase angle.
Can anyone help me calculate t
Since angle_b is practically zero (max 0.00016 degrees at perihelion and ISS directly overhead), you could just do angle_a = math.pi - angle_c.
angle_b
angle_a = math.pi - angle_c
This has about the same accuracy, since you already have made some minor errors with distance a.
a