File html/i5b.html, c=[ 0.2 , 0.4, 0.6, 0.8, 1, 1]

AUBV (scaled)

0.333333 0 -0 -0 0 0
-0.0499118 0.333333 -0 -0 0 0
0.565215 -0.349134 0.333333 -0 0 0
0.535747 -1.71315 0.883819 0.333333 0 0
0.973119 -3.2072 1.57123 0.104471 0.333333 0
6.56081 -9.45014 1.08954 6.83006 -3.46669 0.333333
1 -0.133333 -0.0933333 -0.032 -0.00906667 -0.00234667
1 0.116578 -0.0867019 -0.0900106 -0.0581362 -0.0320274
1 0.0505857 0.0132212 -0.0442414 -0.0871085 -0.0980726
1 0.760254 0.202308 -0.3245 -0.615263 -0.712703
1 1.22504 0.457212 -0.474835 -1.11493 -1.49605
1 -0.896918 -0.0330515 -0.141593 0.813371 3.12991
6.56081 -9.45014 1.08954 6.83006 -3.46669 0.333333
0 0 -0 0 0 1
3.15637 -7.8322 2.55608 -3.12908 3.53191 -0.174677
-2.60131 7.66182 -1.04218 -8.97901 5.01739 0.283553
0.132881 0.0450398 1.75874 -5.78241 3.30753 -0.157463
1.01517 -2.26648 0.998732 -0.541188 0.516769 -0.161805
1 -0.896918 -0.0330515 -0.141593 0.813371 3.12991
0 0 0 0 0 0
0 1.8916 1.22798 -0.443736 -1.32498 -1.05687
0 -0.340269 -0.0737892 0.0969065 0.207498 1.59972
0 0.695683 0.752049 -0.284981 -1.29323 -0.0544567
0 0.43881 0.364615 -0.13835 -0.626674 -0.0316581

Maximum coefficient scaled is 9.450143 .

irks(0.333333,[0.2,0.4,0.6,0.8,1,1], [-0.0575,0.0472,-0.1525,-0.1387],[ 0.0819494 0.057375 1 0; 1 0 0 0; -0.302922 0.412735 -0.0146192 1; -0.78609 1 0 0; ])

[alpha,beta,gamma,stageerror]

-1.71468e-05 0.0435988 -0.0432975 0.00249112
-0 0 0 0.0021684
-0.0575 0.100943 -0.108918 -0.0001842
0.0472 0.0669665 -0.11792 -0.00460993
-0.1525 0.472942 -0.446894 -0.0104878
-0.1387 1.10315 -1.27489 -1.71468e-05

Error estimation

phi=[0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ]
phi1=[15625.000000 -31250.000000 31250.000000 -15625.000000 3570.131940 -445.131940 ] phi0=-3125.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000