Bæreevne-bestand2-NN

import matplotlib.pyplot as plt
import numpy as np

Reproduksjon = 200      #Kor mange prosent populasjonen aukar kvart år
Tid = 50                #Kor mange år som blir simulerte
Sykdom = 50             #Kor mange prosent som døyr av sjukdomsutbrot
Rovdyr = 15             #Kor mange prosent som blir etne av rovdyr
Miljøendring = 70       #Kor mange prosent som døyr ved ei stor miljøendring
Mattilgang = 1000       #Kor mange dyr det er mat til i økosystemet
Pop = 30                #Kor mange dyr som er i området når simuleringa startar
Populasjon = [Pop]
År = [0]
for i in range (1,Tid):
    S = np.random.randint(1,5)  #Vel eit tal mellom 1 og 5. Viss 1 blir valt,vil sjukdom inntreffe (20 % sjanse).
    if S == 1:
        S = 1
    else:
        S = 0
    M = np.random.randint(1,15)
    if M == 1:
        M = 1
    else:
        M = 0
    Pop = Pop*(Reproduksjon/100)*((100-Rovdyr)/100)*(1-S*(Sjukdom/100))*(1-M*(Miljøendring/100))
    if Pop > 1000:
        Pop = Pop*0.1+Mattilgang    #Viss Pop>Mattilgang, vil dei 1000 det er nok mat til + 10 % av bestanden overleve.
    Populasjon.append(Pop)
    År.append(i)

plt.plot(År,Populasjon)
plt.xlabel("År")
plt.ylabel("Populasjon")
plt.show()