牛顿插值法（C语言）

最新推荐文章于 2022-12-25 21:36:09 发布

原创最新推荐文章于 2022-12-25 21:36:09 发布 · 936 阅读

13 ·

CC 4.0 BY-SA版权

文章标签：

#c语言

C算法专栏收录该内容

6 篇文章

订阅专栏

这篇博客介绍了如何用C++编程来表示和描绘一个三次多项式函数y=x^3-5*x^2+6的图像，并通过二阶差商法计算给定点的函数值，同时展示真实值与计算值之间的误差。代码实现中涉及到了输入处理、数学运算和误差分析。

一、函数

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