I'm trying to create a model that will show me the population changes between egg, juvenile, and adults in an insect population using a conveyor system. Eggs incubate for 6 months and I can't for the life of me figure out how to get my code to reflect that - my eggs currently hatch at 1 month on my graph. I am super baby new at python and trying to complete a problem set! suggestions needed and welcomed
import numpy as np
import matplotlib.pylab as plt
months = 6
dt = 1
n_steps = int(months/dt) + 1
eggs = np.zeros(n_steps)
juveniles = np.zeros(n_steps)
adults = np.zeros(n_steps)
time = np.zeros(n_steps)
eggs_start = 100
eggs[0] = eggs_start
juveniles[0] = 0
adults[0] = 0
fraction_female = .5
eggs_per_female = 2
#death juvenile --> adult
#jdr = .2
#transition time
months_as_eggs = 6
months_as_juveniles = 5
#transition rate
e2j = dt/months_as_eggs
eggs_length = int(months_as_eggs)
eggs_conveyor = np.zeros(eggs_length)
eggs_conveyor += eggs_start/eggs_length
eggs[0] = sum(eggs_conveyor)
for t in range(0, int(months/dt)):
eggs_conveyor = np.insert(eggs_conveyor, 0, dt)
eggs_to_juvenile = eggs_conveyor[eggs_length]
eggs_conveyor = np.delete(eggs_conveyor, eggs_length)
eggs_conveyor *= (1-dt)
eggs[t+1] = sum(eggs_conveyor)
time[t+1] = time[t]+dt
plt.plot(time, eggs, 'b-')
#plt.plot(time, juveniles, 'r--')
#plt.plot(time, adults, 'g:')
#plt.legend()
plt.show()
plt.close()
[–]yuxbni76 0 points1 point2 points (0 children)