%Guia 1. Problema 1.4. G. Uicich
%Ecuacion de Van der Pol
clc;
clear;
figure;
hold on;
grid on;
axis([-3 3 -6 6]);
%Intervalo de integracion Tmax
Tmax=100;
%Vector de condiciones iniciales
for i=1:10
%r=3;
x1=[0.05 -1 1 2 -2.5 -0.05 2.5 -2.5 -0.5 0.5];
x2=[0.05 -3 1 -2 2.5 -0.05 4 -4 0.5 -0.5];
color=['b' 'r' 'k' 'g' 'y' 'm' 'y' 'k' 'g' 'b'];
options=odeset('RelTol',1e-3);
[T,X]=ode15s('ecdif2',[0 Tmax],[x1(i) x2(i)],options);
plot(X(:,1),X(:,2),color(i));
end
title('Espacio de Estados');xlabel('X1');ylabel('X2');
hold off;
end