## 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