Метод градиентного спуска

#include <iostream>
#include <math.h>
#include <conio.h>
#include <iomanip>
#include <fstream>
#include <stdlib.h>
using namespace std;
 
double f(double x, double y) // функцию на свою 
{
    return  y*y + cos(x + 1.23*y);
};
 
double gradX(double x, double y)
{
    return -sin(x + 1.23*y);
};
 
double gradY(double x, double y)
{
    return  2*y - 1.23*sin(x + 1.23*y);
};
 
double norm(double x, double y) //наша норма
{
    return sqrt((-sin(x + 1.23*y))*(-sin(x + 1.23*y)) + (2*y - 1.23*sin(x + 1.23*y))*(2*y - 1.23*sin(x + 1.23*y)));
};
 
 
void output(double x, double y)
    {
        cout << setw(4) << fixed << setprecision(5) << x;
        cout << setw(14) << fixed << setprecision(5) << y;
        cout << setw(14) << fixed << setprecision(5) << f(x,y) << endl;
    };
 
int main()
{
    double x0 = 7;
    double y0 = 8;
    double h = 0.0001;
   // double M = 6;
    double x = x0;
    double y = y0;
    output(x,y);
    for(int k = 0; k < 4; ++k)
    {   
        double gradx = gradX(x,y);
        double grady = gradY(x,y);
        double Norm = norm(x,y);
        x = x - h*gradx/Norm;
        y = y - h*grady/Norm;
        
        output(x,y);
    };
    _getch();
};