拉格朗日插值法（C语言）

目录

一、函数

二、函数图像

三、程序结构

四、代码

五、运行结果

一、函数

                    𝑦=𝑥^3−5𝑥^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 n,i,j,k;
    double test_Y,P,X=0,Y=0;
    double x[10],y[10];
    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);
    for(k=0;k<n;k++)
    {
        P=1;
        for(j=0;j<n;j++)
        {
            if(j!=k)
                P=P*(X-x[j])/(x[k]-x[j]);
        }
        Y=Y+P*y[k];
    }
    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;
}

五、运行结果