Definamos la función:

(%i2) f( v) : = matrix([ v[ 1] + v[ 2] - 2 * v[ 4], v[ 1] - 3 * v[ 4]],
   [ v[ 2] + 2 * v[ 3] - 2 * v[ 4], - v[ 1] + 2 * v[ 3] - v[ 2] - v[ 4]]) $
f([ a, b, c, d]) ;
(%o2) [ 2 d + b + a a 3 d 2 d + 2 c + b d + 2 c b a ]

Construimos el isomorfismo que nos lleva las matrices 2x2 a vectores de cuatro componentes:

(%i3) define( Isom( m), flatten( makelist( makelist( m[ i, j], j, 1, 2), i, 1, 2))) $

El núcleo estará compuestos por los vectores cuya imagen es la matriz cero:

(%i5) eq : Isom( f([ x, y, z, t])) $
print( "Ker f={", makelist( eq[ i] = 0, i, 1, 4), "}") $
Ker f={ [ y + x 2 t = 0 , x 3 t = 0 , 2 z + y 2 t = 0 , 2 z y x t = 0 ] }

Por tanto, el núcleo será el subespacio vectorial solución del sistema:

(%i6) sol : linsolve( eq,[ x, y, z, t]) ;
solve: dependent equations eliminated: (4) (sol) [ x = 3 %r1 , y = %r1 , z = 3 %r1 2 , t = %r1 ]

Como tenemos un parámetro el subespacio está generado por un solo vector:

(%i8) v : ev([ x, y, z, t], ev( sol, %rnum_list[ 1] = 1)) $
v / sqrt( v. v) ; /* lo normalizamos */
(%o8) [ 6 53 , 2 53 , 3 53 , 2 53 ]

La respuesta a la pregunta es:

(%i9) round( abs( sum( %[ i], i, 1, 4)) * 100) / 100, numer ;
(%o9) 1.24

Created with wxMaxima.