# Primero me defino las variables
# Ji: cantidad de manteles usados de los comprados en dia i
# K1i: cantidad de manteles usados de los lavados en 1 dia del dia i
# K2i: cantidad de manteles usados de los lavados en 2 dias del dia i

# Podria haber definido variables totalizadoras pero se me hacia mas largo
# vector = (J1,J2,J3,J4,J5,J6,J7,K11,K12,K13,K14,K15,K21,K22,K23,K24)

obj <- c(250,250,250,250,250,250,250,150,150,150,150,150,100,100,100,100)
res <- matrix(c(1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,
                0,0,0,1,0,0,0,0,1,0,0,0,1,0,0,0,
                0,0,0,0,1,0,0,0,0,1,0,0,0,1,0,0,
                0,0,0,0,0,1,0,0,0,0,1,0,0,0,1,0,
                0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,1,
                0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,
                0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,
                0,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,
                0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,1,
                0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0), nrow=12, byrow=TRUE)
sig <- c("=","=","=","=","=","=","=","<=","<=","<=","<=","<=")
rhs <- c(5,6,7,8,7,9,10,5,6,7,8,7)

lpSolve::lp("min", obj, res, sig, rhs)

lpSolve::lp("min", obj, res, sig, rhs)$Solution

# asi como esta no me devuelve los X. Re edito

Z <- lpSolve::lp("min", obj, res, sig, rhs)
X <- Z$solution

print(Z)
print(X)

# Z=8800
# X: 5  6  7  3  0  0  0  0  0  1  3  5   5   6  6   5