import?math
print（‘大地緯度按度分秒輸入‘）
d1?=?float（input（‘d1?=?‘））
f1?=?float（input（‘f1?=?‘））
m1?=?float（input（‘m1?=?‘））
print（‘大地經度按度分秒輸入‘）
d2?=?float（input（‘d2?=?‘））
f2?=?float（input（‘f2?=?‘））
m2?=?float（input（‘m2?=?‘））
H?=?float（input（‘H?=?‘））
a?=?6378137
b?=?6356752.3142
e?=?pow（a*a?-?b*b1/2）/a
e1?=?pow（a*a?