subroutine low_pass(alpha, y_in, y_out, y_out_last)

use bmad
implicit none

real(rp) y_out, y_in
real(rp) y_out_last
real(rp) alpha

    y_out = (1-alpha)*y_out_last+alpha * y_in

end