%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%         2.gyakorlat         %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Nemlin. din. rsz-ek alapjai %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%          2016.09.21.        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Paraméterek:

kezd1 =[0.5 0.5];
kezd2 =[0.1 0.1];
kezd3 =[2.5 2.5];
tmax = 100;
mu = 5;
for mu = 0:10:100 %0:0.1:5    % 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)];                    % Adjuk meg az egyenletet!
    % több kezdeti értékkel kiszámolom az egyenletet
    
    %[t, fi] = ode45(VanDerPol_egyenlet, [0,tmax], kezd1);
    [t, fi] = ode15s(VanDerPol_egyenlet, [0,tmax], kezd1);
    %[t2, fi2] = ode45(VanDerPol_egyenlet, [0,tmax], kezd2);
    [t2, fi2] = ode15s(VanDerPol_egyenlet, [0,tmax], kezd2);
    %[t3, fi3] = ode45(VanDerPol_egyenlet, [0,tmax], kezd3);
    [t3, fi3] = ode15s(VanDerPol_egyenlet, [0,tmax], kezd3);


subplot(2,2,1)
plot(fi(:,1), fi(:,2), 'b') %különbözõ színnel ábrázolom (x,y)
hold on
plot(fi2(:,1), fi2(:,2), 'r')
plot(fi3(:,1), fi3(:,2), 'g')
xlim([-4, 4]);
%ylim([-8, 8]);
xlabel('x');
ylabel('y=dx');
title(['mu=',num2str(mu)]); %aktuális paramétert kiírom az ábra fölé
hold off

subplot(2,2,2)
plot(t, fi(:,2), 'b'); %(t,y) ábrázolás
%plot(t(1:i), fi(1:i,2), 'b') %mindig a adott idõpillanatig rajzolok ki
hold on %a hold on-t minden subplotra külön ki kell adni
plot(t2, fi2(:,2), 'r');
plot(t3, fi3(:,2), 'g');
xlim([0 100]);
%ylim([-8, 8]);
xlabel('t');
ylabel('y=dx');
hold off

subplot(2,2,3)
plot(fi(:,1), t, 'b'); %(x,t) ábrázolás
%plot(fi(1:i,1), t(1:i), 'b') %mindig a adott idõpillanatig rajzolok ki
hold on
plot(fi2(:,1), t2, 'r');
plot(fi3(:,1), t3, 'g');
xlim([-4, 4]);
ylim([0 100]);
xlabel('x');
ylabel('t');
hold off
        
waitforbuttonpress();
pause(0.01);

end