def concatenate_matrices(*matrices):
M = numpy.identity(<NUM_LIT:4>)<EOL>for i in matrices:<EOL><INDENT>M = numpy.dot(M, i)<EOL><DEDENT>return M<EOL>
Return concatenation of series of transformation matrices.

    >>> M = numpy.random.rand(16).reshape((4, 4)) - 0.5
    >>> numpy.allclose(M, concatenate_matrices(M))
    True
    >>> numpy.allclose(numpy.dot(M, M.T), concatenate_matrices(M, M.T))
    True