import numpy as np
import matplotlib.pyplot as plt
p_oud = np.array([1,1,1])
p_nieuw = np.array([.8,.1,.1]) #de start proporties (P)
r = np.array([[0,1,2],[2,1,1],[1,2,0]]) #de uitbetalingsmatrix (R)
iteraties = 50
p1s = []
p2s = []
p3s = []

iteratie = 0
while iteratie < iteraties :
    iteratie += 1
    print("************************")
    print("iteratie:", iteratie)
    print("************************")
    s = np.dot(r,p_nieuw)
    q = np.multiply(p_nieuw,s)
    p_oud = p_nieuw
    p_nieuw = (1/sum(q))*q
    p1s.append(p_nieuw[0])
    p2s.append(p_nieuw[1])
    p3s.append(p_nieuw[2])
    print("S:", s)    
    print("Q:", q)
    print("oude P:", p_oud)
    print("nieuwe P:", p_nieuw)
    print("")

x = np.arange(iteraties)
plt.plot(x,p1s)
plt.plot(x,p2s)
plt.plot(x,p3s)
plt.show()