module: resonant_bandpass
parameters: double fo, double q
inputs: double in
outputs:  double out
static_variables: double wo_warp
classes: Filter filt("wo/Q*s","wo*wo + wo/Q*s + s^2","wo,Q,Ts",2*pi*fo,q,Ts)
init:
//calculate pre-warped center frequency
wo_warp = (2/Ts)*tan(2*pi*fo*Ts/2);
filt.set("wo/Q*s","wo*wo + wo/Q*s + s^2","wo,Q,Ts",wo_warp,q,Ts); 
code:
filt.inp(in);
out = filt.out;
