O pacote lpSolveAPI
aceita restrições de três tipos, "<="
, ">="
e "="
. Creio que na pergunta a demanda deve ser entendida com demanda máxima, não como demanda fixa. Então as restrições a serem usadas serão de menoridade, "<="
#x1 = 50, pelo x2 = 50, pelo x3 = 80 e pelo x4 = 120
add.constraint(lprec, 1, "<=", 50, indices = 1) #restricao 4a
add.constraint(lprec, 1, "<=", 50, indices = 2) #restricao 4b
add.constraint(lprec, 1, "<=", 80, indices = 3) #restricao 4c
add.constraint(lprec, 1, "<=", 120, indices = 4) #restricao 4d
Agora resolve-se o problema tal como está na pergunta. As soluções serão dadas por get.variables(lprec)
.
solve(lprec) #resolver ppl
#[1] 0
get.objective(lprec) #funcao objetivo resolvida
#[1] 1948.62
get.variables(lprec) #qto de cada variavel otimiza o sistema
#[1] 2.6654343 0.5221156 0.0000000 0.2073184
Verificar estas soluções multiplicando os coeficientes da função objetivo pelo vetor solução.
t(c(600, 550, 400, 300)) %*% get.variables(lprec)
# [,1]
#[1,] 1948.62