I have a large model with many variables (about a hundred) and many differential equations (non-linear). What is the best way to handle a system this size in Mathematica?
Note that I am not asking about performance (100 variables is not that much). My concern is with handling the code in a manageable way. Each variable and each equation has a meaning and some particular structure.
If I were doing this in another language, I would create structs or types to hold different groups of properties together. But in Mathematica I have no idea how to handle the complexity.
Is there a recommended way to handle large systems of ODEs in Mathematica in an structured fashion?