how to construct a transfer function H(z)

How to use matlab to construct a transfer function H(Z) with four poles at 0.99<(+88,-88),0.98<(+92,-92),the numerator is unity, so the system function is all pole. also, how to plot the magnitude square of the transfer function 0.99,0.98 is the radius of unit circle in z domain 88,92 is the angle of pole thanks