牛顿插值法（C语言）

目录

一、函数 

二、函数图像

三、代码

四、运行结果

---

一、函数 

                        y=x^3-5*x^2+6

二、函数图像

三、代码

#include <stdio.h>
#include <math.h>
/*y=x^3-5*x^2+6
x   y
1   2
2   -6
3   -12
4   -10
5   6
*/

int main()//二阶差商
{
    int i,j,n;
    double P=1,X,Y,x[10],y[10];
    double test_Y;
    printf("输入已知坐标个数n\n");
    scanf("%d",&n);
    printf("输入已知坐标（X,Y）\n");
    for(i=0;i<n;i++)
    {
        scanf("%lf%lf",&x[i],&y[i]);
    }
    printf("输入需要求的横坐标X\n");
    scanf("%lf",&X);
    Y=y[0];
    for(i=1;i<=n;i++)
    {

        for(j=1;j<=n-i+1;j++)
        {
            y[j-1]=(y[j]-y[j-1])/(x[j+i-1]-x[j-1]);
        }
        P=P*(X-x[i-1]);
        Y=Y+y[0]*P;
    }
    printf("X=%lf\tY=%lf\n",X,Y);
    test_Y=X*X*X-5*X*X+6;
    printf("真实值=%lf\t误差绝对值=%lf\n",test_Y,fabs(test_Y-Y));
    return 0;
}

四、运行结果