Numerical integration by using regular Gaussian quadrature scheme (Fun defined only with double(s)).

function QuadGauss(Fun: TRealFunction; lb: double; ub: double; const Constants: TVec; Const ObjConst: Array of TObject; const BasePoints: TVec; const Weights: TVec; Parts: Integer): double; overload;