经纬度计算他们距离的算法（！！！！！！！！！求助啊）

工作中要跟据GPS得到的经纬度计算两点间距离，在网上找了几个公式和程序想对比，下面这个好像比较专业，但里面有几个问题请高手指点啊，急啊，拜托了

问题1：里面用的e取值是多少啊，科学计数中的10，自然数2.718282，还是一个讲墨卡托航法纬度渐长率的的计算公式为：

double CChartCtrl::CalcltDstns(float fStarPtx, float fStarPty, float fEndPtx, float fEndPty)
{
//已知起始点坐标(fStartPtx, fStartPty)及到达点坐标(fEndPtx,fEndPty)
//计算航程dbDstns.
//起始点,到达点坐标:经纬度
//航程:海里(1852米)
//Created by zhl
//2002.7.3
//precision:0.0001 海里
//check param
double dbDir=CalcltDirct(fStarPtx,fStarPty,fEndPtx,fEndPty);

double delta_fy=fEndPtx-fStarPtx;
double delta_lnmg=fEndPty-fStarPty;

int mk=(int)fEndPtx*(int)fStarPtx;
double fy_m,dbDstns;//
if(mk>=0)
{//不跨赤道航行
fy_m=(fStarPtx+fEndPtx)/2;