Have you tried that Math::Trig function great_circle_distance? Even if your task was to demonstrate code to calculate great circle distance, it would let you compare your results to those from reliable code.
Information about American English usage here and here. Floating point issues? Please read this before posting. — emc