<<math4ti2`

sys={x-y≤2,-3 x+y≤1,x+y≥1,y≥0}

{x-y≤2,-3 x+y≤1,x+y≥1,y≥0}

zsolve[sys]