#include<iostream>
#include<cmath>
using namespace std;
double fx(double x)
int main()
{
const double eps=1.e-5;
double L,R,M;
double fL,fM;
int Count=0;
cout<<"单调区间左端点L=";cin>>L;
cout<<"单调区间右端点R=";cin>>R;
cout<<"f(L)="<<fx(L)<<",f(R)"<<fx(R)<<endl;
if(fx(L)*fx(R)>0)
{
cout<<"非有解区间:\n";system("pause");return 0;
}
while(fabs(R-L)>eps)
{
cout<<++Count<<"次:";
M=(L+R)/2;
fL=fx(L);fM=fx(M);
if(fL*fM<0)R=M;
else if(fL*fM>0)L=M;
else break;
cout<<L<<","<<R<<endl;
}
cout<<"方程的解:"<<"x="<<M<<"\n函数值:"<<fx(M)<<endl;
system("pause");return 0;
}
#include<iostream>
using namespace std;
const double x[7]={1.0,1.5,2.3,3.5,5.0,6.2,7.0}
const double fx[7]={0.8,1.0,2.0,2.5,1.6,1.8,1.4}
double Linear(double xx,int L)
{
double y;
y=fx[L]*(xx-x[L+1])/(x[L]-x[L+1])+fx[L+1]*(xx-x[L])/(x[L+1]-x[L]);
return y;
}
double Parabola(double xx,int L)
{
double y;
y=fx[L]*(xx-x[L+1])*(xx-x[L+2])/(x[L]-x[L+1])*(x[L]-x[L+2])+
fx[L+1]*(xx-x[L])*(xx-x[L+2])/((x[L+1]-x[L])*(x[L+1]-x[L+2]))+
fx[L+2]*(xx-x[L])*(xx-x[L+1])/((x[L+2]-x[L])*(x[L+2]-x[L+1]));
return y;
}
int main()
{
double xx=0;
cout<<" 输入-1结束"<<endl;
while(xx!=-1)
{
cout<<" 输入插值点:";
cin>>xx;
if(xx==-1)
{
cout<<"结束\n";
system("pause")
;
return 0;
}
if(xx<1||xx>7){
cout<<"插值点越界!\n";continue;
}
if(xx>x[4])cout<<"f(x)="<<Parabola(xx,4)<<endl;
else if(xx>x[3])cout<<"f(x)="<<Linear(xx,3)<<endl;
else if(xx>x[2])cout<<"f(x)="<<Linear(xx,2)<<endl;
else cout<<"f(x)="<<Parabola(xx,0)<<endl;
}
system("pause");
return 0;
}