计算方法插值法牛顿插值拉格朗日插值

------------------------------Р#include<iostream>Р#include<cmath>Рusing namespace std;Рint main()Р{Р double m=-2;Р int i, j, k , p, y=0;Р while(m<3)Р {Р double y=0; Р double f[100] ;Р double x[100] ;Р for(i=-5;i<6;i++)Р {Р x[i+5]=i;Р }Р //求出对于x的y值Р for(i=-5;i<6;i++)Р {Р f[i+5]=5.0/(1+i*i);Р }Р double l1[100] ;Р double l2[100] ;Р for(i=0;i<11;i++)Р {Р l1[i]=1;Р l2[i]=1;Р }Р Р for(k=0;k<11;k++)Р {Р for(p=0;p<11;p++)Р {Р if(p==k)Р continue;Р else{Р l1[k]=l1[k]*(m-x[p]);Р l2[k]=l2[k]*(x[k]-x[p]);Р }Р }Р }Р //多项式的每一项相加Р y = 0;Р for(j=0;j<11;j++)Р {Р y=y+l1[j]*f[j]/l2[j];Р }Р double h=5.0/(1+m*m); //精确值Р double r=fabs(y-h);Р cout<<"近似值y="<<y<<"\t精确值h="<<h<<"\t误差r="<<r<<endl;Р //for(int i=0;i<11;i++)Р //System.out.println("l1["+(i-5)+"]="+f[i]);Р m= m+0.2;Р Р }Р return 0;Р}Рy = 0Р输入xi,yi,及n,x,k=0