
global RK_s RK_a RK_b RK_c RK_sh RK_bh RK_CFL RK_em_CFL;


rho = 1.40;
rho2 = 1.40;

RK_CFL = rho;
RK_em_CFL = 1;

RK_s = 6;
RK_sh = 5;


RK_a = [0                0              0              0        0      0; ...
	1                0              0              0        0      0; ...
	1/4              1/4            0              0        0      0; ...
	2046/15625      -454/15625      1533/15625     0        0      0; ...
       -739/5625         511/5625      -566/16875      20/27    0      0; ...
	11822/21875     -6928/21875    -4269/21875    -4/7      54/35  0];

RK_c = [0  1  1/2  1/5  2/3  1];

RK_bh = [1/6     1/6    2/3    0      0      0 ];

RK_b = [1/24    0      0      125/336    27/56    5/48];
