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

AUBV (scaled)

0.272727 0 0 -0 0
-1.37006 0.272727 0 0 0
1.15639 0.363734 0.272727 -0 0
-1.40391 1.40166 0.253643 0.272727 0
1.35248 -0.0571837 0.276579 -0.171017 0.272727
1 -0.0227273 -0.0738636 -0.0355114 -0.0131392
1 1.59734 0.662305 0.177342 0.0117654
1 -1.04285 -0.78852 -0.527976 -0.397963
1 0.475885 -0.62562 -1.03421 -1.13202
1 -0.673583 -0.237343 0.017442 0.0704959
1.35248 -0.0571837 0.276579 -0.171017 0.272727
0 0 0 0 1
-1.4554 -0.775621 -1.53558 0.696452 1.54112
-0.619158 -1.18882 -1.97884 0.829075 1.10769
0.0121999 -0.420177 -0.655617 0.262169 0.283811
1 -0.673583 -0.237343 0.017442 0.0704959
0 0 0 0 0
0 1.52902 0.33154 -0.266829 0.119768
0 1.85005 0.593129 -0.463297 0.225337
0 0.517613 0.305541 -0.218743 0.131757

Maximum coefficient scaled is 1.978838 .

irks(0.272727,[0.25,0.5,0.75,1,1], [0.0234,0.1864,0.703],[ -0.0829033 0.254099 1; -0.196813 1 0; 1 0 0; ])

[alpha,beta,gamma,stageerror]

0.000207802 0.000268293 -0.000454921 -0.00238854
0 -0 -0 0.0133762
0.0234 -0.0113625 -0.00508768 -0.0400385
0.1864 -0.0869652 -0.020674 -0.0824038
0.703 -0.281056 -0.0267534 0.000207802

Error estimation

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