Laborator TS

function [rez] = lab(A, B, C)
    s = sym('s');
    I = eye(2);
    D = s*I - A;
    D1 = inv(D);
    T = C*D1*B;
    T
    h = sym('h');
    Ad = ilaplace(D1, s, h);
    Ad
    t = sym('t');
    Ad1 = ilaplace(D1, s, t);
    S = Ad1*B;
    Bd = int(S, t, 0, h);
    Bd
    Cd = C
    z = sym('z');
    raport = (z-1) / z;
    Hs = T / s;
    k = sym('k');
    Fk = ilaplace(T, s, k*h);
    Hd = raport*ztrans(Fk, k, z);
    Hd
    yfz = Hd*(z / (z-1));
    yfz
    yfk = iztrans(yfz, z, k);
    yfk
end