BFGS

procedure BFGS(

  Func: TFuncNVar;

  Gradient: TGradient;

  X: TVector;

  Lb: Integer;

  Ub: Integer;

  MaxIter: Integer;

  Tol: Float;

  var F_min: Float;

  G: TVector;

  H_inv: TMatrix

);