(%i4) A : matrix([ 1, 2],[ 0, - 1]) $
B : matrix([ 0, - 1],[ 1, 3]) $
AmB : mat_trace( transpose(( B - A)).( B - A)) $
print( "dist(", A, ",", B, ")=", sqrt( AmB)) $
0 errores, 0 advertencias dist( [ 1 2 0 1 ] , [ 0 1 1 3 ] )= 3 3 2
(%i5) round( % * 1000) / 1000, numer ;
(%o5) 5.196

Created with wxMaxima.