clear all;  clf;
K=1;
s=tf('s');

T1=8;
figure()
T2=T1;

numerator = K;
denominator = [T1*T2 T1+T2 1];
sys = tf(numerator, denominator);
step(sys)
xlim([0 60]);

% Odczytane:
Tk=27; %27
Tk0=3.5;
Kk=1;

axis([-1 100 -0.05 1.05]);
hold on; grid on;
Sk=(Kk/(Tk*s+1))*exp(-s*Tk0);
step(Sk);
legend('Czlon dla T2=T1','Model K�pfm�llera')