OD mana

clear all
close all
maxpoolarr = [600:200:1400];
maxhits = 100;
numsims = 2000;
cdf=zeros(numel(maxpoolarr),maxhits);
for j = 1:numel(maxpoolarr)
    maxpool=maxpoolarr(j);
    endval = zeros(1,numsims);
    for i = 1:numsims
        pool=maxpool;
        for endval(i) = 1:maxhits
            pool=pool-100;
            if (pool<100)
                break
            end
            if (rand < 0.4)
                pool=pool+maxpool*0.25;
            end
            if (pool>maxpool)
                pool=maxpool;
            end
        end
    end
    pdf = hist(endval,1:maxhits,100);
    cdf(j,:) = cumsum(pdf);
end
figh=figure
set(figh,'Position',[100, 100, 1280,800])
hold all
grid on
grid minor on
for i = 1:numel(maxpoolarr)
    plot(cdf(i,(1:maxhits-1)));
end
legend(strcat(num2str(maxpoolarr'),' mana'))
title('CDF for running out of mana on OD');
xlabel('Maximum amount of consecutive hits')
ylabel('Chance you are out of mana (%)')