mult:=proc(A::Matrix,B::Matrix)local i,j,k,q,n,C;C:=Matrix(RowDimension(A),ColumnDimension(B));for i from 1 to RowDimension(A) dofor j from 1 to ColumnDimension(B) dofor q from 1 to ColumnDimension(A) do