(%i3) eq1 : 3 * x - y - z = 2 $
eq2 : x + y - 2 * z = 1 $
sol : linsolve([ eq1, eq2],[ x, y, z]) ;
(sol) [ x = 3 %r1 + 3 4 , y = 5 %r1 + 1 4 , z = %r1 ]

Observemos que nos aparecen los parámetros %rn, estos los tenemos en %rnum_list, así podremos hacer mención de ellos aunque estos cambien en cada ejecución.

Expresemos las soluciones en forma matricial.
Un punto de la variedad lineal solución del sistema sería el dado con los parámetros igual a cero:

(%i4) P : ev([ x, y, z], ev( sol, %rnum_list[ 1] = 0)) ;
(P) [ 3 4 , 1 4 , 0 ]

Una base del espacio director está formada por el vector:

(%i5) v : ev([ x, y, z], ev( sol, %rnum_list[ 1] = 1)) - P ;
(v) [ 3 4 , 5 4 , 1 ]

La forma matricial es: (utilizo una representación del espacio director con vector l.d. que evite las fracciones)

(%i6) print( matrix([ x],[ y],[ z]), "=", transpose( matrix( P)), "+ λ",
   transpose( 4 * matrix( v))) $ ;
[ x y z ] = [ 3 4 1 4 0 ] + λ [ 3 5 4 ]

Nos preguntan por la suma, en valor absoluto, de las coordenadas del vector director unitario de la recta afín; es decir,

(%i8) v / sqrt( v. v) ;
round( abs( sum( %[ i], i, 1, 3)) * 100) / 100, numer ;
(%o7) [ 3 5 2 , 1 2 , 2 3 2 5 ] (%o8) 1.7

Created with wxMaxima.