Frequency and Phase response of an IIR system whose system function is given H(z)

My system function is of the form below,

I want to plot its magnitude and phase response. w should be between -pi and pi. How can I plot?

Answer 1

[Paul Hoffrichter]

Here is an example to plot the magnitude and phase frequency response

>> b = [0.0563   -0.0009   -0.0009    0.0563];
>> a = [1.0000   -2.1291    1.7834   -0.5435];
>> freqz(b,a)

https://www.mathworks.com/help/signal/ref/freqz.html