非线性方程二分法求解

1. // binarySearch.cpp : Defines the entry point for the console application.  
2. //  
3.   
4. #include "stdafx.h"  
5.   
6. #include <CMATH>  
7. #include <IOSTREAM>  
8.   
9. using namespace std;  
10.   
11. // 非线性方程  
12. double func(double x)  
13. {  
14.     return 2*x*x*x-5*x-1;     
15. }  
16.   
17. // 二分法  
18. double BinarySearch(double a, double b, double err)  
19. {  
20.     err = fabs(err);  
21.     double c = (a + b) / 2;             // 取中间值  
22.   
23.     while(fabs(func(c)) > err && func(a-b) > err)  
24.     {  
25.         if ((func(c) * func(b)) < 0)  
26.         {  
27.             a = c;  
28.         }  
29.         if((func(c) * func(a)) < 0)  
30.         {  
31.             b = c;  
32.         }  
33.         c = (a + b) / 2;                // 二分法确定新区间  
34.     }  
35.     return c;  
36. }  
37.   
38. int main(int argc, char* argv[])  
39. {  
40.     double a = 1.0, b = 2.0;  
41.     double err = 1e-10;  
42.     double result;  
43.   
44.     result = BinarySearch(a ,b, err);  
45.   
46.     cout<<result<<endl;  
47.   
48.     return 0;  
49. }