File html/e4b.html, c=[ 0.25 , 0.5, 0.75, 1, 1]

AUBV (scaled)

0 0 0 0 0
0.320112 0 -0 0 0
-0.475669 0.537635 0 -0 -0
-0.456266 0.421881 0.247582 -0 0
-0.463912 -0.125869 -0.410284 0.701687 -0
1 0.25 0.0625 0.015625 0.00390625
1 0.179888 0.0899439 0.064979 0.042493
1 0.688034 0.2627 0.107837 0.0773183
1 0.786803 0.434879 0.351344 0.399781
1 1.29838 0.569877 -0.231322 -1.02247
-0.463912 -0.125869 -0.410284 0.701687 0
0 0 0 0 1
-1.39531 -1.09018 -0.513546 1.53777 0.0820588
0.934241 -1.06952 -0.940365 0.603287 0.735
0.413923 -0.177612 -0.170126 -0.0282253 0.245098
1 1.29838 0.569877 -0.231322 -1.02247
0 -0 -0 -0 -0
0 1.37921 0.318494 0.0863759 1.01959
0 -0.262646 -0.66363 -0.801027 0.710088
0 -0.283058 -0.207906 -0.307933 0.482532

Maximum coefficient scaled is 1.537770 .

irks(0,[0.25,0.5,0.75,1,1], [0.17,0.0279,0.7497],[ -0.324197 0.435914 1; 1 0 0; 0.876408 1 0; ])

[alpha,beta,gamma,stageerror]

0.000333333 -0.0741312 0.0755201 0.00551532
-0 0 0 0.00948308
0.17 -0.228215 0.100766 0.025901
0.0279 0.286427 -0.137343 0.0551316
0.7497 -0.130257 -0.403893 0.000333333

Error estimation

phi=[0.000000 0.000000 0.000000 0.000000 0.000000 ]
phi1=[-1024.000000 1536.000000 -1024.000000 319.700178 -63.700178 ] phi0=256.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000