Задача по методам оптимизации

#include "stdafx.h"
#include <cmath>
 
bool sravn(double x[], double y[], double eps, int n)
{
    int i;
    int k = 0;
    bool s = true;
    for (i = 0; i<n; i++)
        if (fabs(x[i] - y[i]) <= eps)
        {
            k++;
        }
    if (k == n) s = false;
    return s;
}
 
double pr(double (*f) (double*) , double eps, double *x, int i, int n)
{
    double f1 = 0, f2 = 0;
    double *x1 = new double[n];
    double *x2 = new double[n];
    for (int j = 0; j < n; j++)
    {
        x1[j] = x2[j] = x[j];
    }
    double dx = 0.1;
    x1[i] += dx;
    x2[i] -= dx;
    f1 = (f(x1) - f(x2)) / (2 * dx);
    dx = dx / 10;
    for (int j = 0; j < n; j++)
    {
        x1[j] = x2[j] = x[j];
    }
    x1[i] += dx;
    x2[i] -= dx;
    f2 = (f(x1) - f(x2)) / (2 * dx);
    while (fabs(f2 - f1) > eps)
    {
        f1 = f2;
        dx /= 10;
        for (int j = 0; j < n; j++)
        {
            x1[j] = x2[j] = x[j];
        }
        x1[i] += dx;
        x2[i] -= dx;
 
        f2 = (f(x1) - f(x2)) / (2 * dx);
        if (f2 == 0)
            break;
    }
    return f1;
}
 
//double *proizv(double *x, int n, double eps, double(*fu) (double*))
//{
//  double *f = new double[n];
//  for (int i = 0; i < n; i++)
//  {
//      f[i] = pr(fu, eps, x, i, n);
//  }
//  return f;
//}
 
double rosenbr(double *x)
{
    return (1 - x[0])*(1 - x[0]) + 100 * (x[1] - x[0] * x[0])*(x[1] - x[0] * x[0]);
}
 
double fo(double *x)
{
    return pow(x[0], 3) + 2 * pow(x[1], 2) - 3 * x[0] - 4 * x[1];
}
double fo1(double *x)
{
    return x[0] * x[0]-4;
}
double fo2(double *x)
{
    return 3 * x[0] * x[0] + x[0] * x[1] + 2 * x[1] * x[1] - x[0] - 4 * x[1];
}
double fo3(double *x)
{
    return -x[0] * x[0] - 5 * x[1] * x[1] - 3 * x[2] * x[2] + x[0] * x[1] - 2 * x[0] * x[2] + 2 * x[1] * x[2] + 11 * x[0] + 2 * x[1] + 18 * x[2] + 10;
}
 
double pr11(double x)
{
    return 3 * x*x - 3;
}
 
double pr22(double x)
{
    return 4 * x - 4;
}
int main()
{
    double(*f) (double*);
    f = &fo;
    double ff;
    double eps = 0.00000001;
    int n=2;
    int k=0, p=0, iter=0;
    double a=0.5, b=0.2;
    double *T = new double[n];
    double *x0 = new double[n];
    for (int i = 0; i < n; i++)
    {
        scanf("%lf", &x0[i]);
    }
    double *x1 = new double[n];
    for (int i = 0; i < n; i++)
    {
        scanf("%lf", &x1[i]);
    }
    double *x = new double[n];
    double f0, f1;
    f0 = f(x0);
    f1 = f(x1);
    //printf("%lf %lf", f0, f1);
    if (f1 > f0)
    {
        T = x0;
        x0 = x1;
        x1 = T;
    }
 
    do
    {
        if (k == 1)
        {
            x0 = x1;
            x1 = x;
        }
        k = 0;
        iter++;
        for (int i = 0; i < n; i++)
        {
            x[i] = x1[i] - a*pr(f, 0.000001, x1, i, n)+b*(x1[i]-x0[i]);
        }
        //x[0]=x1[0]-a*pr11(x1[0])+b*(x1[0] - x0[0]);
        //x[1] = x1[1] - a*pr11(x1[1]) + b*(x1[1] - x0[1]);
        ff = f(x);
        if (ff > f1 && p == 0 && b > 0.1)
            b -= 0.1;
        else
        {
            if (b <= 0.1) p = 1;
            if (ff > f1 && p == 1 && a >= eps) a /= 2;
            if (ff<f1) k = 1;
        }
    } while (sravn(x1,x,eps,n) && a>eps);
    printf("%lf %lf ", x[0], x[1]);
    printf("%lf   iter %d  a %lf  b %lf", f(x), iter, a ,b);
 
 
 
 
    system("pause");
    return 0;
}