t0=0;
figure(1)
hold off

gamma=0.1;
sella=[-1 0]
V1=[-gamma/2-sqrt(gamma^2+8)/2 1]
V2=[-gamma/2+sqrt(gamma^2+8)/2 1]

tf=20
x0=sella + V2*1e-3
[t,x]=ode45('esesep',[t0 tf],x0);
plot(x(:,1),x(:,2),'r')
hold on

tf=5
x0=sella + V1*1e-3
[t,x]=ode45('esesepm',[t0 tf],x0);
plot(x(:,1),x(:,2),'g')

tf=5.5
x0=sella - V2*1e-3
[t,x]=ode45('esesep',[t0 tf],x0);
plot(x(:,1),x(:,2),'r')

tf=9.5
x0=sella - V1*1e-3
[t,x]=ode45('esesepm',[t0 tf],x0);
plot(x(:,1),x(:,2),'g')

axis([-3 3 -4 4])
warning=' stabile, poi inst, stab, inst'
gtext('stabile','color','green')
gtext('instabile','color','red')
gtext('stabile','color','green')
gtext('instabile','color','red')
warning=' 2 A poi 2 B'
gtext('A')
gtext('A')
gtext('B')
gtext('B')
pause
xlabel('x')
ylabel('y=dx/dt')
title(' gamma=0.1')
print -dpsc colesesepar.ps
print -dps bwesesepar.ps
