syms n z;
% note that sym(1.2) is used for 1.2delta(n)
HZ = sym(1.2) + ztrans(0.5*(-0.5)^n*heaviside(n) -0.6*(-0.6)^n*heaviside(n) );
H_inverse = 1/HZ;



fprintf('The inverse of H(z) is \n')
pretty(simplify(H_inverse))

fprintf('h inverse as function of n\n')
iztrans(H_inverse)

