Impulse response invariant discretization of BICO (Bode's Ideal Cut-Off) transfer function

Impulse response invariant discretization of BICO (Bode's Ideal Cut-Off) transfer function

Impulse response invariant discretization of BICO (Bode's Ideal Cut-Off) transfer function
BICO TF= 1./((s./w0+sqrt((s./w0).^2+1)).^r);
r is a positive real number
IN:
Ts: The sampling period
w0: Cutoff frequency
r: order
q : Denominator and Numerator order of the approximated z-transfer function
TEST CODE:
BICO_irid(0.01,1,1.2,5);