Mercurial > repos > public > sbplib

comparison +time/+rk/rungekutta_4.m @ 887:50d5a3843099 feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Rename package rk4 to rk
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 15 Nov 2018 16:42:58 -0800
parents +time/+rk4/rungekutta_4.m@48b6fb693025
children
comparison
equal deleted inserted replaced
886:8894e9c49e40 887:50d5a3843099
1 % Takes one time step of size k using the rungekutta method
2 % starting from v_0 and where the function F(v,t) gives the
3 % time derivatives.
4 function v = rungekutta_4(v, t , k, F)
5 k1 = F(v ,t );
6 k2 = F(v+0.5*k*k1,t+0.5*k);
7 k3 = F(v+0.5*k*k2,t+0.5*k);
8 k4 = F(v+ k*k3,t+ k);
9 v = v + (1/6)*(k1+2*(k2+k3)+k4)*k;
10 end