13 lines
427 B
Text
13 lines
427 B
Text
DEF Haversine (lat1, long1, lat2, long2)
|
|
OPTION ANGLE RADIANS
|
|
LET R = 6372.8 !radio terrestre en km.
|
|
LET dLat = RAD(lat2-lat1)
|
|
LET dLong = RAD(long2-long1)
|
|
LET lat1 = RAD(lat1)
|
|
LET lat2 = RAD(lat2)
|
|
LET Haversine = R *2 * ASIN(SQR(SIN(dLat/2)^2 + SIN(dLong/2)^2 *COS(lat1) * COS(lat2)))
|
|
END DEF
|
|
PRINT
|
|
|
|
PRINT "Distancia de Haversine:"; Haversine(36.12, -86.67, 33.94, -118.4); "km"
|
|
END
|