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.DeltaSigma