c =；

& gt& gtb =；

& gt& gtVLB =；

& gt& gtVUB =[]；

& gt& gt[x，fval]=linprog(c，A，b，Aeq，beq，VLB，VUB)

Optimization terminated successfully.

x =

1.0e+004 *

5.75 13

1.5036

0.0000

4.5000

4.0000

3.0000

4.2487

0.0000

2.6 140

0.0000

0.0000

fval =

- 1.4375e+005

& gt& gt Because it is minimization and your topic is maximization, I changed the coefficient of your objective function to a negative value, so I only need to add a negative sign before the final function value, that is, fval= 1.4375e+005.