polcoe.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 31 行

#define NRANSI
#include "nrutil.h"

void polcoe(float x[], float y[], int n, float cof[])
{
	int k,j,i;
	float phi,ff,b,*s;

	s=vector(0,n);
	for (i=0;i<=n;i++) s[i]=cof[i]=0.0;
	s[n] = -x[0];
	for (i=1;i<=n;i++) {
		for (j=n-i;j<=n-1;j++)
			s[j] -= x[i]*s[j+1];
		s[n] -= x[i];
	}
	for (j=0;j<=n;j++) {
		phi=n+1;
		for (k=n;k>=1;k--)
			phi=k*s[k]+x[j]*phi;
		ff=y[j]/phi;
		b=1.0;
		for (k=n;k>=0;k--) {
			cof[k] += b*ff;
			b=s[k]+x[j]*b;
		}
	}
	free_vector(s,0,n);
}
#undef NRANSI