Mercurial > repos > public > sbplib

diff +time/+rk4/rungekutta_4.m @ 0:48b6fb693025

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Initial commit.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 17 Sep 2015 10:12:50 +0200
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/+time/+rk4/rungekutta_4.m	Thu Sep 17 10:12:50 2015 +0200
@@ -0,0 +1,10 @@
+% Takes one time step of size k using the rungekutta method
+% starting from v_0 and where the function F(v,t) gives the
+% time derivatives.
+function v = rungekutta_4(v, t , k, F)
+    k1 = F(v         ,t      );
+    k2 = F(v+0.5*k*k1,t+0.5*k);
+    k3 = F(v+0.5*k*k2,t+0.5*k);
+    k4 = F(v+    k*k3,t+    k);
+    v = v + (1/6)*(k1+2*(k2+k3)+k4)*k;
+end
\ No newline at end of file