#Objetivo : maximizar beneficio = ingreso - costo 
#XE: litros de whisky escoses
#XK:litros de whisky kilt
#XT:litros de whisky tartán
#XAi: litros de licor A destinados al whisky i
#XBi:litros de licor B destinados al whisky i
#XCi:litros de licor C destinados al whisky i
#i = E,K,T

library(lpSolve)

#Vector Objetivo 

OBJ<-c(6.8,5.7,4.5,-7,-7,-7,-5,-5,-5,-4,-4,-4)

#Matriz restricciones

REST<-c(-1,0,0,0,0,0,-0.60,-0.20,0,0,0,
        0,0,-1,0,0,0,0,0,-0.15,-0.60,0,
        0,-1,0,0,0,0,0,0,0,0,-0.50,
        1,0,0,1,0,0,1,0,0,0,0,
        0,0,1,1,0,0,0,0,1,0,0,
        0,1,0,1,0,0,0,0,0,0,0,
        1,0,0,0,1,0,0,0,0,0,0,
        0,0,1,0,1,0,0,0,0,0,0,
        0,1,0,0,1,0,0,0,0,0,0,
        1,0,0,0,0,1,0,1,0,0,0,
        0,0,1,0,0,1,0,0,0,1,0,
        0,1,0,0,0,1,0,0,0,0,1)

MAT<-matrix(REST, nrow=11, ncol=12)
print(MAT)

DIR<-c("=","=","=","<=","<=","<=",">=","<=",">=","<=","<=")
print(DIR)

RHS<-c(0,0,0,2000,2500,1200,0,0,0,0,0)

#RESOLUCION

SOL<-lp(direction="max", objective.in = OBJ, const.mat = MAT, const.dir = DIR, const.rhs = RHS)
print(SOL)

VECTORSOL<-SOL$solution
print(VECTORSOL)