Use the following Matlab script to compute the coefficients for an FIR digital filter.
order $=40$
wc $=0.491$
type $=$ high';
b1 = fir1(order, wc, type);
order $=40$;
$w c=0.491 ;$
type $=$ "high";
$b 1=$ fir1(order, wc, type);
Use the Matlab function freqz to plot the magnitude and phase plots for the filter designed using this script which we will call $H(z)$ for this problem.