I am looking to solve large systems of non-linear ODEs. There appears to be a very large list of methods available varying in complexity, and I have a hard time searching through them and picking one. Are any of these methods preferred for large systems? Both speed and accuracy play large roles, so I'm hoping that there are methods that are in general considered to be better for large systems. Some additional information: The systems usually consist of ~100 ODEs that are quite heavily linked, usually consisting of a lot of quartic terms. (2-loop renormalization group equations)Thanks
Numerically solving systems of about 100 ODEs
ode
100 equations is not a particular large system. There are certainly many good integrators for this out there -- starting with Matlab's ode45 which should have no problems with a system of 100 equations.The challenge with ODEs is not typically the size, but the character. For example, is your system stiff? If so, you may want to look at CVODE. Do you need to preserve certain invariants? If so, you may want to look into symplectic integrators.