[求助]卡方分布函数？？

noahhu

以下是我参照数理统计实用算法编写的卡方分布函数，这个程序在df为偶数时正确，但df为奇数时计算结果不正确，请各位高手帮忙，看看什么原因。另外，论坛中<font size="-1"><font color="#cc0033">Handbook Of Mathematical Functions找不到，如果哪位能下载请帮我下载一下， 我的邮箱是：dreams_wing@163.com 谢谢了！</font></font><br>double CDFChi(double X2, int df){<br>    int n=df, IAI;<br>    double F, u, pp;<br>    if (X2==0) return 0.0;<br>    double PIS=sqrt(PI());<br>    double x=X2/2;<br>    double CHS=sqrt(x);<br><br>    if (n%2==0) {<br>        u=x*exp(-x);<br>        F=1-exp(-x);<br>        IAI=2;<br>    }<br>    else {<br>        u=sqrt(x)*exp(-x)/PIS;<br>        pp=CDFNor(CHS);<br>        F=2*(pp-0.5);<br>        IAI=1;<br>    };<br><br>    if (IAI==n) return F;<br>    double n2=n-2;<br>    for (double i=IAI; i&lt;=n2; i+=2) {<br>        F=F-2*u/i;<br>        u=X2*u/i;<br>    };<br>    return F;<br>};

johhan

C语言编的,俺们看不懂!!!!!!!!!!!