Nlsolve for find root for a system of ODE for an implicit runge kutta method

No, I haven’t. It’s too much work for me to figure out what you mean without a mathematical specification of the problem and code that doesn’t run. Maybe I can give you a complete working version of something similar for reference.