Mercurial > repos > public > sbplib

view +time/+rkparameters/rk4.m @ 885:18e10217dca9 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Generalize Staggered1DAcoustics to handle variable coefficients defined by both function handles and vectors.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 04 Nov 2018 12:36:30 -0800
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