function filt=Rprzn(w,wc)


   
   if w <= wc/2.
      
      filt=w - 6*w*(w*w/(wc*wc))*(1.-w/wc);
      
   
elseif  ( w<=wc ) & (w >(wc/2)) 
   
   filt= 2.* w * (1.-w/wc) * (1.-w/wc) * (1.-w/wc);
   
else 
   
   filt=0.;
   
end
