For the differential equation

(a) Determine all critical points of the given system of equations.

(b) Find the corresponding linear system near each critical point.

(c) Find the eigenvalues of each linear system. What conclusions can you then
draw about the nonlinear system?

(d) Draw a phase portrait of the nonlinear system to confirm your conclusions,
or to extend them in those cases where the linear system does not provide
definite information about the nonlinear system.

The critical points of this system are (0,0) and (-1.19345, -1.4797), and thus they are both isolated critical points.

For functions *x' = 2x + y + xy ^{3} = F(x,y)* and

The corresponding linear system near the origin is

The eigenvalues and eigenvectors for this linear system are found using MATLAB's 'eig' command:

>> A=[2 1;1 -2]; >> [V,D]=eig(A) V = 0.2298 -0.9732 -0.9732 -0.2298 D = -2.2361 0 0 2.2361 >>

Thus, the eigenvectors and eigenvalues for the linear system about the origin are

The eigenvalues are real, distinct and have opposite signs, so the origin is an unstable saddle point.

The linear system that approximates the nonlinear system near the critical
point (*x*_{0}, *y*_{0}) = (-1.9345, -1.4797) is

We can be glad we have MATLAB to find the eigenvalues and eigenvectors of this system:

>> B=[-1.23982 -6.83922;2.4797 -0.80655]; >> [V2,D2]=eig(B) V2 = -0.8567 -0.8567 0.0271 + 0.5151i 0.0271 - 0.5151i D2 = -1.0232 + 4.1125i 0 0 -1.0232 - 4.1125i >>

Thus, we see that the eigenvalues and eigenvectors are

Because the eigenvectors are complex with negative real part, the critical point at (-1.19345, -1.4797) is a stable spiral point.

The module dirfieldAuto.m will plot a direction field for any 2x2 system of autonomous differential equations, and is particularly useful for nonlinear systems such as this one.

>> 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' => 2*x + y + x*y^3 Enter y' as a function of x and y" y' => x - 2*y - x*y Enter the smallest x value => -2.5 Enter the largest x value => 2.5 Enter the smallest y value => -2.5 Enter the largest y value => 2.5 Warning: Divide by zero. >In dirfieldAuto at 88Warning: Divide by zero. >In dirfieldAuto at 89>> Would you like to plot a solution curve on this direction field? (y or n) => y Enter initial condition: x(0) => 2 y(0) => 1 Would you like to plot another solution curve? (y or n) => y

We see that the origin is indeed a saddle point, and that there is an inward spiral around the point (-1.19345, -1.4797), as we predicted in the analysis above.