overimpose_the_structure(premium, XOL, loss, lossPick, commission)

function loss = overimpose_the_structure(premium, XOL, loss, lossPick, commission)

%Adjusting for Excess of Loss contract
loss(find(loss > XOL)) = XOL;

%Each column represents one simulation
loss = sum(loss);

%Adjusting losses for commission - 5% for loss less than losspick

for i=1:length(loss)
    if loss(i) < lossPick
        loss(i) = loss(i) + commission*premium;
    end
end

%Looking at losses after overimposing the structure
figure;
subplot(2,1,1); hist(loss);
title('Distribution of losses after adjustments');
xlabel('Dollars ($)');
ylabel('#of observations');

subplot(2,1,2); plot(loss);
xlabel('Multiple Simulations');
ylabel('Dollars ($)');