-- перевод в доли градуса по формуле: 
   -- Длина или Широта = Градусы + Минуты/60 + Секунды/3600, где Градусы,Минуты,Секунды: целочисленные

   -- код тестировался на PL/SQL для Oracle 9i

   FUNCTION DistanceInMeter(LaengeA in float, -- долгота A  в градусах
                            BreiteA in float, -- широта A  в градусах
                            LaengeB in float, -- долгота B  в градусах
                            BreiteB in float) -- широта B  в градусах
   RETURN INTEGER 
      Result   integer;

      fPhimean float;
      fdLambda float;
      fdPhi    float;
      fAlpha   float;
      fRho     float;
      fNu      float;
      fR       float;
      fz       float;
      fTemp    float;

      D2R      float := 0.017453;
      a        float := 6378137.0;
      e2       float := 0.006739496742337;

   BEGIN
      if not((LaengeA is NULL) or (BreiteA is NULL) or (LaengeB is NULL) or (BreiteB is NULL)) then
         fdLambda := (LaengeA - LaengeB) * D2R;
         fdPhi := (BreiteA - BreiteB) * D2R;
         fPhimean := ((BreiteA + BreiteB) / 2.0) * D2R;

         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 :=
            Sqrt(Power(Sin(fdPhi / 2.0), 2) + Cos(BreiteB * D2R) * Cos(BreiteA * D2R) *
               Power(Sin(fdLambda / 2.0), 2));

         fz := 2 * ASin(fz);

         if fz=0 then
            Result:=0;
         else
            fAlpha := Cos(BreiteB * D2R) * Sin(fdLambda) * 1 / Sin(fz);
            fAlpha := ASin(fAlpha);

            fR := (fRho * fNu) / ((fRho * Power(Sin(fAlpha), 2)) + (fNu *
               Power(Cos(fAlpha), 2)));

            Result := (fz * fR);
         end if;
      else
         Result:=NULL;
      end if;

      return(Result);
   END DistanceInMeter;