czas = 0:0.1:15;
k_t = [0.7, 1.6];
T_t = [1.3, 0.9];
kolor = ['-g', '-r'];
tlumienie1 = 3;
tlumienie2 = 0.3;

for i=1:2
    k = k_t(i);
    T = T_t(i);
    
    licznik = [0, k];
    mianownik = [T, 1];
    
    [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
    [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
    
    figure(1)
    subplot(1, 2, 1)
    plot(czas_s, y_s, kolor(i))
    hold on
    
    subplot(1, 2, 2)
    plot(czas_i, y_i, kolor(i))
    hold on
end

for i=1:2
    k = k_t(i);
    T1 = T_t(i);
    T2 = T1 + 1;
    
    licznik = [0, 0, k];
    mianownik = [T1*T2, T1+T2, 1];
    
    [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
    [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
    
    figure(2)
    subplot(1, 2, 1)
    plot(czas_s, y_s, kolor(i))
    hold on
    
    subplot(1, 2, 2)
    plot(czas_i, y_i, kolor(i))
    hold on
end

for i=1:2
    k = k_t(i);
    T = T_t(i);
    
    licznik = [0, 0, k];
    mianownik = [T^2, 2*tlumienie1*T, 1];
    
    [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
    [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
    
    figure(3)
    subplot(1, 2, 1)
    plot(czas_s, y_s, kolor(i))
    hold on
    
    subplot(1, 2, 2)
    plot(czas_i, y_i, kolor(i))
    hold on
end

for i=1:2
    k = k_t(i);
    T = T_t(i);
    
    licznik = [0, 0, k];
    mianownik = [T^2, 2*tlumienie2*T, 1];
    
    [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
    [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
    
    figure(4)
    subplot(1, 2, 1)
    plot(czas_s, y_s, kolor(i))
    hold on
    
    subplot(1, 2, 2)
    plot(czas_i, y_i, kolor(i))
    hold on
end