close all;
%3. gyakorlat: Van der Pol oszcillátor
%% Idõ inicializálás:
N = 10000; %eddig fut
h = 0.01; %lépésköz
t =0:h:N*h; %idõ: szebb lesz itt a görbe ,ha kézzel állítom be a lépésközt

for mu = 0:0.1:5 % mu=0:0.01:1 %különbözõ mu paraméterekre számolom ki az egyenletet
    
    VanDerPol_egyenlet=(@(t, xy)[xy(2); mu*(1-xy(1)^2)*xy(2)-xy(1)]); 
%több kezdeti értékkel kiszámolom az egyenletet
[t, fi] = ode45(VanDerPol_egyenlet, t, [0.5 0.5]);
plot(fi(:,1), fi(:,2), 'b') %különbözõ színnel ábrázolom
hold on %egymásra rajzolom a görbéket
[t2, fi2] = ode45(VanDerPol_egyenlet, t, [0.1 0.1]);
plot(fi2(:,1), fi2(:,2), 'r')
[t3, fi3] = ode45(VanDerPol_egyenlet, t, [2.5 2.5]);
plot(fi3(:,1), fi3(:,2), 'g')
xlim([-4, 4]);
ylim([-8, 8]);
title(['mu=',num2str(mu)]); %aktuális paramétert kiírom az ábra fölé
hold off

   pause(0.01); %késlsltetés a rajzolásnál, animáció lesz
  % waitforbuttonpress; %kattintásra vált
end