s = tf('s');


% ODPOWIEDZI NA SKOK - ROZNE ROZMIESZCZENIE BIEGUNOW
B_Array = [15, -2, 1, -1];
C_Array = [2, -8, 15, 6];

for i = 3

    B = B_Array(i);
    C = C_Array(i);
    sys = (1/(s^2+B*s+C));
    sys_txt = sprintf('$\\frac{1}{s^2 %+.ds %+.d}$', B, C);

    figure(i);
    [y, t] = step(sys);
    f = fit(t,y,'exp1')
    grid on;
    title(['K(s) = ', sys_txt], 'FontSize', 25, 'interpreter', 'latex');
    hold on;

    expp = a*exp(-b*t)+c;
    plot(t, f,'Color',[0,0,0]);