For a set of under-determined linear equations, I was wondering if there is any closed form for all non-negative solutions? Is there a way to analytically characterize the feasibility set of such equations?
More formally, for given $M\in \mathbb{R}^{n\times m}_{+}$ and $y\in \mathbb{R}^{m\times 1}$, find all solution $x\geq 0$ of $Mx=y$. Note that $M$ and $y$ are non-negative as well.