void bandPassInit() {
  w1=0 ;
  w2=0 ;
  float r=0.99 ;
  float f=57e3/fSample ;
  a1=-2*r*cos(f*2*PI) ;
  a2=r*r ;
  b0=1 ;
  b1=0 ;
  b2=-1 ;
  }
  
float bandPass(float x){
  float y ;
  w=x-a1*w1-a2*w2 ;
  y=b0*w+b1*w1+b2*w2 ;
  w2=w1 ;
  w1=w ;
  return y ;
  }
