Howdy Folks,

I'm trying to help someone convert a ODE based model from one software
package to another. The other software package uses the LSODA solver
[1]. Her system is pretty stiff, and I could get it to solve for one
set of parameters using ode15s, but I had to switch to ode23s to get
it to solve with a different set of parameters. Based on her
experience with the LSODA solver, it seems more robust.

I was wondering if there was anyone out there using this solver in
Matlab with an implementation I could use for comparison purposes.








[1] https://computation.llnl.gov/casc/odepack/odepack_home.html