Implements a time discrete integrator with saturation of the output between min__ and max__
[sim, y] = ld_limited_integrator(sim, ev, u, min__, max__, Ta)
u * - input
y * - output
y(k+1) = sat( y(k) + Ta*u , min__, max__ )