やはり、RのコードもWolframのコードもないクズ書き込みのみ。
Fラン確定！

Clear[{R,f,a,b,g}]

(* 球の半径 *)
R = 20; (* 半径を任意の値に設定 *)

f[a_,b_]:=(
{lat1,lon1}={a,0} Degree;
{lat2,lon2}={b,90} Degree;
cosc=Sin[lat1]Sin[lat2]+Cos[lat1]Cos[lat2]Cos[lon1-lon2];
R ArcCos[cosc]
)
f[a,b]

(* 点の座標 *)
latLonToCartesian[R_, φ_, λ_] := R * {
Cos[φ Degree] * Cos[λ Degree],
Cos[φ Degree] * Sin[λ Degree],
Sin[φ Degree]
};

g[a_,b_]:=(
(* 点 pC, pA, pB を設定 *)
pC = latLonToCartesian[R, 90, 0] ;(* 点 P (90°, 経度0°) *)
pA = latLonToCartesian[R, a, 0] ; (* 点 A (緯度a°, 経度0°) *)
pB = latLonToCartesian[R, b , 90] ; (* 点 B (緯度0°, 経度90°) *)

(* 大円の法線ベクトル *)
n = Cross[pA, pB];

(* 点 P から大円への最短距離（球面上） *)
distance = ArcSin[Abs[Dot[n, pC]] / Norm[n]] * R
)
g[a,b]