% Grafica i risultati
close all

th = 0.1;

maxr = max(r(find(r(:,1)<ts/2),2));
minr = min(r(find(r(:,1)<ts/2),2));

figure, plot(r(:,1),r(:,2),'-',...
        r(:,1),maxr*(1+th*sign(maxr))*ones(size(r(:,1))),'--',...
        r(:,1),minr*(1-th*sign(minr))*ones(size(r(:,1))),'--'), 
        title('Residuo con guasto su y(t), f_y(t)'), xlabel('Tempo (s.)')