CLEAR >SET DECIMALS TO 5 > >m.lat1 = 36.12 >m.long1 = -86.67 >?"lat1",m.lat1,"long1",m.long1 > >m.lat2 = 33.94 >m.long2 = -118.40 >?"lat2",m.lat2,"long2",m.long2 > >m.DeltaSigma = MyDistance(m.lat1,m.long1,m.lat2,m.long2) > >?"delta-sigma",m.DeltaSigma > >m.Dist = m.DeltaSigma * 6372.795 >?m.Dist >RETURN > > >FUNCTION myDistance >LPARAMETERS m.lat1, m.long1, m.lat2, m.long2 >LOCAL m.rlat1, m.rlong1, m.rlat2, m.rlong2 > >m.rlat1 = DTOR(m.lat1) >m.rlong1 = DTOR(m.long1) >?"rlat1",m.rlat1,"rlong1",m.rlong1 > >m.rlat2 = DTOR(m.lat2) >m.rlong2 = DTOR(m.long2) >?"rlat2",m.rlat2,"rlong2",m.rlong2 > >m.DeltaLong = m.long2 - m.long1 > >m.rDeltaLong = DTOR(m.DeltaLong) > >m.DeltaSigma = ; > atan( ; > SQRT(; > (COS(m.rlat2) * SIN(m.rDeltaLong))^2 + ; > (COS(m.rlat1) * SIN(m.rlat2) - SIN(m.rlat1) * COS(m.rlat2) * COS(m.rDeltaLong))^2; > ) ; > / ; > (SIN(m.rlat1) * SIN(m.rlat2)) + (COS(m.rlat1) * COS(m.rlat2) * COS(m.rDeltaLong)) ; > ) > >return m.DeltaSigmaMike,
? GeoDist(36.12,-86.67,33.94,-118.40) Function GeoDist parameters lat1,lng1,lat2,lng2 * * Degrees to radian * lat1 = Dtor(lat1) lat2 = Dtor(lat2) lng1 = Dtor(lng1) lng2 = Dtor(lng2) return 6372.795 * acos(sin(lat1) * sin(lat2) +; cos(lat1) * cos(lat2) * cos(lng2 - lng1)) && kilometersCetin