Which matlab ode solver to use

Select the China site in Chinese or English for best site performance. Other MathWorks country sites are not optimized for visits from your location. Toggle Main Navigation. Search Answers Clear Filters. Answers Support MathWorks. Search Support Clear Filters. Support Answers MathWorks. Search MathWorks. MathWorks Answers Support. Open Mobile Search. Trial software. You are now following this question You will see updates in your activity feed.

You may receive emails, depending on your notification preferences. What are the general guidelines for choosing ODE solvers for simulating a model using Simulink? Show older comments. Documentation Help Center Documentation. This page contains two examples of solving nonstiff ordinary differential equations using ode For most nonstiff problems, ode45 performs best.

However, ode23 is recommended for problems that permit a slightly cruder error tolerance or in the presence of moderate stiffness. Likewise, ode can be more efficient than ode45 for problems with more stringent error tolerances or when the ODE function is computationally expensive to evaluate. If the nonstiff solvers take a long time to solve the problem or consistently fail the integration, then the problem might be stiff.

Rewrite this equation as a system of first-order ODEs by making the substitution. The resulting system of first-order ODEs is. The general functional signature of an ODE function is. That is, the function must accept both t and y as inputs, even if it does not use t for any computations. As an aside, here is an interesting fact about higher order Runge-Kutta methods. Classical Runge-Kutta required four function evaluations per step to get order four.

Dormand-Prince requires six function evaluations per step to get order five. You can't get order five with just five function evaluations. And then, if we were to try and achieve higher order, it would take even more function evaluations per step. Let's use ODE45 to compute e to the t. We can ask for output by supplying an argument called tspan. Zero and steps of 0. If we supply that as the input argument to solve this differential equation and get the output at those points, we get that back as the output.

And now here's the approximations to the solution to that differential equation at those points. If we plot it, here's the solution at those points. And to see how accurate it is, we see that we're actually getting this result to nine digits. ODE45 is very accurate. Let's look at step size choice on our problem with near singularity, is a quarter. The differential equation is y prime is 2 a-t y squared.

We let ODE45 choose its own step size by indicating we just want to integrate from 0 to 1. We capture the output in t and y and plot it. Advanced event location — restricted three body problem. Nonstiff problem — Euler equations of a rigid body without external forces. Freeman, San Francisco, Malcolm, and C. Moler, and S. Choose a web site to get translated content where available and see local events and offers.

Based on your location, we recommend that you select:. Select the China site in Chinese or English for best site performance. Other MathWorks country sites are not optimized for visits from your location.

Toggle Main Navigation. Search MathWorks. Open Mobile Search. Off-Canvas Navigation Menu Toggle. Main Content. Choose an ODE Solver Ordinary Differential Equations An ordinary differential equation ODE contains one or more derivatives of a dependent variable, y , with respect to a single independent variable, t , usually referred to as time. If there is a mass matrix, it must be constant. References [1] Shampine, L. Select a Web Site Choose a web site to get translated content where available and see local events and offers.