//預(yù)先定義插值節(jié)點的個數(shù)為1000個，根據(jù)控制臺輸入的個數(shù)num從而確定插值節(jié)點的個數(shù)const int N=1000;// 拉格朗日插值算法float lglr(float x[], float y[],int n,float t){float yResult=0.0;//LValue[N]存放的是每次求解的插值基函數(shù)的通項float LValue[N];//循環(huán)變量k,mint k,m;//插值基函數(shù)中的上下累乘temp1,temp2float temp1,temp2;for(k=0;k { temp1=1.0; temp2=1.0;for(m=0;m {if(m==k) {continue; } temp1 *= (t-x[m]); temp2 *= (x[k]-x[m]); } LValue[k]=temp1/temp2; }for(int i=0;i { yResult += y[i]*LValue[i]; }return yResult;}