Аппроксимация с alglib

#include <vcl.h>
#pragma hdrstop
#include <iostream.h>
#include <fstream.h>
#include <math.h>
#include <iomanip>
#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
//---------------------------------------------------------------------------
#include "interpolation.h"  
#include "stdafx.h"
#include "optimization.h"
 
using namespace alglib;
 
double x[31]={0};
double y[31]={0};
real_1d_array c="[1,0]";
 
double Get_Array() //здесь я читаю из файла свои данные, всё ок, работает
{
 
    char buf[1000]={0};
    FILE *tab;
 
 
    float num1;
    float num2;
 
    strcpy(buf,"D:\\pub\\Bagrova\\Proc_Fit\\sGauss.dat");
    tab=fopen(buf,"r");
 
    if (tab)
    {
        fgets(buf,1000,tab);
        for (int i=0; i <= 30; i++)
        {
            fscanf(tab,"%f%f\n", &num1,&num2);
            x[i]=num1;
            y[i]=num2;
        }
    fclose(tab);
    }
    else printf("Error open");
    printf("\tX\ty\n");
    for (int i=0; i <= 30; i++)
    {
        printf("%f\t%f\n",x[i],y[i]);
    }
    getche();
    return 1;
}
void  gauss_func(const real_1d_array &c, const real_1d_array &x, double &func, void *ptr)
{
    func=exp(0-pow(x[0]-c[1],2)/(2*pow(c[0],2)))/pow(2*M_PI,0.5); //c[0] - sigma, c[1] - mu
}
 
void gauss_grad(const real_1d_array &x, double &func, real_1d_array &grad, void *ptr)
{
    func=exp(0-pow(x[0]-c[1],2)/(2*pow(c[0],2)))/pow(2*M_PI,0.5);
    grad[0]=pow(x-c[1],2)*exp(0-pow(x[0]-c[1],2)/(2*pow(c[0],2)))/(pow(2*M_PI,0.5)*pow(c[0],4))-exp(0-pow(x[0]-c[1],2)/(2*pow(c[0],2)))/(pow(2*M_PI,0.5)*pow(c[0],2));
    grad[1]=(x-c[1])*exp(0-pow(x[0]-c[1],2)/(2*pow(c[0],2)))/(pow(2*M_PI,0.5)*pow(c[0],3));
 
}
 
int lsFit_FG()
{
    real_2d_array X1 = "[[-3],[-2.8],[-2.6],[-2.4],[-2.2],[-2],"
    "[-1.8],[-1.6],[-1.4],[-1.2],[-1],[-0.8],[-0.6],[-0.4],[-0.2],[0],[0.2],[0.4],"
    "[0.6],[0.8],[1],[1.2],[1.4],[1.6],[1.8],[2.0],[2.2],[2.4],[2.6],[2.8],[3]]";
    real_1d_array Y1;
 
    Y1.setcontent(31,y);
 
    double epsf = 0;
    double epsx = 0.000001;
    ae_int_t maxits = 0;
    ae_int_t info;
    lsfitstate state;
    lsfitreport rep;
 
    lsfitcreatefg(x, y, c, true, state);
    lsfitsetcond(state, epsf, epsx, maxits);
    alglib::lsfitfit(state, gauss_func, gauss_grad);
    lsfitresults(state, info, c, rep);
    printf("%d\n", int(info)); // EXPECTED: 2
    printf("%s\n", c.tostring(1).c_str()); 
 
    return 0;
}

void gauss_grad(const real_1d_array &c,const real_1d_array &x, double &func, real_1d_array &grad, void *ptr)