Rabbits and Foxes
unknown
python
a year ago
2.3 kB
6
Indexable
# Number of rabbits and foxes r_initial = 500 f_initial = 1 alpha = 0.2 beta = 0.005 gamma = 0.001 delta = 0.2 r_t = r_initial f_t = f_initial for t in range(0, 15): print("At time t =", str(round(t)) + ",", "there are", round(r_t), "rabbits, and", round(f_t), "foxes") r_next = r_t + (alpha * r_t) - (beta * r_t * f_t) f_next = f_t + (gamma * f_t * r_t) - (delta * f_t) r_t = r_next f_t = f_next if r_t < 0: r_t = 0 if f_t < 0: f_t = 0 # How long can 100 rabbits survive? max_years = 0 max_f_initial = 0 for f_initial in range(1, 1001): r_initial = 100 alpha = 0.2 beta = 0.005 gamma = 0.001 delta = 0.2 t = 0 r_t = r_initial f_t = f_initial while r_t != 0: r_next = r_t + (alpha * r_t) - (beta * r_t * f_t) f_next = f_t + (gamma * f_t * r_t) - (delta * f_t) r_t = r_next f_t = f_next if r_t < 0: r_t = 0 if f_t < 0: f_t = 0 t = t + 1 years_survived = t - 1 if years_survived > max_years: max_years = years_survived max_f_initial = f_initial print("If there are", max_f_initial, "foxes at t = 0, the rabbit species will survive for", max_years, "years, which is the maximum possible number of years the rabbits can survive") # Saving rabbits from extinction r_initial = 1 found = False r_equilibrium = 0 f_equilibrium = 0 while True: for f_initial in range(1, 201): alpha = 0.2 beta = 0.005 gamma = 0.001 delta = 0.2 r_t = r_initial f_t = f_initial while r_t != 0: r_next = r_t + (alpha * r_t) - (beta * r_t * f_t) f_next = f_t + (gamma * f_t * r_t) - (delta * f_t) # Break out of for loop since either rabbits or foxes went extinct if r_t <= 0 or f_t <= 0: break if r_t == r_next and f_t == f_next: found = True r_equilibrium = r_initial f_equilibrium = f_initial break r_t = r_next f_t = f_next if found: break r_initial = r_initial + 1 print("For", f_equilibrium, "foxes, and", r_equilibrium, "rabbits at t = 0, the fox and rabbit species will never be extinct.")
Editor is loading...
Leave a Comment