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)