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

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

Step 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_srid(0.01,1,1.2,5);