I assume new y0 = y0 + x instead of new y0 = y0 + x · u with u the facilitation term. this because It results a very good and more realistic approach to study of the GABAA dynamics. than just one mono-exponential.