procedure Eckerle4(const B: TVec; const X: TVecList; const y: TVec);
begin
  //B[0]/B[1] * Exp(-0.5*Sqr((X-B[2])/B[1]));

  y.Normalize(X[0], B[2], B[1]);
  y.Sqr;
  y.Scale(-0.5);
  y.Exp;
  y.Scale(B[0]/B[1]);
end;

void __fastcall Eckerle4(TVec* const b, TVecList* const x, TVec * const y)
{
  double* B = b->PValues1D(0);

  y.Normalize(X[0], B[2], B[1]);
  y.Sqr();
  y.Scale(-0.5);
  y.Exp();
  y.Scale(B[0]/B[1]);
}