function Boros()
    Aeq=[0.1 0.5 0.3 0.2 0.4; 0.4 0.3 0.3 0.2 0.1; 0.2 0.02 0.1 0.05 0.05; 0.3 0.18 0.3 0.55 0.45]
    f=[45; 20;50;90; 70]
    beq=[0.3 ; 0.2 ;0.1;0.4]
    x=linprog(f,[],[],Aeq,beq,zeros(4,1))
    
    options = optimset('LargeScale', 'off', 'Simplex', 'on');
    x=linprog(f,[],[],Aeq,beq,zeros(4,1),ones(4,1),zeros(5,1),options)
    
