%%Q2_mk3
profit=[80 50 40];
components =   [1    1   0
                1    0   0
                1    1   0
                1    1   0
                2    2   1];
max_profit=0;
another=[];
P=[];
p1lower=0;p1inc=5;p1upper=250;
p2lower=0;p2inc=5;p2upper=450;
p3lower=0;p3inc=5;p3upper=600;
%For-loop for finding maximum profit
for p1=p1lower:p1inc:p1upper
    for p2=p2lower:p2inc:p2upper
        for p3=p3lower:p3inc:p3upper
            p=[p1,p2,p3];
            if components*p<=[450,250,800,450,600]
                profit=unit_profit*p;
            end
            if profit>max_profit
                max_profit=profit;
                production=[p1,p2,p3];
            end
        end
    end
end
disp('The maximum profit is:')
disp(max_profit)