function res = trivial(alpha, epsilon, T, draw) x(1) = alpha; c(1) = 0; for k=1:T-1 x(k+1) = x(k) + epsilon * ( round(x(k) - c(k)) - x(k)); c(k+1) = c(k) + round(x(k)- c(k)) - x(k); end res.x = x; res.c = c; res.min_x = min(x); res.max_x = max(x); res.dev = max(abs(x - alpha)); if draw figure subplot(2,1,1); plot(x,'.'); subplot(2,1,2); plot(c,'.'); end