File html/i2b.html, c=[ 0.5 , 1, 1]

AUBV (scaled)

0.444444 -0 0
0.885561 0.444444 0
0.604175 -0.262731 0.444444
1 0.0555556 -0.194444
1 -0.330006 -0.77445
1 0.214111 0.0323975
0.604175 -0.262731 0.444444
0 0 1
-0.88556 0.251269 0.691511
1 0.214111 0.0323975
0 0 0
0 -0.05722 0

Maximum coefficient scaled is 0.885561 .

irks(0.444444,[0.5,1,1], [0.27861],[ 1; ])

[alpha,beta,gamma,stageerror]

0.0048011 -0.00211258 0.0364633 -0.0618093
0 0 -0 -0.274136
0.27861 -0.0559717 -0.0216519 0.0048011

Error estimation

phi=[0.000000 0.000000 0.000000 ]
phi1=[-8.000000 1.841561 2.158439 ] phi0=4.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000