Оптимизация функции двух переменных

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include<iostream>
#define N 2// количество переменных
double x_k[N];//неизвестные
double x_k_1[N];//значение последвательного приближения
double lambda_k;//шаг
 
double y(double *x, int n)
{
return (-2*pow(x[0],2)+3*x[1]*x[0] - 5*x[1]*x[1]+3*x[1]); //функция
}
 
double dy_dx0(double *x, int n) // первая частная производная по Х
{
return (-4*x[0]+3*x[1]-3);
}
 
double dy_dx1(double *x, int n) // первая частная производная по Y
{
return (3*x[0]-10*x[1]+3);
}
 
double dy2_dx0(double *x, int n)// 2-я частная производная по X
{
return (-4);
}
 
double dy2_dx1(double *x, int n)// 2-я частная производная по Y
{
return (-10);
}
 
int main(void)
{ 
setlocale(LC_CTYPE, "Russian");
lambda_k = 0.001;//шаг
x_k[0] = 0;//начальное
x_k[1] = 5;//приближение
while(1)//будет длится до конца интервала
{
x_k_1[0] = x_k[0]- lambda_k*dy_dx0(x_k, N) ;//последвательное
x_k_1[1] = x_k[1]-lambda_k*dy_dx1(x_k, N);//приближение
if(fabs(dy_dx0(x_k_1, N))>lambda_k && fabs(dy_dx1(x_k_1, N))<lambda_k) break;//условие при котором меняются значения меняется начальное и конечное значение промежутка
x_k[0] = x_k_1[0];//присвоение начальному значению
x_k[1] = x_k_1[1];//значение приближения
}
x_k[0] = x_k_1[0];
x_k[1] = x_k_1[1];
printf("\tГрадиентный метод:\n");
printf("\tМинимум найден на x1 =%.3lf, x2 =%.3lf, Y(X1,X2) =%3.3f\n", x_k[0], x_k[1], y(x_k, N));//Вывод минимальных точек и значения функции в этой точке
getch();
return 0;
}