Great article. I have used PROC MODEL for years to solve ODEs (especially SEIR models), and to fit the models to data. I will have to play with odeqn to see how it works with some difficult problems. I also know there is a ODE solver in PROC MCMC for Bayesian analysis, but I have not tried that feature yet.