Mercurial > repos > public > sbplib

comparison +time/+rk/butcherTableau.m @ 929:0344fff87139 feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve error msg in rk.butcherTableau
author Martin Almquist <malmquist@stanford.edu>
date Mon, 03 Dec 2018 15:40:26 -0800
parents d1c1615bd1a5
children
comparison
equal deleted inserted replaced
918:679f4ddd982f 929:0344fff87139
45 a(6,1) = 71/2400 - 587/12000*sqrt(5); a(6,2) = 187/480 - 391/2400*sqrt(5); a(6,3) = -38/75 + 26/375*sqrt(5); a(6,4) = 27/80 - 3/400*sqrt(5); a(6,5) = (1+sqrt(5))/4; 45 a(6,1) = 71/2400 - 587/12000*sqrt(5); a(6,2) = 187/480 - 391/2400*sqrt(5); a(6,3) = -38/75 + 26/375*sqrt(5); a(6,4) = 27/80 - 3/400*sqrt(5); a(6,5) = (1+sqrt(5))/4;
46 a(7,1) = -49/480 + 43/160*sqrt(5); a(7,2) = -425/96 + 51/32*sqrt(5); a(7,3) = 52/15 - 4/5*sqrt(5); a(7,4) = -27/16 + 3/16*sqrt(5); a(7,5) = 5/4 - 3/4*sqrt(5); a(7,6) = 5/2 - 1/2*sqrt(5); 46 a(7,1) = -49/480 + 43/160*sqrt(5); a(7,2) = -425/96 + 51/32*sqrt(5); a(7,3) = 52/15 - 4/5*sqrt(5); a(7,4) = -27/16 + 3/16*sqrt(5); a(7,5) = 5/4 - 3/4*sqrt(5); a(7,6) = 5/2 - 1/2*sqrt(5);
47 b = [1/12 0 0 0 5/12 5/12 1/12]; 47 b = [1/12 0 0 0 5/12 5/12 1/12];
48 c = [0, 4/7, 5/7, 6/7, (5-sqrt(5))/10, (5+sqrt(5))/10, 1]; 48 c = [0, 4/7, 5/7, 6/7, (5-sqrt(5))/10, (5+sqrt(5))/10, 1];
49 otherwise 49 otherwise
50 error('That Runge-Kutta method is not implemented', method) 50 error('Runge-Kutta method %s is not implemented', method)
51 end 51 end