BEGIN{
  fs=16384*1.64676
}
!/x/{
  e=$2-cos($1/65536.0*3.1415926)*fs
  if (e>emax) emax=e
  if (e<emin) emin=e
  e2sum = e2sum + e*e
  n = n+1
}
END{
  printf "test covers %d points, full scale is %d bits\n", n, fs
  if (-emin > emax) emax=-emin
  rms=sqrt(e2sum/n)
  printf "peak error %6.2f bits, %5.3f %%\n", emax, emax*100/fs
  printf "rms  error %6.2f bits, %5.3f %%\n", rms,  rms*100/fs
}