syms n

delta = heaviside(n) - heaviside(n-1);
hn = 1.2.*delta + 0.5 * (-0.5).^n .* heaviside(n) - 0.6 * (0.2).^n .* heaviside(n);

ztransform = ztrans(hn);
ztransform_inverse = 1/ztransform;

hn_inverse = iztrans(ztransform_inverse)

