clc;
l = 1;
m = conv([0.5 1], [1 2.5 2.5]);
t = 0:0.1:10;
yo= step(l,m,t);
plot(t,yo),grid;
[M i] = max(diff(yn));
yn(12);
[lp mp] = pade(0.29, 8);
la = 0.4 * lp;
ma = conv([0.58 1], [0.58 1]);
ya = step(la,ma,t);
plot(t,yo,t,ya),grid