version 1 of get matrix function

getMatrix.py

>>> import numpy as np
>>> def getMatrix(matrix, turns):
...     newMatrix = np.identity(matrix.shape[0])
...     for i in range(turns):
...             newMatrix = newMatrix*matrix
...     return newMatrix

>>> twoTurnMatrix = getMatrix(transitionMatrix,2)
>>> twoTurnMatrix
matrix([
        [2, 1, 0, 1, 0, 0, 0, 1, 0, 1],
        [1, 2, 0, 1, 0, 0, 0, 1, 0, 0],
        [0, 0, 2, 0, 1, 0, 1, 0, 0, 0],
        [1, 1, 0, 2, 0, 0, 0, 0, 0, 1],
        [0, 0, 1, 0, 3, 0, 1, 0, 1, 0],
        [0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
        [0, 0, 1, 0, 1, 0, 3, 0, 1, 0],
        [1, 1, 0, 0, 0, 0, 0, 2, 0, 1],
        [0, 0, 0, 0, 1, 0, 1, 0, 2, 0],
        [1, 0, 0, 1, 0, 0, 0, 1, 0, 2]])