MatAlgebra( F, n )
returns the full matrix algebra of n by n matrices over the field F.
gap> a:= MatAlgebra( GF(2), 2 ); UnitalAlgebra( GF(2), [ [ [ Z(2)^0, 0*Z(2) ], [ 0*Z(2), 0*Z(2) ] ], [ [ 0*Z(2), Z(2)^0 ], [ Z(2)^0, 0*Z(2) ] ] ] ) gap> Size( a ); 16
GAP 3.4.4