%% L04_Invztransform
%% ==============
%% Script to find inverse z-transform using symbolic toolbox (maple).
%% Refer to: Lecture04, Example 1, slide 16-18
%%
%% Created  : 27 July 2005
%% Modified : 27 July 2005
%%
%% Copyright (c) 2005 Salman Durrani .
%% ------------------------------------------------------------

%% clear matlab memory 
clear all;
clc

syms w a c n z

%% take inverse z transform
H = (z+ 0.4)/((z-0.2)*(z-0.5));
Ans8 = maple('invztrans',H,z,n)
pretty(Ans8)

%% plot h[n]
n=[0:1:10];
h = -10.*((1/5).^n) + 6.*((1/2).^n);
h(1) = h(1)+4;
stem(n,h,'filled')
xlabel('Sample n')
ylabel('Amplitude')