calcs creates the decaying rally values till the rally falls below 0.005 (i.e. too small to be significant)

plot.m

vals = calcs(5);
hours = size(vals, 1);

x_vals = [1:hours]

plot(x_vals, vals)
xlabel('hours')
ylabel('rally')
title('Rally over time')

function vals = calcs(start)
    
    vec = eye(1,1);
    vec(1,1) = start;
    
    i = 0;
    while (start > 0.005)
        i = i + 1;
        if (start > 1)
            start = start - start*0.04;
        else
            start = start - start*0.02;
        end
        vec(i,1) = start;
    end
    vals = vec;
    
end