JuliaLang/julia — pull request #615 by pao (2012-03-20 04:53:57+00:00)
Title: Stiff ODE solver + small fix

Finally figured out the error in the Rosenbrock solver after carefully auditing vs. Numerical Recipes. The Jacobian matrix estimator is not industrial-strength but seems to get the job done. Suggestions for a good way to handle providing an explicit function for the Jacobian matrix are welcome.

I also cleaned up a small fix I made earlier in ode23() since the promotion behavior in min/max has changed to make it unneeded.