Compute the DC gain of the filter (assuming 0 if unstable) and use that to set the initial conditions of a filter as follows: x(t<0)=x(0) y(t<0)=x(0) × dc_gain