## Problem #14

For the Van der Pol equation,

(a) Find all the critical points (equilibrium solutions).

(b) Use a computer to draw a direction field and phase portrait for the
system.

(c) From the plot(s) in part (b), determine whether each critical point is
asymptotically stable, stable, or unstable, and classify it as to type.

## Part (a)

By solving *x'* = 0 and *y'* = 0, we find that the only
critical point of this system is (0, 0).

## Part (b)

The module dirfieldAuto.m will plot a
direction field for any autonomous system (meaning, without explicit use of
'*t*') of two equations in *x* and
*y* and gives the user the option to plot solution curves to
initial value problems. The difference between this module and 'dirfield2d'
is that 'dirfield2d' solves linear systems exactly, where 'dirfieldAuto' uses
one of MATLAB's built-in ODE solvers to calculate a solution. It is slower and less exact than dirfield2d, but for non-linear systems it's the best we
can do.

` >> dirfieldAuto
*-----------------------------------------------------------*
| This module will plot a direction field for an autonomous |
| system of equations x' = f(x,y), y' = g(x,y). |
*-----------------------------------------------------------*
Enter x' as a function of x and y: x' => y
Enter y' as a function of x and y: y' => (1-x^2)*y-x
Enter the smallest x value => -3
Enter the largest x value => 3
Enter the smallest y value => -3
Enter the largest y value => 3
Warning: Divide by zero.
> In `__dirfieldAuto at 86__
Warning: Divide by zero.
> In __dirfieldAuto at 87__
>> Would you like to plot a solution curve on this direction field? (y or n) => n
>>

Looking at the direction field, it seems that the origin is an unstable
critical point, as it appears that the trajectories all leave the origin
in a spiraling motion. However, we'll need to look at some actual
trajectories to be sure that this is the case.

Now, we can run 'dirfieldAuto' again, answering 'y' to the final question
to plot some trajectories. the figure below was created using 8 sets of initial conditions.

**Warning:** This module takes awhile to calculate solutions.
Be patient.

`
`

## Part (c)

From the plot in part (b), we see that the origin is an unstable
critical point. In fact, all trajectories approach the pillow-shaped
closed curve shown on the phase portrait, whether they start inside or outside
of that curve.