Euler's Method in Python: Predator-Prey Example

Hi everyone! This video is about how to implement Euler's method for numerical ODE solving in Python, using the Lotka-Volterra Predator-Prey model as an example.
All the code from my videos is available on my Github:
github.com/mikesaint-antoine/...
Thanks for watching and let me know if you have any questions!

[@MikeSaintAntoine]

Hi everyone, thanks for watching! Here's the link to the code on my Github in case you want to try it for yourself (also feel free to add me as a friend on Github if you want 🙂) : github.com/mikesaint-antoine/Comp_Bio_Tutorials/blob/main/numerical_methods/eulers_method_predator_prey.py

[@shabihah]

Thank you for your explanation. It really helps me to understand my uni module!

[@MikeSaintAntoine]

Thanks for watching! 🙂

[@achrafBadiry]

I had problems understanding the Euler numerical approx. for a "prey-predator" uni. programming project. but you explained perfectly for this matter. thank you !

[@MikeSaintAntoine]

Thanks for watching! 🙂

[@anikamittal6721]

Hi Mike, I have no experience coding, but I wanted to calculate the parameters to create a Prey-Predator model for data I have collected. Is there any simple way to do so using python?

[@MikeSaintAntoine]

Hi Anika, I think that will be a bit difficult with no experience coding, but maybe I can try to help. Can you email me at [e-mail] and tell me about what you're trying to do, and attach the dataset? Then I'll see if there's an easy way to do it. Thanks!

[@cameronl1434]

Hi mike, I am running the simulation with the same variables as you, with a slightly modified code, but I am unable to figure out why when I plot the results, the prey population immediately dumps to 0 rather than increases? here is my code below, any pointers would be greatly appreciated # model 1 - exponential prey growth with type I functional response while x > 0 and y > 0 and t < tmax: dx = expo_growth(alpha, x) - f_type1(beta, x, y) dy = pred_growth_t1(delta, x, y) - pred_death(gamma, y) #update values x += dx * dt y += dy * dt t += dt # append values to lists prey_pop_m1.append(x) pred_pop_m1.append(y) t_m1.append(t)

[@MikeSaintAntoine]

Hi Cameron, yeah I can try to help with that! Could you email me your full code at [e-mail]? Then I'll try to run it and take a look. It's hard to tell based on the short code snippet because the problem might be in one of the functions that you're using for the different ODE terms.

[@cameronl1434]

Thank you so much, I've popped over an email