%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%This program will take in wc and      %
%Output a low pass filter              %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [a,num,den] =first_low_pass(wc)

a=(1-sin(wc))/cos(wc);          %calculate a
num = (1-a)*[1,1];              %num vec. og trans
den = 2*[1,-a];                 %den vec of trans
freqz(num,den);                 %returns the frequency response 
                                %vector h and the corresponding 
                                %frequency vector w for the quantized filter
                                %num, plots frequency response.