由于此地板涉嫌违规，因此它将折叠在系统中以隐藏此地板并显示此地板
ex_P2=ex_P1dydhh/(L1 )；
ex_P3=ex_P2dydhh*y/(L1 )/(L2 )；
dex_px=ex_p3-ex_p2；
ex _ px=ex _ P3-dex _ px * dex _ px/(ex _ P3-2 * ex _ p2ex _ P1 )；
}
return ex_px 1；
}
是双精度LNP (双精度)
{
if(x==1)返回0；
ELSEif(x2 )返回LNP )1/x；
ELSEif(x.1
{
双精度n=-1；
双精度a；
德奥
{
n=n-.6；
a=x/exp(n；
}
wile(A2|||A1 )；
返回LNP (a ) n；
}
double y=x-1，ln_p1=0，ln_p2=0，ln_p3=0，ln_p=0，ln_tmp=1，dln_px=1，TMN
int l；
for(L=1，tmp=1； (ln_px-ln_tmp(1e-10|| ) ln_px-ln_tmp () ) ) )。
{
ln_tmp=ln_px；
if(L==1) tmp=tmp/l；
else tmp=tmp/-l；
ln_p1=tmp；
ln_P2=ln_P1-1*tmp*y*L/(L1 )；
ln_P3=ln_P2dydhh*y*L/(L2 )；
dln_px=ln_p3-ln_p2；
ln _ px=ln _ P3-dln _ px * dln _ px/(ln _ P3-2 * ln _ p2ln _ P1 )；
}
return ln_px；
}
双精度(双精度，双精度y ) /计算x^y
{
if(x==0y！=0)返回0；
elseif(x==0y==0) return 1；
elseif(y0 )返回1/powp (x ) x，-y； //指数小于0时改为1/x^-y计算
Elseif(x0y-int ) y！=0)返回0； //x为负，y不是整数时出错，返回0
如果elseif(x0y-int(y )==0)/x为负且y为整数，则循环计算
{
双pow int=1；
int i；
for(I=1； i=y； I ) powint*=x；
返回pow int；
}
returnexpp(y*LNP ) x )；
}