For woking with inequalities, Maxima also has code for the simplex algorithm. I don't know a great deal about this package; here are two examples: (%i7) load(simplex)$ (%i13) minimize_lp(x,[1 < x, x < y, y < 23]); (%o13) [1,[y=23,x=1]] (%i14) maximize_lp(x,[1 < x, x < y, y < 23]); (%o14) [23,[y=23,x=23]] --Barton