% DCS 2021
% z Transform

z = tf( 'z' )

% 2nd-order

A = [ -2 4  ; 1 0 ];

M1 = z*eye(2) - A 

M2 = inv( M1 )

M3 = M2 *[ 1 ; 0  ]

M4 = [ 1 3 ] * M3


% 3th-order

A = [ -2 4 -6 ; 1 0 0 0; 0 1 0 0; 0 0 1 0 ];

M1 = z*eye(4) - A 

M2 = inv( M1 )

M3 = M2 *[ 1 ; 0 ; 0 ; 0 ]

M4 = [ 1 3 5 7] * M3


% 4th-order

A = [ -2 -4 6 8 ; 1 0 0 0; 0 1 0 0; 0 0 1 0 ];

M1 = z*eye(4) - A 

M2 = inv( M1 )

M3 = M2 *[ 1 ; 0 ; 0 ; 0 ]

M4 = [ 1 3 5 7] * M3
