Vežbe 05 (18.11.2011)

Metod Runge-Kutta

LHO

  • Primer sa predavanja: Fortran / C++ (isto rešenje u oba programska jezika)
  • C++ drugi način:
    • Jednostavan algoritam (direktna primena jednačina): rk-0.cpp
    • Modifikovan kod uz primenu funkcije rk4_old (ova f-ja radi samo u ovom i sličnim jednostavnim slučajevima): rk-1.cpp
    • Konačna verzija rk4 funkcije za simultano rešavanje jednačina: rk-2.cpp
  • Mathematica 8 (notebook)

LHO na koji deluje viskozna sila

  • Rešavanje jednačine oblika frac {d^2 x} {d t^2} + 2 beta frac {dx} {dt} + omega_0^2 = 0
  • Rešenje viskozna-sila.cpp (korišćenjem rk4 funkcije)
  • Pogrešno rešenje viskozna-sila-1.cpp (korišćenjem rk4_old funkcije). Primer pokazuje razliku između ove dve funkcije i zbog čega rk4_old algoritam ne može da se koristi.
  • Mathematica 8 (notebook)

Funkcija koja realizuje Runge-Kuta metod

Imate pitanje ili komentar?