pccheb.c

来自「适合大型数值计算代码 现在网络上已经找不到了 购买需要20$」· C语言 代码 · 共 21 行

void pccheb(float d[], float c[], int n)
{
	int j,jm,jp,k;
	float fac,pow;

	pow=1.0;
	c[0]=2.0*d[0];
	for (k=1;k<n;k++) {
		c[k]=0.0;
		fac=d[k]/pow;
		jm=k;
		jp=1;
		for (j=k;j>=0;j-=2,jm--,jp++) {
			c[j] += fac;
			fac *= ((float)jm)/((float)jp);
		}
		pow += pow;
	}
}