num_runs = 100;
delivery_quantity = 4000; %number to have delivered in each shipment (X in notes)


for j=1:6
    for i = 1:25
    delivery_interval = i; 
    DAILY_COST(:,i,j) = simulate(num_runs, delivery_interval, delivery_quantity);
    end
    delivery_quantity = delivery_quantity + 4000
end

for j=1:6
    for i = 1:25
    AVES(i,j) = sum(DAILY_COST(:,i,j))/num_runs
    end
end



semilogy(1:25,AVES(:,1),'.',1:25,AVES(:,2),'x',1:25,AVES(:,3),'o',1:25,AVES(:,4),'+',1:25,AVES(:,5),'v',1:25,AVES(:,6),'^')
legend('delivery quantity 4000','delivery quantity 8000','delivery quantity 12000','delivery quantity 16000','delivery quantity 20000','delivery quantity 24000')
ylabel('average daily cost')
xlabel('delivery interval in days')