function filt=Rhamm(w,wc)


   
   if w <= wc 
      
      filt=w   *.54   +  .46*w*cos(pi*w/wc);
      
   else 
      
       filt=0.;

   
end
