R中不等式组的解

我有以下不等式系统：Ay>=0，其中A是9x3矩阵，

y=（y1，y2，y3）

是3个元素的向量。不等式的解是一个区域，但我想返回一个可能的元组（y1，y2，y3）来解这个不等式。请注意，y的所有元素都必须大于0

有没有办法在R中实现这一点

谢谢

 A = matrix(runif(27, -0.5, 0.5), nc = 3)

 ### The system to solve is:    A %*% y >= as.matrix(rep(0,9))    Solve for y. Do not return region, but rather any possible tuple that solves the inequalities

---

试试线性规划：

library(Rglpk)
rhs <- rep(c(0, 1e-3), c(9, 3))
ge <- rep(">=", 12)

# an example with no feasible solution (status = 1)
set.seed(123)
A <- rbind(matrix(runif(27, -0.5, 0.5), nc = 3), diag(3))
Rglpk_solve_LP(obj = numeric(3), mat = A, dir = ge, rhs = rhs)

# an example with a feasible solution (status = 0)
A2 <- rbind(cbind(numeric(9), 1, -1), diag(3))
Rglpk_solve_LP(obj = numeric(3), mat = A2, dir = ge, rhs = rhs)

库（Rglpk）
rhs
A%*%y
是一个向量。如何确定“>0”？
optim
的问题是，由于我不知道可行区域是什么，如何提供落在可行区域内的起始值？@Carl%*%y的每一行都必须大于零。我已经在帖子中调整了这个问题。。尺寸也不正确，我也重新调整了。
c（0，0，0）
是问题中所述问题的解决方案，与
a
无关。很抱歉，我没有指定y>0谢谢！请告诉我如何安装此软件包好吗？
install.packages（“Rglpk”）
我收到以下消息：将软件包安装到“C:/Users”（由于“lib”未指定）软件包“Rglpk”可用作源软件包，但不是二进制警告消息：软件包“Rglpk”不可用（对于R版本2.15.3）使用最新版本的R。你的已经过时多年了。好的，明白了！我会试试的。谢谢！