Erratum
Eq.(23) should read y(n)=f_1 y(n-1) + f_2 x(n-1) + f_3 x(n)
This was already correctly implemented in the Fortran code provided as Supplementary material.