【无标题】

#C++实现二分法代码

#include
#include
#include
using namespace std;

double gold(function<double(double)> func, double a, double b, double eps) {
double r = (sqrt(5.0) - 1.0) / 2.0; // 黄金比例
double c1, c2;
double fc1, fc2;
do {
c1 = b - r * (b - a);
c2 = a + r * (b - a);
fc1 = func(c1);
fc2 = func(c2);
if (fc1 < fc2) {
b = c2;
}
else {
a = c1;
}
} while (b - a > eps);
return (a + b) / 2.0; // 返回近似最优解
}

int main() {
auto f = [](double x) { return x * x - 2 * x + 1; }; //在[0, 2]区间内求最小值
double a = 0.0, b = 2.0, eps = 1e-7;
cin >> a >> b;
cout << endl;
double result = gold(f, a, b, eps);
cout<< result << std::endl;
return 0;
}