Well, your code matches the results from Matlab/Simulink very closely… But it does not look like a discrete time implementation…
How can I implement the discrete time transfer function:
\frac{1}{T} \frac{z-1}{z+T_s/T-1}
with T being the time constant and T_s the sampling time in Julia directly?