Saturation between lowerlimit and upperlimit
[sim,y] = ld_sat(sim, events, inp_list, lowerlimit, upperlimit)