Problema 2

686 days ago by lop08413

RealNumber=float 
       
Integer=int 
       
from cvxopt.base import matrix 
       
from cvxopt import solvers 
       
c = matrix([0.75, 0.5,0.15]) 
       
A = matrix([[35,60,30], [0.5,300,20],[0.5,10,10]]) 
       
b = matrix([0.5,15,4]) 
       
sol = solvers.lp(-c,A,b) 
        
     pcost       dcost       gap    pres   dres   k/t
 0: -7.0511e-02 -7.0511e-02  2e-02  1e-01  3e-16  1e+00
 1: -7.0511e-02 -7.0511e-02  2e-04  1e-03  2e-16  1e-02
 2: -7.0511e-02 -7.0511e-02  2e-06  1e-05  2e-16  1e-04
 3: -7.0511e-02 -7.0511e-02  2e-08  1e-07  2e-16  1e-06
 4: -7.0511e-02 -7.0511e-02  2e-10  1e-09  2e-16  1e-08
     pcost       dcost       gap    pres   dres   k/t
 0: -7.0511e-02 -7.0511e-02  2e-02  1e-01  3e-16  1e+00
 1: -7.0511e-02 -7.0511e-02  2e-04  1e-03  2e-16  1e-02
 2: -7.0511e-02 -7.0511e-02  2e-06  1e-05  2e-16  1e-04
 3: -7.0511e-02 -7.0511e-02  2e-08  1e-07  2e-16  1e-06
 4: -7.0511e-02 -7.0511e-02  2e-10  1e-09  2e-16  1e-08
print sol['x'] 
        
   9.5263e-03
   3.8265e-02
   2.9489e-01
   9.5263e-03
   3.8265e-02
   2.9489e-01