Mercurial > repos > public > sbplib

view +time/+rkparameters/rk4.m @ 737:521076a6ffa5 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update Runge Kutta for discrete data to return stage rates if needed.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 25 Apr 2018 16:58:45 -0700
parents 9c3c180c7ed3
children
line wrap: on
line source

function [a,b,c,s] = rk4()

% Butcher tableau for classical RK$
s = 4;
a = sparse(s,s);
a(2,1) = 1/2;
a(3,2) = 1/2;
a(4,3) = 1;
b = 1/6*[1; 2; 2; 1];
c = [0; 1/2; 1/2; 1];

end