Passing data to the ODE integrator in user-defined functions for double integration

before you get a double integration working, please get a single integration working which is called from inside a user function. I don’t have an immediate example which I can refer you to, but this more simpler exercise should let things sink in, I hope.