%phase_focus_lin.m

clear all
clc 

%Stable
mu=1;
%Unstable
%mu=-1;

%Equilibrium 
x1=0;
x2=0;

A=[-mu*(x1^2+x2^2)-2*mu*x1^2 -1-mu*2*x1*x2; 1-2*mu*x1*x2 -mu*(x1^2+x2^2)-2*mu*x2^2]

[V,D]=eig(A)