void main()
{
cout<<" ***************插值法******************"<<endl;
cout<<"请选择您要完成的功能:"<<endl;
cout<<"<0>.退出 <1>.分段线性插值法 <2>.分段二次插值法 <3>.全区间上的拉格朗日插值法"<<endl;
int choice;
cin>>choice;
int number,a;
while(choice)
{
cout<<"请问您要输入几组数据?";
cin>>number;
cout<<"请输入您的数据:"<<endl;
float *px=new float[number];
float *py=new float[number];
for(a=0;a<number;a++)
cout<<"x"<<a+1<<" ";
cout<<endl;
for(a=0;a<number;a++)
cin>>px[a];
for(a=0;a<number;a++)
cout<<"y"<<a+1<<" ";
cout<<endl;
for(a=0;a<number;a++)
cin>>py[a];
switch(choice)
{
case 1: linearity(px,py,number);break;
case 2: squared(px,py,number);break;
case 3: lagrange(px,py,number);break;
default: cout<<"您的输入有误，请重试!"<<endl;
}
cout<<"继续吗?"<<endl
<<"<0>.退出 <1>.分段线性插值法 <2>.分段二次插值法 <3>.全区间上的拉格朗日插值法"<<endl;
cin>>choice;
delete px,py;
}
}
void linearity(float *px,float *py,int number)
{
float x,x0,x1,y0,y1,f;
int a;
cout<<"请输入您要计算的x的值:";
cin>>x;
float *pa=new float[number];
for(a=0;a<number;a++)
{
pa[a]=px[a]-x;
if(pa[a]<0)
pa[a]=-pa[a];
}
float *p=new float[number];
for(a=0;a<number;a++)
p[a]=pa[a];
float temp;
for(int pass=1;pass<number;pass++)
{
