klimavektor = [randi([18,28],365,1), zeros(365,1)];
[zk, sk] = size(klimavektor);
bakterium = 0;
tot = 3;
for k=1:1:zk
    if klimavektor(k,1) >=20
        bakterium = bakterium + 1;
        tot = 3;
    else
        tot = tot - 1;
        if tot == 0
            bakterium = 0;
        end
    end
    klimavektor(k,2) = bakterium;
end

plot(klimavektor)
axis([0 365 0 40])