By default, the ODE solver passes all components of the solution to the output function. odeprint displays the solution components on the screen. odephas2 and odephas3 produce two- and three-dimnesional phase plane plots, respectively. When you call the solver with no output arguments, it calls the default odeplot to plot the solution as it is computed. Four output functions are provided: odeplot, odephas2, odephas3, odeprint. If you specify an output function as the value of the OutputFcn property, the solver calls it with the computed solution after each time step. The values indicate which event the solver detected. Indices into the vector returned by the function specified in the Events option. Solutions that correspond to events in sol.xe.

sol.xe(end) contains the exact point of a terminal event, if any. Points at which events, if any, occurred. If you specify the Events option and events are detected, sol also includes these fields: The structure sol always includes these fields:Įach column sol.y(:,i) contains the solution at sol.x(i). You must pass odefun as a function handle. Returns a structure that you can use with deval to evaluate the solution at any point on the interval.

direction(i) = 0 if all zeros are to be computed (the default), +1 if only the zeros where the event function increases, and -1 if only the zeros where the event function decreases.Ĭorresponding entries in TE, YE, and IE return, respectively, the time at which an event occurs, the solution at the time of the event, and the index i of the event function that vanishes.

isterminal(i) = 1 if the integration is to terminate at a zero of this event function and 0 otherwise.

Do this by setting the 'Events' property to a function, e.g., events or and creating a function = events( t, y). For each event function, you specify whether the integration is to terminate at a zero and whether the direction of the zero crossing matters. Solves as above while also finding where functions of, called event functions, are zero. Use options = as a place holder if no options are set. to the function odefun, whenever it is called. Solves as above, passing the additional parameters p1,p2. Commonly used properties include a scalar relative error tolerance RelTol ( 1e-3 by default) and a vector of absolute error tolerances AbsTol (all componen ts are 1e-6 by default). Solves as above with default integration parameters replaced by property values specified in options, an argument created with the odeset function. To obtain solutions at the specific times t0, t1.,tf (all increasing or all decreasing), use tspan = . Each row in the solution array Y corresponds to a time returned in column vector T. Function f = odefun(t,y), for a scalar t and a column vector y, must return a column vector f corresponding to. With tspan = integrates the system of differential equations from time t0 to tf with initial conditions y0. Optional parameters that the solver passes to odefun and all the functions specified in options. Optional integration argument created using the odeset function.

To obtain solutions at specific times (all increasing or all decreasing), use tspan =. ode15s and ode23t can solve problems with a mass matrix that is singular, i.e., differential-algebraic equations (DAEs).Ī vector specifying the interval of integration. The ode23s solver can solve only equations with constant mass matrices.

All solvers solve systems of equations in the form or problems that involve a mass matrix. Where solver is one of ode45, ode23, ode113, ode15s, ode23s, ode23t, or ode23tb.Ī function that evaluates the right-hand side of the differential equations. = solver(odefun,tspan,y0,options,p1,p2.) Solve initial value problems for ordinary differential equations (ODEs) Ode45, ode23, ode113, ode15s, ode23s, ode23t, ode23tb (MATLAB Functions) MATLAB Function Reference

Matlab function handle ode