>> n
n =
     4

>> m=8
m =
     8

>> Gp=tf((100000+n)/(100000-n), conv([1 0], conv([1 n], [1 n m])))
 
Transfer function:
             1
---------------------------
s^4 + 8 s^3 + 24 s^2 + 32 s
 
>> 
>> Gp_zpk=zpk(Gp)
 
Zero/pole/gain:
        1.0001
----------------------
s (s+4) (s^2 + 4s + 8)
 
>> sisotool(Gp)