Coding a Fourth-Order Runge-Kutta Integrator in Python and Matlab

Поділитися
Вставка
  • Опубліковано 21 лип 2024
  • In this video, I code up a 4th-order accurate Runge-Kutta integrator in Python and Matlab, and then I use this integrator to simulate the chaotic Lorenz 1963 system.
    Playlist: • Engineering Math: Diff...
    Course Website: faculty.washington.edu/sbrunto...
    @eigensteve on Twitter
    eigensteve.com
    databookuw.com
    This video was produced at the University of Washington
    %%% CHAPTERS %%%
    0:00 Problem setup and Lorenz 1963 example
    6:16 Matlab code example
    24:48 Python code example
  • Наука та технологія

КОМЕНТАРІ • 14

  • @GoobNoob
    @GoobNoob Рік тому +4

    Wow. What a fantastic video and teacher. You really nailed it. This should be the most recommended R-K Python/Matlab tutorial and explanation on UA-cam. The others don't match your level.

  • @nesslange1833
    @nesslange1833 Рік тому +1

    Your Lecture always has been great and full of detail, but coding code from scratch makes it all even more approachable!

  • @stavrosoik129
    @stavrosoik129 Рік тому

    This video came right on time where I need to translate a pre existing epidemiological model from matlab to R where we are using ode45 to solve the ODEs! Thanks!

  • @MrHaggyy
    @MrHaggyy Рік тому +1

    Great lecture and code to play around and work with. You mentioned stiff ode for a view videos, and google got me an okisch answer. But a video with some intuition what things might brake ode45 or potential potholes might come in handy.

  • @hoseinzahedifar1562
    @hoseinzahedifar1562 Рік тому +1

    Thank you very much...🌺🌺🌺
    comet3 is a really cool function...🤩🤩

  • @deralaw
    @deralaw 9 місяців тому

    An excellent lecture! Please do a session on Implicit RK4 including the matlab coding. Thank you

  • @lucynowacki3327
    @lucynowacki3327 Рік тому +2

    Can we have implicit method based on a fixed-point theorem?

  • @tabishbadar1320
    @tabishbadar1320 Рік тому

    Firstly, thanks. Secondly, I computed the trajectories using three methods, 1. Forward Euler, 2. RK2, and 3. RK4 of the chaotic system. Third, using the reference output from MATLAB's ode45, the RK2 had the most minor error value (used MATLAB's immse function to compute errors). Fourth, why so?

  • @user-hu3gn1gi5f
    @user-hu3gn1gi5f Рік тому

    Could you tell me somebody what pen is used in this video? Awesome!!

  • @whdaffer1
    @whdaffer1 5 місяців тому

    A quick question: how is climate a chaotic system? I understand the reference to Weather but I'm a little unclear on why climate would be chaotic.
    They are, after all, two different sorts of problems. Weather is an initial value problem, climate is a boundary value problem and one of the boundary values is pretty well understood: namely the insolation.

  • @TheTimeRay
    @TheTimeRay Рік тому +3

    brightness\contrast 😕

  • @charlesspringer4709
    @charlesspringer4709 Рік тому

    Is there any info on how these videos are made?
    Also, the derivation in Numerical Recipes is wrong but the code is correct. I assumed when I noticed this in 1990 that it is an error copied from book to book.