I need to solve a mix of equations/inequalities:
a1 + a2 + a3 >= 50
b1 + b2 + b3 <= 40
c1 + c2 + c3 <= 10
a1 + b1 + c1 = 50
a2 + b2 + c2 = 35
a3 + b3 + c3 = 15
There are nine variables and six equations/inequalities, but I will be zeroing three of the variables at a time (I will have to check all combinations).
I tried a couple of packages in R (limSolve, matlib) with no success.
#The matrix representation:
X1 <- c(1,1,1,0,0,0,0,0,0)
X2 <- c(0,0,0,1,1,1,0,0,0)
X3 <- c(0,0,0,0,0,0,1,1,1)
Y1 <- c(1,0,0,1,0,0,1,0,0)
Y2 <- c(0,1,0,0,1,0,0,1,0)
Y3 <- c(0,0,1,0,0,1,0,0,1)
A1 <- matrix(c(X1,X2,X3,Y1,Y2,Y3),c(9,9))
A2 <- t(A1); colnames(A2) <- c("a1","a2","a3","b1","b2","b3","c1","c2","c3")
b <- c(50,40,10,50,35,15)
Any help will be appreciated.
This can be solved using linear programming. There is nothing to optimize, you are just searching for a feasible solution.
library(lpSolve)
res=lp(
"min",
rep(0,ncol(A2)),
A2,
c(">=","<=","<=","==","==","=="),
b
)
res$solution
[1] 50 35 15 0 0 0 0 0 0