function filt=R(m,wc)


   
   if  m==0 
      
   filt=wc*wc/2./pi; 
   
   
 
   elseif rem(m,2)==0
    
    filt=0;
    
  
    else 
    
      filt=-2*wc*wc/m/m/pi/pi/pi;
   
   
end