In this video, the difference equation of a causal LTI discrete-time system is used to find the transfer function H(z) then the factored form of the transfer function using partial fractions expansion (residuez command).
Finally the zeros and poles were found(tf2zpk command), plotted (zplane command) and the stability assessed.