弦截法是一种求方程根的基本方法,在计算机编程中常用。他的思路是这样的:任取两个数x1、x2,求得对应的函数值f(x1)、f(x2)。如果两函数值同号,则重新取数,直到这两个函数值异号为止。连接(x1,f(x1))与(x2,f(x2))这两点形成的直线与x轴相交于一点x,求得对应的f(x),判断其与f(x1)、f(x2)中的哪个值同号。如f(x)与f(x1)同号,则f(x)为新的f(x1)。将新的f(x1)与f(x2)连接,如此循环直到f(x)小于某个确定的精度为止。
编一程序,计算 x 3 + 2 ∗ x 2 + 5 ∗ x − 1 = 0 x^3 + 2*x^2+5*x-1=0 x3+2∗x2+5∗x−1=0在区间【-1,1】之间的根。
#include <iostream>
using namespace std;
int main()
{
double x, x1 = -1, x2 = 1, f1, f2, f, epsilon;
cout << "请输入精度:";
cin >> epsilon;
do {
f1 = x1 * x1 * x1 + 2 * x1 * x1 + 5 * x1 - 1;//计算f(x1)
f2 = x2 * x2 * x2 + 2 * x2 * x2 + 5 * x2 - 1;//计算f(x2)
x = (x1 * f2 - x2 * f1) / (f2 - f1);
f = x * x * x + 2 * x * x + 5 * x - 1;
if (f * f1 > 0) x1 = x;
else x2 = x;
} while (fabs(f) > epsilon);//误差在允许范围内
cout << "方程的根是:" << x << endl;
return 0;
}
运行结果: