global u
global mi step
step = 0.01;
u = input('Zvolte budiaci signal u: ');
mi = input('Zvolte parameter mi: ');
t0 = input('Zadajte pociatocny cas t0: ');
t1 = input('Zadajte konecny cas t1: ');
rTemp = roots([1 -2*mi 1]);
s1 = rTemp(1,1);
s2 = rTemp(2,1);

result1 = @(t,x)[x(2); mi*(1 - x(1)*x(1))*x(2) - x(1)+u];
result2 = @(t,x)[x(2); -x(1)+mi*x(2)+u];

f = figure;
tab = uitabgroup(f);

tab1 = uitab(tab, 'Title', 'Priebehy-nelinearne');
hax1 = axes('Parent', tab1);

[t,x]=ode45(result1, [t0 t1], [1 0]);
[t_rk,x_rk]=rk4(result1, [0 20], [1 0], 0.01);
subplot(2,1,1);
plot(t, x);
title('VanDerPol Nelineárny - ode45');
legend('x1(t)','x2(t)');
xlabel('t(s)');
ylabel('x1(t) x2(t)');
subplot(2,1,2);
plot(t_rk, x_rk);
title('VanDerPol Nelineárny - Runge Kutta');
legend('x1(t)','x2(t)');
xlabel('t(s)');
ylabel('x1(t) x2(t)');

tab2 = uitab(tab, 'Title', 'Priebehy-linearne');
hax2 = axes('Parent', tab2);
[t,x]=ode45(result2, [t0 t1], [1 0]);
[t_rk,x_rk]=rk4(result2, [0 20], [1 0], 0.01);
subplot(2,1,1);
plot(t, x);
title('VanDerPol Lineárny - ode45');
legend('x1(t)','x2(t)');
xlabel('t(s)');
ylabel('x1(t) x2(t)');
subplot(2,1,2);
plot(t_rk, x_rk);
title('VanDerPol Lineárny - Runge Kutta');
legend('x1(t)','x2(t)');
xlabel('t(s)');
ylabel('x1(t) x2(t)');

tab3 = uitab(tab, 'Title', 'FP-Van Der Pol - Linearny');
hax3 = axes('Parent', tab3);
fp(result2, [0 50], [-3 3], [-3 3], 10, 10);
title('Van Der Pol - Linearny System');
axis([-4 4 -4 4]);

tab4 = uitab(tab, 'Title', 'FP-Van Der Pol - Nelinearny');
hax4 = axes('Parent', tab4);
fp(result1, [0 50], [-3 3], [-3 3], 10, 10);
title('Van Der Pol - Nelinearny System');

%Fázovy portrét
% reálne císla
if (isreal(s1) == 1) && (isreal(s2) == 1)
    % nestabilný uzol
    if (s1 > 0) && (s2 > 0)
        title('Van Der Pol - Lineárny Systém - Nestabilný uzol');
    end

    % stabilný uzol
    if (s1 < 0) && (s2 < 0)
        title('Van Der Pol - Lineárny Systém - Stabilný uzol');
    end

    % sedlo
    if (s1 < 0) && (s2 > 0)
        title('Van Der Pol - Lineárny Systém - Sedlo');
    end

    % sedlo
    if (s1 > 0) && (s2 < 0)
        title('Van Der Pol - Lineárny Systém - Sedlo');
    end
end

% komplexné cisla
if (isreal(s1) == 0) && (isreal(s2) == 0)
    % stabilné ohnisko
    if (real(s1) < 0) && (real(s2) < 0)
        title('Van Der Pol - Lineárny Systém - Stabilné ohnisko');
    end

    % nestabilné ohnisko
    if (real(s1) > 0) && (real(s2) > 0)
        title('Van Der Pol - Lineárny Systém - Nestabilné ohnisko');
    end

    % stred
    if (real(s1) == 0) && (real(s2) == 0)
        title('Van Der Pol - Lineárny Systém - Stred');
    end
end