Идентификатор min не определен

a = min(x - h, x + h);

#include <iostream>
#include <math.h>
 
 
 
using namespace std;
 
/* -------------------------------------------
исследуемая функция
------------------------------------------- */
double function(double x) {
    return x*x - 2 * x; // тут можно написать свою функцию, например "return x*x*x;" для икс в кубе
}
 
/* -------------------------------------------
Свенн
------------------------------------------- */
void Svenn(double(*f)(double), double x, double& a, double& b, double h, int& k) {
    k = 0;
 
    if (f(x) < f(x + h))
        h *= -1;
 
    for (k = 0; f(x + h) < f(x); k++) {
        h *= 2;
        x += h;
    }
 
    a = min(x - h, x + h);
    b = max(x - h, x + h);
}
 
/* -------------------------------------------
золотое извращение 2
------------------------------------------- */
double GoldenSection2(double(*f)(double), double a, double b, int &k, double eps = 0.001) {
    double x1 = a + 0.38196 * (b - a);
    k = 0;
    while (abs(b - a) > eps) {
        double x2 = a + b - x1;
 
        if (x1 < x2)
            if (f(x1) < f(x2))
                b = x2;
            else {
                a = x1;
                x1 = x2;
            }
            if (x1 > x2)
                if (f(x1) < f(x2))
                    a = x2;
                else {
                    b = x1;
                    x1 = x2;
                }
                k++;
    }
    return (a + b) / 2;
}
 
/* -------------------------------------------
Пауэлл и компания
------------------------------------------- */
double Powell(double(*f)(double), double a, double b, int &k, double eps = 0.001) {
    k = 1;
    double c = b;
    b = (a + c) / 2;
 
    // фиг его знает что это, но, вроде, работает:
    double d = 0.5*(f(a)*(b*b - c*c) + f(b)*(c*c - a*a) + f(c)*(a*a - b*b)) / (f(a)*(b - c) + f(b)*(c - a) + f(c)*(a - b));
 
    while (abs((d - b) / b) > eps || abs((f(d) - f(b)) / f(b)) > eps) {
        if (k > 1)
            d = ((a + b) / 2) + (0.5*(f(a) - f(b))*(f(a) - f(b))*(b - c)*(c - a)) / (f(a)*(b - c) + f(b)*(c - a) + f(c)*(a - b));
 
        if (f(b) > f(d)) {
            c = b;
            b = d;
        }
        else {
            a = b;
            b = d;
        }
        k++;
    }
    return (b + d) / 2;
}
 
int main(void) {
 
    int k0, k1, k2;    // итерации для каждого из методов
    double a, b,       // здесь будут начальные А и Б после свенна
        x_m, min,      // минимум после золотого и минимум после пауэлла.
        h = 0.01,    // ----------------- тут задаёте шаг для свенна
        eps = 0.0001,  // ----------------- тут задаёте точность эпсилон
        x0 = 1.0;     // ----------------- тут начальную точку x0
 
    system("cls");
 
    Svenn(function, x0, a, b, 0.001, k0);
    x_m = GoldenSection2(function, a, b, k1, eps);
    // поскольку пауэлл самый тяжёлый, в качестве А и Б мы берём минимум из предыдущего
    // метода плюс минус эпсилон.
    min = Powell(function, x_m - eps, x_m + eps, k2, eps);
    // ну и типа вывод: 
    // пауэлл более жёсткий, поэтому нужно для него более точно указывать промежуток, все
    // дела (короче, по подобию нашего вывода отсюда: [url]http://7bit.net.ru/MO/462[/url])
 
    cout << "Epsilon: " << eps << endl;
    cout << "x0:      " << x0 << endl;
    cout << "------------------------" << endl;
    cout << "GS-2:    " << x_m << " / Iterations: " << k1 << endl;
    cout << "Powell:  " << min << " / Iterations: " << k2 << endl;
    cout << endl;
 
    system("pause");
 
    return 0;
}

#include <algorithm>