23/07/22 05:46:22.35 54NVUurG0.net
fn <- function(
ai=45, # 北緯°
ak=145, # 東経°
bi=45, # 北緯°(南緯はー°)
bk=-125){ # 東経°(西経はー°)
R=6378 # 地球の半径km
u=pi/180
A=R*c(cos(ai*u)*cos(ak*u),cos(ai*u)*sin(ak*u),sin(ai*u))
B=R*c(cos(bi*u)*cos(bk*u),cos(bi*u)*sin(bk*u),sin(bi*u))
AB=sqrt(sum((A-B)^2))
asin(AB/(2*R))*2*R
}
fn()