I am looking for some advice regarding using the built-in ode solver in EJS.

I am attempting to simulate cold matter interacting with light. My simulation is producing results which are of a good approximation for a small number of atoms (n=3 for example) however when I increase the value to n=15 for example the ode solver breaks down and does not produce any results.

I am just wondering if there is a limit to the complexity of the equations for which EJS can handle.

Ryan-- There is no inherent limit, other than memory and processor speed, on the complexity or number of ODE equations that EJS can solve. The rate functions must, of course, be without discontinuities. I have solved systems with a few hundred ODEs without difficulty. See the Fermi Pasta Ulm model and the falling chain model.