import numpy as np

A = np.array([[0.86,0.08],
             [-0.12,1.14]])

# Antall ?r
n = 100

# Antall ulver og antall harer
x0 = np.array([100,51])
x = x0
print('?r ' + str(0) +  ': ' + str(np.rint(x[0])) + ' ulver og ' + str(np.rint(x[1])) + ' harer.')
h = True
u = True

# Dynamisk system
for j in range(1,n+1):
  x = np.matmul(A,x)
  if (x[0] >= 1 and x[1] >= 1):
     print('?r ' + str(j) +  ': ' + str(np.rint(x[0])) + ' ulver og ' + str(np.rint(x[1])) + ' harer.')
  if (x[1] < 1 and h):
    h = False
    print('Harene ble spist opp i ?r ' + str(j) + '. :-(')
  if (x[0] < 1 and u):
    u = False
    print('Ulvene d?de ut i ?r ' + str(j) + '. :-(')
	
# Uttrykk x0 ved hjelp av egenvektor-basisen
B = np.array([[2,1],
              [1,3]])
a = (np.linalg.solve(B,x0))[0]
b = (np.linalg.solve(B,x0))[1]

print("Her er a=" + str(a) + " og b=" + str(b) + ".")