For comparison with other CAS choose from: Derive Macsyma Maple Mathematica Reduce
polys := [45*p + 35*s - 165*b - 36, _
35*p + 40*z + 25*t - 27*s, _
15*w + 25*p*s + 30*z - 18*t - 165*b**2, _
- 9*w + 15*p*t + 20*z*s, _
w*p + 2*z*t - 11*b**3, _
99*w - 11*s*b + 3*b**2, _
b**2 + 33/50*b + 2673/10000]
[35s + 45p - 165b - 36, 40z + 25t - 27s + 35p, 2 3 30z + 15w - 18t + 25p s - 165b , 20s z - 9w + 15p t, 2t z + p w - 11b , 2 2 33 2673 99w - 11b s + 3b , b + -- b + -----] 50 10000 Type: List Polynomial Fraction Integer
vars := [w, p, z, t, s, b]
[w,p,z,t,s,b] Type: List Symbol
groebner(polys)
49 1143 19 1323 37 27 5 9 [z + -- b + ----, w + --- b + -----, t - -- b + ---, s - - b - ---, 36 2000 120 20000 15 250 2 200 31 153 2 33 2673 p - -- b - ---, b + -- b + -----] 18 200 50 10000 Type: List Polynomial Fraction Integer
solve(polys, vars)
[ - 9500b - 3969 3100b + 1377 - 24500b - 10287 1850b - 81 [w= --------------, p= ------------, z= ----------------, t= ----------, 60000 1800 18000 750 500b + 9 2 s= --------, 10000b + 6600b + 2673= 0] 200 ] Type: List List Equation Fraction Polynomial Integer