%%
clc; clear all; close all;

%%

f= [-2 -4 -5.2 0.3*60 0.3*40 0.3*10 0.4*60 0.4*40 0.4*10 0.3*60 0.3*40 0.3*10].*-1;
a = [-1 0 0 8 6 1 0 0 0 0 0 0;
     0 -1 0 4 2 1.5 0 0 0 0 0 0; 
     0 0 -1 2 1.5 0.5 0 0 0 0 0 0;
     -1 0 0 0 0 0 8 6 1 0 0 0 ;
     0 -1 0 0 0 0 4 2 1.5 0 0 0; 
     0 0 -1 0 0 0 2 1.5 0.5 0 0 0;
     -1 0 0 0 0 0 0 0 0 8 6 1 ;
     0 -1 0 0 0 0 0 0 0 4 2 1.5 ; 
     0 0 -1 0 0 0 0 0 0 2 1.5 0.5;
    ];
b= zeros(9,1);
lb = zeros(12,1);
ub = [inf inf inf 50 20 200 150 110 225 250 250 500]';

intcon = [1 2 3 4 5 6 7 8 9 10 11 12]';

%%

[sol fval] = intlinprog(f,intcon,a,b,[],[],lb,ub);