ratval.c

来自「Numerical Recipes Software 提供的算法子程序集」· C语言 代码 · 共 13 行

double ratval(x,cof,mm,kk)
double cof[],x;
int kk,mm;
{
	int j;
	double sumd,sumn;

	for (sumn=cof[mm],j=mm-1;j>=0;j--) sumn=sumn*x+cof[j];
	for (sumd=0.0,j=mm+kk;j>=mm+1;j--) sumd=(sumd+cof[j])*x;
	return sumn/(1.0+sumd);
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */