function show_dig(m0,time,x,s)
for i=1:size(m0,1)
    dig(i,:) = m0(i,i,:);
end
shadedErrorBar(time',mean(dig,2),std(dig,0,2)./sqrt(size(dig,2)),s,0.5);
ylim([-x, x]);
end