python编程作业9

[2019302857-陈春宇]

import math
def sphere_distance(p1,p2):
    a1=p1[0]
    b1=p1[1]
    a2=p2[0]
    b2=p2[1]
    r = 6371
    if (a1<0) or (b1<0) or (a1>90) or (b1>180) or (a2<0) or (b2<0) or (a2>90) or (b2>180):
        print("Parameter Error")
    else:
        c1 = a1 * math.pi / 180
        c2 = a2 * math.pi / 180
        d1 = b1 * math.pi / 180
        d2 = b2 * math.pi / 180
        m = 2 * r * math.asin(math.sqrt(math.sin((c2-c1)/2)**2+math.cos(c1)*math.cos(c2)*math.sin((d2-d1)/2)**2))
        m == round(m,2)
        print(m)
sphere_distance(p1=(45,45),p2=(45,90))