The following ROOT page gives examples to implement numerical integration in ROOT:
I would like to implement my own integrator function and compare to the results from the various gsl,ROOT integrator functions(ADAPTIVE,SINGULAR,etc…)
The function I have is of the form:
double dgauss(double (*f)(double&), double &a, double &b, double &eps).
Is there a way to wrap this similar to the way Functor works?
I would want something like: