弦截法求方程的根

蒲奎民 — Wed, 23 Dec 2009 14:14:00 GMT

弦截法求方程的根是一种解方程得基本方法，在计算机编程中常用。他的思路是这样的：任取两个数，判断这两个数的函数值，如果函数值是同号，换两个数再试，直到两个数x1，x2对应的函数值为异号时为止，这时方程的解肯定在这两个数x1，x2之间。连接这两点所对应的函数值，连线与x轴的交点为新的x，若f(x)与f(x1)同号，则把x当作新的x1，将新的x1与x2连接，如此循环……如果f(x)与f(x1)异号，则把把x当作新的x2，将x1与新的x2连接，循环……

编程使用弦截法求方程f(x)=x^3-5x^2+16x-80=0的根的近似值(保留4位小数)，当fabs(f(x))<0.0001时就停止计算。5.0000

例如输入：

2 6

输出：5.0000

#include
#include
#include
using namespace std;
#define f(x) (x*x*x-5*x*x+16*x-80)
#define jd 0.000001
#define max 100
double aa;
bool gg(double x0, double x1)
{
    double xk, xk1, xk2;
    xk = x0;
    xk1 = x1;
    for (int i=0; i    {
        xk2 = xk1-f(xk1)*(xk1-xk)/(f(xk1)-f(xk));
        if (fabs(xk2-xk1) < jd)
        {
            aa = xk1;
            return 1;
        }
        else
        {
            xk = xk1;
            xk1 = xk2;
        }
    }
    return 0;
}
int main()
{
    double x0, x1;
    cin>>x0>>x1;
    if(gg(x0, x1))
    {
        cout<    }
    return 0;
}

蒲奎民 2009-12-23 22:14 发表评论

C++博客-真爱无限的C++博客-随笔分类-hhtc-C++

弦截法求方程的根