PRINT "Distance = " ; FNhaversine(36.12, -86.67, 33.94, -118.4) " km"
      END

      DEF FNhaversine(n1, e1, n2, e2)
      LOCAL d() : DIM d(2)
      d() = COSRAD(e1-e2) * COSRAD(n1) - COSRAD(n2), \
      \     SINRAD(e1-e2) * COSRAD(n1), \
      \     SINRAD(n1) - SINRAD(n2)
      = ASN(MOD(d()) / 2) * 6372.8 * 2