四牛顿插值多项式

四牛顿插值多项式牛顿插值多项式1、实验目的:按下列数据作五次插值，并求x1=0.46,x2=0.55,x3=0.60时的函数近似值Xi0.300.420.500.580.660.72Yi1.044031.084621.118031.156031.198171.232232、实验步骤:添加库函数#include"stdafx.h"#include"conio.h"3、代码:#include"stdafx.h"#include"conio.h"#defineN5voiddiffer(floatx[],floaty[],intn){float*f=newfloat[n+1];intk,i;for(k=1;k<=n;k++){f[0]=y[k];for(i=0;i<k;i++)f[i+1]=(f[i]-y[i])/(x[k]-x[i]);y[k]=f[k];}}voidmain(){inti,m;floatx[N+1]={0.30,0.42,0.50,0.58,0.66,0.72};floaty[N+1]={1.04403,1.08462,1.11803,1.15603,1.19817,1.23223};floatb,varx[3]={0.46,0.55,0.60};differ(x,y,N);b=y[N];for(m=0;m<=2;m++){for(i=N-1;i>=0;i--)b=b*(varx[m]-x[i])+y[i];printf("f(%f)=%f\n",varx[m],b);}}4、运行结果:5、:熟悉了牛顿插值法的算法