clc
clear
b10 =1000*[4.6281;0.0058;0.3957;0.1203;0.0281;0.1763 ;0.0065;0.0169;0.0895]
Aeq=[50.7467 51.2556 46.6971 56.8828 58.3879 55.1321 0 0 0 0
0.0630 0.0636 0.0300 0.0097 0.0877 0.1122 0 0 0 0
2.9956 3.0257 4.6033 4.5767 5.5857 6.7706 0 0 0 0
0.0085 0.0086 0.0243 0.1496 0.6327 5.8034 0 0 0 0
0.0851 0.0860 0.0324 0.2200 0.8586 0.9672 0 0 0 0
1.7021 1.7191 2.2692 1.9803 0.7773 2.4181 0 0 0 0
0.0919 0.0928 0.0673 0.0194 0.0208 0.0648 0 0 0 0
0.1021 0.1031 0.1216 0.6601 0.1356 0.3095 0 0 0 0
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0 0 0 0]
lb=[5;5;3;3;3;19]
ub=[16;45;45;6;6;20]
f=[9.3300;10.3100;9.3300; 9.8400;5.5000;9.820;0;0;0;0]
[x,fval,exitflag,output]=linprog(f,[],[],Aeq,b10,lb,ub)