The way you are doing it looks pretty good to me. Your residuals are down in the region of 10^-8 to10^-9, and given that your starting point is a numerical solution from odesolve anything better would be meaningless. The only thing I would do differently is that I would set up a vector of times, and then have one vectorized equation in the solve block rather than 12 separate equations.
↧
