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.