牛顿插值公式：牛顿前插（C语言）

题目：

算法：

#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);
}