aa

function [Ad, Bd, Cd, Tt, Hdz, Yfz ] = lab4(A, B, C)

syms s;
Ts = -2/((s+1)*(s+2));
syms h;
eA = simplify(inv(s*sym(eye(size(A))) - A));
Ad = simplify(ilaplace(eA, s, h));
Bd = int((Ad*B), h, 0, h);
Cd = C;
syms z;
syms t;
syms k;
Tt = ilaplace((Ts/s), k*h);
Hdz = simplify((z-1)/z * ztrans(Tt, z));
Yfz = simplify(Hdz * (z / (z-1)));

end