Несовместимость параметра double

#include <iostream>
#include<stdio.h>
#include<math.h>
#include <iomanip>
#include<stdlib.h>
using namespace std;
 
double MyFxDelegate(int nNumVars, double fX[], double fParam[]);
int nNumVars = 2;
 
double MyFx(int N, double X[], double  fParam[])
{
    return 10 + pow(X[0] - 2, 2) + pow(X[1] + 5, 2);
}
double MyFxEx(int nNumVars, double fX[], double fParam[], double fDeltaX[], double fLambda)
{
 
    int i;
    double fXX[2];
 
    for (i = 0; i < nNumVars; i++)
    {
        fXX[i] = fX[i] + fLambda * fDeltaX[i];
    }
 
    return MyFxDelegate(nNumVars, fXX, fParam);
}
void GetGradients(int nNumVars, double fX[], double fParam[], double fDeriv[], double fDerivNorm)
{
 
    int i;
    double fXX, H, Fp, Fm;
 
    fDerivNorm = 0;
    for (i = 0; i < nNumVars; i++)
    {
        fXX = fX[i];
        H = 0.01 *(1 + abs(fXX));
        fX[i] = fXX + H;
        Fp = MyFxDelegate(nNumVars, fX, fParam);
        fX[i] = fXX - H;
        Fm = MyFxDelegate(nNumVars, fX, fParam);
        fX[i] = fXX;
        fDeriv[i] = (Fp - Fm) / 2 / H;
        fDerivNorm += pow(fDeriv[i], 2);
    }
    fDerivNorm = sqrt(fDerivNorm);
}
 
bool LinSearch_DirectSearch(int nNumVars, double fX[], double fParam[], double fLambda, double fDeltaX[], double InitStep, double MinStep)
{
    double F1, F2;
 
    F1 = MyFxEx(nNumVars, fX, fParam, fDeltaX, fLambda);
 
    do
    {
        F2 = MyFxEx(nNumVars, fX,  fParam, fDeltaX, fLambda + InitStep);
        if (F2 < F1)
        {
            F1 = F2;
            fLambda += InitStep;
        }
        else
        {
            F2 = MyFxEx(nNumVars, fX, fParam, fDeltaX, fLambda - InitStep);
            if (F2 < F1)
            {
                F1 = F2;
                fLambda -= InitStep;
            }
            else
            {
                InitStep /= 10;
            }
        }
    } while (!(InitStep < MinStep));
 
    return true;
 
}
 
double CalcOptim(int nNumVars, double fX[], double fParam[], double fEpsFx, int nMaxIter, int nIter,string sErrorMsg,double MyFxDelegate )
{
 
    int i;
    double fDeriv[2];
    double fDerivOld [2];
    double F, fDFNormOld, fLambda, fLastF, fDFNorm = 0;
 
    //MyFxDelegate = MyFx;
 
    // расчёт функции в определённой точке
    fLastF = MyFx(nNumVars, fX, fParam);
 
    GetGradients(nNumVars, fX,  fParam, fDeriv,fDFNorm);
 
    fLambda = 0.1;
    if (LinSearch_DirectSearch(nNumVars,  fX, fParam, fLambda,  fDeriv, 0.1, 0.000001))
    {
        for (i = 0; i < nNumVars; i++)
        {
            fX[i] += fLambda * fDeriv[i];
        }
    }
    else
    {
        sErrorMsg = "Ошибка! Не удалось выполнить операцию!";
        return fLastF;
    }
 
    nIter = 1;
    do
    {
        nIter++;
        if (nIter > nMaxIter)
        {
            sErrorMsg = "Достигнуто максимальное значение итераций";
            break;
        }
        fDFNormOld = fDFNorm;
        for (i = 0; i < nNumVars; i++)
        {
            fDerivOld[i] = fDeriv[i]; // Сохранить старый градиент
        }
        GetGradients(nNumVars, fX, fParam, fDeriv, fDFNorm);
        for (i = 0; i < nNumVars; i++)
        {
            fDeriv[i] = pow((fDFNorm / fDFNormOld), 2)* fDerivOld[i] - fDeriv[i];
        }
        if (fDFNorm <= fEpsFx)
        {
            sErrorMsg = "Градиент подходит под условия";
            break;
        }
 
        fLambda = 0;
 
        if (LinSearch_DirectSearch(nNumVars, fX, fParam, fLambda, fDeriv, 0.1, 0.000001))
        {
            for (i = 0; i < nNumVars; i++)
            {
                fX[i] += fLambda * fDeriv[i];
            }
            F = MyFx(nNumVars, fX, fParam);
            if (abs(F - fLastF) < fEpsFx)
            {
                sErrorMsg = "Последовательные значения функции отвечают критериям конвергенции";
                break;
            }
            else
            {
                fLastF = F;
            }
 
        }
        else
        {
            sErrorMsg = "Ошибка! Не удалось выполнить операцию!";
            break;
        }
    } while (true);
 
    return fLastF;
 
}
 
 
 
int _tmain(int argc, _TCHAR* argv[])
{
    
    double fX[] = {0, 0};
    double fParam[] =  { 0, 0 };
    int nIter = 0;
    int nNumVars = 2;
    int nMaxIter = 100;
    double fEpsFx = 0.0000001;
    int i;
    double fBestF;
    string sAnswer;
    string sErrorMsg = "";
    cout<<"Метод сопряжённых градиентов(Метод Флетчера-Ривса)"<<endl;
    fBestF = CalcOptim(nNumVars, fX, fParam, fEpsFx, nMaxIter, nIter, sErrorMsg, MyFx);
    cout<<"Значение функции = "<<fBestF;
    cout<<"Количество итераций = "<<nIter;
    
    return 0;
}