#!/usr/bin/env python
import ephem
import math
pi=math.pi
r2d=180/pi
R=6371
def gleason(lon,lat):
r=(pi/2-lat)*R
x=r*math.cos(lon)
y=r*math.sin(lon)
return (x,y)
melbourne = ephem.city("Melbourne")
newyork = ephem.city("New York")
santiago= ephem.city("Santiago")
dis=ephem.separation(melbourne,newyork)*R
print("Distancia ortodrómica Melbourne-Nueva York:",dis)
dis=ephem.separation(melbourne,santiago)*R
print("Distancia ortodrómica Melbourne-Santiago:",dis)
dis=ephem.separation(santiago,newyork)*R
print("Distancia ortodrómica Santiago-Nueva York:",dis)
(x1,y1)=gleason(melbourne.lon,melbourne.lat)
(x2,y2)=gleason(newyork.lon,newyork.lat)
(x3,y3)=gleason(santiago.lon,santiago.lat)
dis=math.sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))
print("Distancia en Gleason Melbourne-Nueva York:",dis)
dis=math.sqrt((x1-x3)*(x1-x3)+(y1-y3)*(y1-y3))
print("Distancia en Gleason Melbourne-Santiago:",dis)
dis=math.sqrt((x2-x3)*(x2-x3)+(y2-y3)*(y2-y3))
print("Distancia en Gleason Santiago-Nueva York:",dis)