Bæreevne-bestand-1

import matplotlib.pyplot as plt
import numpy as np

Reproduksjon = 200      #Hvor mange prosent populasjonen øker hvert år
Tid = 50                #Hvor mange år som simuleres
Sykdom = 50             #Hvor mange prosent som dør av sykdomsutbrudd
Rovdyr = 15             #Hvor mange prosent som spises av rovdyr
Miljøendring = 70       #Hvor mange prosent som dør ved en stor miljøendring
Mattilgang = 1000       #Hvor mange dyr det er mat til i økosystemet
Pop = 30                #Hvor mange dyr som er i området når simuleringen starter
Populasjon = [Pop]
År = [0]
for i in range (1,Tid):
    S = np.random.randint(1,5)  #Velger et tall mellom 1 og 5. Hvis 1 velges, vil sykdom 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*(Sykdom/100))*(1-M*(Miljøendring/100))
    if Pop > 1000:
        Pop = Pop*0.1+Mattilgang    #Hvis Pop>Mattilgang, vil de 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()