《C++程序设计(谭浩强)》代码,均按章节分类
源代码在线查看: 例4.9.txt
例4.9 用弦截法求方程f(x)=x3次方-5x2次方+16x-80=0的根。
#include
#include
#include
using namespace std;
double f(double); //函数声明
double xpoint(double, double); //函数声明
double root(double, double); //函数声明
int main( )
{ double x1,x2,f1,f2,x;
do
{cout cin>>x1>>x2;
f1=f(x1);
f2=f(x2);
} while(f1*f2>=0);
x=root(x1,x2);
cout //指定输出7位小数
cout return 0;
}
double f(double x) //定义f函数,以实现f(x)
{double y;
y=x*x*x-5*x*x+16*x-80;
return y;
}
double xpoint(double x1, double x2) //定义xpoint函数,求出弦与x轴交点
{double y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1)); //在xpoint函数中调用f函数
return y;
}
double root(double x1, double x2) //定义root函数,求近似根
{double x,y,y1;
y1=f(x1);
do
{x=xpoint(x1,x2); //在root函数中调用xpoint函数
y=f(x); //在root函数中调用f函数
if (y*y1>0)
{y1=y;
x1=x;
}
else
x2=x;
}while(fabs(y)>=0.00001);
return x;
}