(%i4) eq1 : 8 * x - 2 * y - 1 * z = 1 $
eq2 : 9 * x - 4 * y + b * z = - 2 $
eq3 : 2 * x + 1 * y + a * z = 2 $
sol : linsolve([ eq1, eq2, eq3],[ x, y, z]) ;
(sol) [ x = 5 b + 8 a + 6 12 b + 14 a + 17 , y = 14 b + 25 a + 22 12 b + 14 a + 17 , z = 13 12 b + 14 a + 17 ]

Como la solución del sistema es [1/3,1,-1/3]:

(%i5) linsolve( ev( sol,[ x = 1 / 3, y = 1, z = - 1 / 3]),[ a, b]) ;
solve: dependent equations eliminated: (3) (%o5) [ a = 1 , b = 3 ]

Luego la suma de a+b es

(%i6) ev( a + b, %) ;
(%o6) 2

Created with wxMaxima.