题目:
算法:
#include<stdio.h>
#define N 3
void difference(float y[],float f[4][4],int n)
{
int k,i;
f[0][0]=y[0];
f[1][0]=y[1];
f[2][0]=y[2];
f[3][0]=y[3];
for(k=1;k<=n;k++)
for(i=0;i<=(N-k);i++)
f[i][k]=f[i+1][k-1]-f[i][k-1];
return;
}
void main()
{
int i,k=1;
float a,b=1,m=21.4,t=1.4,f[4][4]={0};
float x[5]={20,21,22,23,24};
float y[5]={1.30103,1.32222,1.34342,1.36173,1.38021};
difference(y,f,N);
a=f[0][0];
for(i=1;i<=N;i++)
{
k=k*i;
b=b*(t-i+1);
a+=b*f[0][i]/k;
}
printf("差分表\nx(k)\n");
for(i=0;i<=4;i++)
printf("%-20f",x[i]);
printf("\ny(k)\n");
for(i=0;i<=4;i++)
printf("%-20f",y[i]);
for(k=1;k<=3;k++)
{
printf("\nF(%d)\n",k);
for(i=0;i<=(3-k);i++)
{
printf("%-20f",f[i][k]);
}
}
printf("\nf(%f)=%-20f\n",m,a);
}