z=tf('z')
delay=z^-15
licz=0.0242*z^2 + 0.0233*z
mian=(z-0.8356 - 0.1846i)*(z- 0.0000 -0.0000i)*(z-0.8356 +0.1846i)
Gz=(licz/mian)*delay
nmian=(z-p1)*(z-p2)*(z-p3)
nGz=(licz/mian)*delay