Note: This page has been translated by MathWorks. Please click here

To view all translated materials including this page, select Japan from the country navigator on the bottom of this page.

To view all translated materials including this page, select Japan from the country navigator on the bottom of this page.

Solve stiff differential equations — low order method

```
[t,y] =
ode23s(odefun,tspan,y0)
```

```
[t,y] =
ode23s(odefun,tspan,y0,options)
```

```
[t,y,te,ye,ie]
= ode23s(odefun,tspan,y0,options)
```

`sol = ode23s(___)`

`[`

,
where `t`

,`y`

] =
ode23s(`odefun`

,`tspan`

,`y0`

)`tspan = [t0 tf]`

, integrates the system of
differential equations $$y\text{'}=f\left(t,y\right)$$ from `t0`

to `tf`

with
initial conditions `y0`

. Each row in the solution
array `y`

corresponds to a value returned in column
vector `t`

.

All MATLAB^{®} ODE solvers can solve systems of equations of
the form $$y\text{'}=f\left(t,y\right)$$,
or problems that involve a mass matrix, $$M\left(t,y\right)y\text{'}=f\left(t,y\right)$$.
The solvers all use similar syntaxes. The `ode23s`

solver
only can solve problems with a mass matrix if the mass matrix is constant. `ode15s`

and `ode23t`

can
solve problems with a mass matrix that is singular, known as differential-algebraic
equations (DAEs). Specify the mass matrix using the `Mass`

option
of `odeset`

.

`[`

additionally
finds where functions of (`t`

,`y`

,`te`

,`ye`

,`ie`

]
= ode23s(`odefun`

,`tspan`

,`y0`

,`options`

)*t*,*y*),
called event functions, are zero. In the output, `te`

is
the time of the event, `ye`

is the solution at the
time of the event, and `ie`

is the index of the triggered
event.

For each event function, specify whether the integration is
to terminate at a zero and whether the direction of the zero crossing
matters. Do this by setting the `'Events'`

property
to a function, such as `myEventFcn`

or `@myEventFcn`

,
and creating a corresponding function: [`value`

,`isterminal`

,`direction`

]
= `myEventFcn`

(`t`

,`y`

).
For more information, see ODE Event Location.

returns
a structure that you can use with `sol`

= ode23s(___)`deval`

to evaluate
the solution at any point on the interval `[t0 tf]`

.
You can use any of the input argument combinations in previous syntaxes.

`ode23s`

is based on a modified Rosenbrock
formula of order 2. Because it is a single-step solver, it may be
more efficient than `ode15s`

at solving problems
that permit crude tolerances or problems with solutions that change
rapidly. It can solve some kinds of stiff problems for which `ode15s`

is
not effective. The `ode23s`

solver evaluates the
Jacobian during each step of the integration, so supplying it with
the Jacobian matrix is critical to its reliability and efficiency [1].

[1] Shampine, L. F. and M. W. Reichelt, “The
MATLAB ODE Suite,” *SIAM Journal on Scientific
Computing*, Vol. 18, 1997, pp. 1–22.

Was this topic helpful?