(%i2) f( v) : =[ v[ 1], v[ 2] - v[ 1], v[ 3] - v[ 1], 2 * v[ 1] - v[ 2]] $
e3 : ident( 3) $

El ejemplo nos dice que la aplicación es:

(%i4) expr : rat( f([ a, b, c]). matrix([ 1],[ X],[ X ^ 2],[ X ^ 3]), X) $
print( "f(a,b,c)=", expr) $
f(a,b,c)= ( b + 2 a ) X 3 + ( c a ) X 2 + ( b a ) X + a

La matriz asociada será la proporcionada por las imágenes de la base canónica de \(\mathbb{R}^3\)

(%i5) transpose( matrix( f( row( e3, 1)[ 1]), f( row( e3, 2)[ 1]), f( row( e3, 3)[ 1]))) ;
(%o5) [ 1 0 0 1 1 0 1 0 1 2 1 0 ]

Verifiquemos el resultado:

(%i6) rat( transpose( %. matrix([ a],[ b],[ c])). matrix([ 1],[ X],[ X ^ 2],[ X ^ 3]), X) ;
(%o6)/R/ ( b + 2 a ) X 3 + ( c a ) X 2 + ( b a ) X + a

Created with wxMaxima.