用拉格郎日插值多项式求函数近似值.doc

文档介绍：用拉格郎日插值多项式求函数近似值
#include<>
#include<>
#include<>
#define max 100
float lagrange(float *x,float *y,float xx,int n)
{
int i,j;
float *a,yy=;
a=(float *)malloc(n*sizeof(float));
for(i=0;i<=n-1;i++)
{
a[i]=y[i];
for(j=0;j<=n-1;j++)
if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]);
yy+=a[i];
}
free(a);
return yy;
}
int main()
{
while(1)
{
float x[max];
float y[max];
float xx,yy;
printf("请输入x的值\n");
scanf("%f",&xx);
int m,i;
printf("请输入个数:\n");
scanf("%d",&m);
printf("请输入x(i):\n");
for(i=0;i<m;i++)
scanf("%f",&x[i]);
printf("请输入y(i):\n");
for(i=0;i<m;i++)
scanf("%f",&y[i]);
yy=lagrange(x,y,xx,m);
//clrscr();
printf("\n结果:\n");
printf("x=%f,y=%f\n",xx,yy);
getch();
}
}
牛顿差值多项式求函数
#include <>
#include<>
#include<>
#define N 4
void Difference(float *x,float *y,int n)
{

float * f;
int k,i;
f=(float *)malloc(n*sizeof(float));
for(k=1;k<=n;k++)
{

f[0]=y[k];
for(i=0;i<k;i++)
f[i+