library(lpSolve)

# Ej 2.4

# Defino los coeficientes de la función objetivo
f.obj <- c(1000, 1500, 1500, 1800)

# Defino la matriz con los coeficientes tecnológicos
f.con <- matrix(c(5, 6, 0, 0,
                  0, 0, 4, 4,
                  1.6, 0, 0, 1.2,
                  0, 1.8, 1.8, 0,
                  0, 1, 1, 0), nrow = 5, byrow = TRUE)

# Defino los signos de las inecuaciones
f.dir <- c("<=", "<=", "<=", "<=", ">=")

# Defino los RHS
f.rhs <- c(80, 80, 20, 36, 10)

# Valor Final (z)
lp("max", f.obj, f.con, f.dir, f.rhs)

# Resultado de las variables
solution <- lp("max", f.obj, f.con, f.dir, f.rhs)$solution
print(solution)