1. 1.
    +2
    d12 = 2.asin( √(sin²(Δφ/2) + cos(φ1).cos(φ2).sin²(Δλ/2)) )
    φ1 = acos( sin(φ2) − sin(φ1).cos(d12) / sin(d12).cos(φ1) )
    φ2 = acos( sin(φ1) − sin(φ2).cos(d12) / sin(d12).cos(φ2) )

    if sin(λ2−λ1) > 0
    θ12 = φ1, θ21 = 2.π − φ2
    else
    θ12 = 2.π − φ1, θ21 = φ2

    α1 = (θ1 − θ12 + π) % 2.π − π
    α2 = (θ21 − θ2 + π) % 2.π − π

    α3 = acos( −cos(α1).cos(α2) + sin(α1).sin(α2).cos(d12) )
    d13 = atan2( sin(d12).sin(α1).sin(α2), cos(α2)+cos(α1).cos(α3) )
    φ3 = asin( sin(φ1).cos(d13) + cos(φ1).sin(d13).cos(θ1) )
    Δλ13 = atan2( sin(θ1).sin(d13).cos(φ1), cos(d13)−sin(φ1).sin(φ3) )
    λ3 = (λ1+Δλ13+π) % 2.π − π
    where

    φ1, λ1, θ1 : 1st point & bearing
    φ2, λ2, θ2 : 2nd point & bearing
    φ3, λ3 : intersection point

    % = mod
    ···
   tümünü göster