Question

C语言编写拉格朗日插值函数，题目，代码如下。哪里错了

设y=1/x,节点 x0=2, x1=2.5, x2=4，求函数y=1/x的插值多项式L（x）。提示：L(3)=0.325。#include"stdio.h"
#include"math.h"
float Lagrange(float array[],int n,float x)
{
float a=1.0,b=1.0,l[3];
int i,j;
for(i=0;i<=n-1;i++)
{
for(j=0;j<=n-1;j++)
{ if(i==j) continue;
a=a*(x-array[j]);
b=b*(array[i]-array[j]);
}
l[i]=a/b;
return l[i];
}
}
main()
{int k,m;
float t,L=0.0,y[3],l[3],x[3]={2.0,2.5,4.0};
printf("please input t=");
scanf("%f",&t);
for(k=0;k<=2;k++)
{ y[k]=1/x[k];}
m=sizeof(x)/sizeof(x[0]);
for(k=0;k<=2;k++)
{l[k]= Lagrange(x,m,t);
L=L+l[k]*y[k];
}
printf("L=%f\n",L);
}

Answer 1

我自己写的，给你个例子：
#include<iostream.h>
void main()
{
char q;
do{
double X;
cout<<"请输入 X=";
cin>>X;
double l=1,L=0;
double x[5]={-3.0,-1.0,1.0,2.0,3.0};
double y[5]={1.0,1.5,2.0,2.0,1.0};
for(int i=0;i<5;i++)
{
for(int j=0;j<5;j++)
{
if(i==j)continue;
l*=(X-x[j])/(x[i]-x[j]);
}
L+=l*y[i];
l=1;
}
cout<<"x="<<X<<"时函数的近似值 L="<<L<<endl;
cout<<"按任意键继续，按q退出";
cin>>q;
}while(q!='q');
}