y="[130.9045354801, 135.9066659041, 131.9036065061, 123.9053434201]";
x="[[32.13088888,128],[46.53452299,126],[49.49980573,128],[45.8625616,132]]";

c="[0.01,1000,0.1]"
cScale[0]=0.01;
cScale[1]=1000;
cScale[2]=0.1;
epsf=1e-10;
epsx=1e-10;
maxits=0;
diffstep=1e-7;
bndl="[-1e8,-1e8,-1e8]";
bndu="[-1e8,-1e8,-1e8]";

void function_cx_1_func(const real_1d_array &c, const real_1d_array & x, double &func, void *ptr) {
   func = c[0] * square(x[0] + c[1]) / square(1.0 + c[2] * x[1]);
}
void function_cx_1_grad(const real_1d_array &c, const real_1d_array & x, double &func, real_1d_array &grad, void *ptr) {
   grad[0] = square(x[0] + c[1]) / square(1 + c[2] * x[1]);
   grad[1] = 2 * c[0] * (x[0] + c[1]) / square(1 + c[2] * x[1]);
   grad[2] = -2 * c[0] * c[1] * square(x[0] + c[1]) / (square(1 + c[2] * x[1])*(1 + c[2] * x[1]));
}

lsfitcreatefg(x, y, c, diffstep, state);
lsfitsetcond(state, epsf, epsx, maxits);
lsfitsetbc(state, bndl, bndu);
lsfitsetscale(state, cScale);
alglib::lsfitfit(state, function_cx_1_func, function_cx_1_grad);
lsfitresults(state, info, c, rep);