Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement