牛顿插值方法程序.doc

#include<stdio.h>

voidmain()
{
inti,j,k,m;floats[10][10];
printf("开始输入\n");
k=0;

do{
printf("请输入点的值");
scanf("%f",&s[k][0]);
printf("请输入点的函数值\n");
scanf("%f",&s[k][1]);
k=k+1;
for(j=2;j<k+1;j++)
for(i=1;i<k;i++)
s[i][j]=(s[i][j-1]-s[i-1][j-1])/(s[i][0]-s[i-j+1][0]);
printf("newton插值多项式\n");
for(i=0;i<k;i++)
printf("第%d次项系数为%f\n",i,s[i][i+1]);
printf("要继续输入吗？是，请输入1；否输入0\n");
scanf("%d",&m);

}while(m!=0&&k<100);

}

#include<stdio.h>
floatnewton(floats){inti,j;for(j=2;j<n+1;j++)for(i=1;i<n;i++){s[i][j]=(a[i][j-1]-a[i-1][j-1])/(a[i][0]-a[i-j+1][0]);
}returns;}
main(){inti,j,n;printf("有多少个点?");scanf("%d",&n);
floats[n][n+1]=0;printf("请输入点的值");for(i=0;i<n;i++)scanf("%f",&s[i][0]);printf("请输入各点的函数值");for(j=0;j<n;j++)scanf("%f",&s[j][1]);
floatnewton(s);
for(i=0;i<n;i++)for(j=i+1;j<n+1;i++)printf("newton插值多项式系数为");printf("%d次系数=%f",i,s[i][j]);
}

#include<stdio.h>
#definen3
voidmain()
{
inti,j;
floats[n][n+1]={0};
printf("请输入点的值\n");
for(i=0;i<n;i++)
scanf("%f",&s[i][0]);
printf("请输入各点的函数值\n");
for(j=0;j<n;j++)
scanf("%f",&s[j][1]);
for(j=2;j<n+1;j++)
for(i=1;i<n;i++)
s[i][j]=(s[i][j-1]-s[i-1][j-1])/(s[i][0]-s[i-j+1][0]);
printf("newton插值多项式系数为");
for(i=0;i<n;i++)
printf("系数=%f\n",s[i][i+1]);
}