##

Plotting of integral curves

The output of `ode23`

for a set of two differential equations as in equation (5.2), say
and
where
is the time vector and
the corresponding states, can also be used for drawing the phase plane. To this end, you need to plot the solutions
and
in the [
,
]-plane. This can be achieved by plotting the second column of
versus the first column of
after you have used the command `ode23`

:

`>> plot(x(:,1),x(:,2))`

This gives a figure with one integral curve. If you want to draw several integral curves with different plot commands in the same figure, then you have to enter the command **hold on** after the first plot command (see Chapter 1). If you draw the integral curves and the direction field in one figure, then the direction vectors will be tangent to the integral curves. You can also draw an integral curve by assigning the property `odephas2`

to `OutputFcn' in the ode-solver. The command:

`>> ode23('deq',[0,1],[2,8],odeset('OutputFcn','odephas2'))`

results in a two-dimensional integral curve from the given initial condition.

Esteur
2010-03-22