该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
ex_p2=ex_p1+tmp*y/(l+1);
ex_p3=ex_p2+tmp*y*y/(l+1)/(l+2);
dex_px=ex_p3-ex_p2;
ex_px=ex_p3-dex_px*dex_px/(ex_p3-2*ex_p2+ex_p1);
}
return ex_px+1;
}
double lnP(double x)
{
if(x==1) return 0;
else if(x>2) return -lnP(1/x);
else if(x<.1>
{
double n=-1;
double a;
do
{
n=n-.6;
a=x/exp(n);
}
while(a>2 || a<1);
return lnP(a)+n;
}
double y=x-1,ln_p1=0,ln_p2=0,ln_p3=0,ln_p=0,ln_px=0,ln_tmp=1,dln_px=1,tmp;
int l;
for(l=1,tmp=1;(ln_px-ln_tmp)>1e-10 || (ln_px-ln_tmp)
{
ln_tmp=ln_px;
tmp*=y;
if(l==1) tmp