Mercurial > repos > public > sbplib

view +time/+rkparameters/rk4.m @ 1275:8ff3a95ad7cc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change metric computation in CurvilinearUpwind so that one can use odd orders for the diffOp
author Martin Almquist <malmquist@stanford.edu>
date Wed, 03 Jun 2020 10:29:53 -0700
parents ee4cfb37534d
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