is there any other command besides intsimp that allows to perform a standard integral but allowing for more inputs in the integrad procedure besides the integrating variable?
The inthp4 function will integrate a function over a specified range. It is a little more difficult to use, because you need to place your extra arguments inside of DS structures and pass those into your function.
Since you have the source code for intsimp (or alternatively intquad1), you could make a new version that passes two inputs to the function being integrated. You would need to change the procedure prototype from this:
to something like:
and then inside the new function intsimp2args, you would add this new argument, extra_arg to the inputs being passed to your function, which is called f inside of intsimp. For example, this line:
sfa = f(xl);
would need to be changed to:
sfa = f(xl, extra_input);