>Алиса С
>Может через ...

Посмотри здесь и (под интерес) здесь


С уважением,
Владимир

>ВМоисеев (дополнительно), если хочешь посмотреть программы.

//-- Вычисление расстояния
//-------------------------------------------
SET @dist := 6371 * 2 * ASIN(SQRT(
POWER(SIN((src_lat - ABS(dst_lat)) * PI()/180 / 2), 2) +
COS(src_lat * PI()/180) *
COS(ABS(dst_lat) * PI()/180) *
POWER(SIN((src_lon - dst_lon) * PI()/180 / 2), 2)
));

//-------------------------------------
А если надо на эллипсоиде WGS84:

function TCore.GetDistance_Ellipse(Lon1, Lat1, Lon2, Lat2: double): double;
var
fdLambda : double;
fdPhi : double;
fPhimean : double;
fTemp : double;
fRho : double;
fNu : double;
fz : double;
fAlpha : double;
fR : double;
const
e2 = 0.0067394967423;
a = 6378137;
begin
Result:=0;

if (abs(Lon1-Lon2)<=0.0000014) and (abs(Lat1-Lat2)<=0.0000008) then begin
Exit;
end else begin
Lat1:=Lat1*Pi/180;
Lon1:=Lon1*Pi/180;
Lat2:=Lat2*Pi/180;
Lon2:=Lon2*Pi/180;

fdLambda:=(Lon1-Lon2);
fdPhi:=(Lat1-Lat2);
fPhimean:=(Lat1+Lat2)/2;
fTemp:=1-e2*Power(Sin(fPhimean),2);
fRho:=a*(1-e2)/Power(fTemp,1.5);
fNu:=a/sqrt(1-e2*Sin(fPhimean)*Sin(fPhimean));
fz:=2*ArcSin(sqrt(Power(Sin(fdPhi/2.0),2)+Cos(Lat2)*Cos(Lat1)*Power(Sin(fdLambda/2.0),2)));
fAlpha:=ArcSin(Cos(Lat2)*Sin(fdLambda)/Sin(fz));
fR:=fRho*fNu/(fRho*Power(Sin(fAlpha),2)+fNu*Power(Cos(fAlpha),2));
Result:=fz*fR;
end;
end;

С уважением,
Владимир.