How to Solve Coupled Differential Equations ODEs in Python

I walk through how to use the scipy odeint method within Python to solve coupled Ordinary Differential Equations (ODEs) and plot the results using matplotlib. Link to Python script created in this video: github.com/vas...
Link to the Google Doc containing the problem statement: docs.google.co...
Content adapted from APMonitor.com: • Simulate Coupled Diffe...

[@IKMCDANIEL]

Vincent - you have real talent. You explained a tremendous amount of technical content in just 10min. Very well done and very helpful. Thank you so much for creating this content.

[@VincentStevenson]

Thank you for your kind words, I appreciate it!

[@ospectrum]

Thank you so much, I am currently working on a paper of population dynamics, and your video literally saved me a ton of time!

[@4amoats]

hey, just checking up on you

[@dotpio8246]

This was so helpful. Straightforward, straight to the point, and really easy to understand. Thank you Vincent!!

[@ammarshahzad9627]

Hello mate, I was searching R_2_score and found your channel. I saw your multi disciplinary videos and I am amazed by your content. You have mastery over alot of fields and have the gift of teaching. Subscribed instantly, hope you will gain the recognition you deserve.

[@chasefoxen]

Fantastic explanation, this helped me a ton with my numerical analysis homework. Straight to the point and clear!

[@pineapple5923]

Thank you so much. Better than my professor by a mile!

[@rajanalexander4949]

Concise and useful; thank you!
I would just comment/say that we're assigning each _vector element_ to a _variable_ @ 3:03. The ODEs themselves are represented by the dXdt assignments.

[@joaocruz9048]

Thank you, Vicent! I was with difficult on working with more than one initial conditions and this video is really helpful.

[@schrodingergaming7]

That's a very good solution in a small time frame! Great job!

[@matts2565]

This video saved my life, thanks boss!!

[@ibrahimaniakh4447]

JUST EXCELLENTLY EXPLAINED. EXCELLENT, EXCELLENT, EXCELLENT. I’m immediately subscribing bruh.

[@user-uu4wg6xq9l]

Thank you for this video, it was extremely well explained. It was incredibly useful to me.

[@yogeshwarthakare8356]

Hi, this was so useful. Would like to see more such videos on solving kinetic model using Python

[@hasatamashi]

Great video, great voice, really helpful

[@michaelhanson5336]

Very helpful. Thank you!

[@animeshsinghal3405]

Thanks a lot. You got me started with this very quickly.

[@fannumber122]

Perfect explanation thank very much.

[@thickeyebrows1981]

Wonderful video!

[@trumpetjonny5783]

Nice, thank you so much 🙂

[@Onlinefizik]

Very good,Thanks

[@saswatipal7280]

Awesome! thank you. you saved me a big time.

[@anandamoymukhopadhyay5158]

Vincent you have explained the code nicely within a short span of time. If the above system contains some arbitrary parameters, then please explain how to deal with it by the help of continuation method? Waiting for your reply.

[@Novoselova.Polina]

thank you! really helpful

[@zishanshaikh3029]

Very helpful , mate , cheers !!

[@omurwaAfrica]

This is helpful. Thank you.

[@edholanda5897]

Professor how would I solve this system of first order edos numerically by plotting the graph for the different values ​​of (n). the derivatives are in relation to ha (r).
a'/r = -e^2*v^2*(g^2 - 1)
g' = - a*g/r
given the boundary conditions
a(0) = n a(inf)=0
g(0) = 0 g(inf)=1
o (n) varies from 1 to 8.
where (e)=0.5 and (v)=1 are constant. please give a helping hand there, I looked for and did not find any problems like this on the python website.
I'm from Brazil.

[@rodrigocesarraimundo3316]

Amazing!

[@XboxGamersAndMore]

Great vid, very helpful, thank you

[@carolina2066]

this is very nice and helpful thanks a lot :))

[@Elpidio1000]

Thank you so much dude.

[@mohamedvadhelebnououmar8779]

Thanks, very helpful

[@user-ng2jh9hv1u]

Hi Vincent, thank you for such a nice video, it's extremely useful. I was wondering if you can make a video to calculate the Lyapunov exponents of coupled nonlinear ODEs.
Thank you

[@jowadulkader9006]

thank you sir

[@minhtu5059]

Thank you so much!!!

[@ponho9715]

Thanks :D

[@raulfuerte5512]

Hi Vincent. This was an amazing demo of using python's IVP solver. For a school project I was wondering if you could make a video showing how to use python to solve BVPs with a combination of Neumann and Dirichlet boundary conditions.

[@mehrdadkazemi3969]

Thank you so much

[@lordyabo9910]

guys I tried to copy his code word for word and run the code, but I got nothing. Can someone give me advice on what should I do. Btw im using pydroid 3

[@bryanlay5215]

What if there’s an error saying list object not callable?

[@VincentStevenson]

I would double check to make sure when you create the class object, you have the () in place before calling the object after.

[@MsJennifervelazquez5]

thanks, you awesome

[@AF-ge4pe]

what if they are coupled odes but with parameters that im asked to define with runge kutta 4th order

[@h.a.3030]

First thank you for this simple well explained video,i'm actually working on the same model dynamics, i am wondering if you can help me for example to see the evolution of just one variable with respect to a changing parameter, i tried to use a loops, but i can't get to the results I'm expecting! TIA for any help.

[@lonemeertens7021]

Hi Vincent, thank you for this video. I have a follow up question :) What if I have (in one of the equations) a parameter (instead of a constant) that is linked to a algebraic equation that itself is also dependend on this parameter. Can some one help me? Thanks!!

[@shivathapa3934]

Hi sir. What if there was a second order derivative in the first equation (d^2A/dt^2)? What modification would be needed in the code?

[@analuciavillelaloaiza8930]

Kinda late, but you would need to do a substitution to make two first order edos instead of one of second order. Like: u=dA/dt & du/dt=d2A/dt2.

[@lunthangify]

p = odeint(odes,r0,θ0,ϕ0,x0,z0,t)
NameError: name 'r0' is not defined