Определить минимум функции

struct Find_func_min_data {
    double(*func)(double, double);  // указатель на исследуемую функцию
    double x0, hx;
    double y0, hy;
    double xmin, ymin;   // здесь будет результат
    size_t nx, ny;
};
 
void find_func_min(Find_func_min_data& d) {
    double z, zmin = DBL_MAX, x, y;
    for (size_t i = 0; i < d.nx; ++i) {
        x = d.x0 + i * d.hx;
        for (size_t j = 0; j < d.ny; ++j) {
            y = d.y0 + j * d.hy;
            z = d.func(x, y);
            if (zmin > z) {
                zmin = z;
                d.xmin = x;
                d.ymin = y;
            }
        }
    }
}

#include <iostream>
#include <iomanip>
 
using namespace std;
 
struct Find_func_min_data {
    double(*func)(double, double);  // указатель на исследуемую функцию
    double x0, hx;
    double y0, hy;
    double xmin, ymin;   // здесь будет результат
    size_t nx, ny;
};
 
void find_func_min(Find_func_min_data& d) {
    double z, zmin = DBL_MAX, x, y;
    for (size_t j = 1; j <= d.ny; ++j) {
        y = d.y0 + j * d.hy;
        for (size_t i = 1; i <= d.nx; ++i) {
            x = d.x0 + i * d.hx;
            z = d.func(x, y);
            //cout << setw(6) << z << ' ';    // тестовая печать
            if (zmin > z) {
                zmin = z;
                d.xmin = x;
                d.ymin = y;
            }
        }
        //cout << endl;                       // тестовая печать
    }
}
 
double func1(double x, double y) {
    return 3 * x * x - 2 * y * y + 4 * x * y - 8 * x;
}
 
double func2(double x, double y) {
    return 6 * x * x + y * y - 4 * x - 6 * y + 1;
}
 
int main() {
    Find_func_min_data d1{
        func1,
        1.0, 0.2,
        2.0, 0.1,
        0.0, 0.0,
        10, 20
    };
    find_func_min(d1);
    cout << "min for function-1 is (" << d1.xmin << ", " << d1.ymin << ")\n";
 
    Find_func_min_data d2{
        func1,
        0.0, 0.1,
        -1.0, 0.1,
        0.0, 0.0,
        15, 10
    };
    find_func_min(d2);
    cout << "min for function-2 is (" << d2.xmin << ", " << d2.ymin << ")\n";
}