Dear Michael, thanks for the code. I will try it on Mathematica directly, however I have a final question,
you say a is not there in the ODE, which is true, however, it is in your code:
Subscript[a, b
iand it is further down given along with b and c, which are in the equation, as a defined value (1). Where is a in the system of the ODEs, if it is included here in your code, and if it is not-related to the ODE, what is the a for?
Thanks!
PS: I take I can replace b and c with some different values.
Cheers