block-computed two-pole FIR filter
This runs a block-computing version of the reson filter.
From the original documentation:
float breson(x,rsnvals,array,num) float *x,*array,*rsnvals;