matinbasis:=function(mat,basi)
   local m,x;
   m:=List(basi,x->x*mat);
   Apply(m,x->SolutionMat(basi,x));
   return m;
end;

# matinbasis([[1,0],[0,2]],[[1,1],[0,2]]);