Marquardt

procedure Marquardt(

  Func: TFuncNVar;

  HessGrad: THessGrad;

  X: TVector;

  Lb: Integer;

  Ub: Integer;

  MaxIter: Integer;

  Tol: Float;

  var F_min: Float;

  G: TVector;

  H_inv: TMatrix;

  var Det: Float

);