function [om,m0,m1]=transfer(h)
%
% given the filter h 
% transfer functions m0 and m1 are computed
%
     om=linspace(0,pi);
     k=length(h);
     in=1:k;
     aq=(-1).^(in-1);
     g=h(k:-1:1).*aq;
     m0=zeros(1,100);
     m1=m0;
     for j=1:k
        m0=m0+h(j)*exp(-i*(j-1)*om);
        m1=m1+g(j)*exp(-i*(j-k+1)*om);
     end
