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

AUBV (scaled)

0 0 0 0 0
0.198404 0 0 0 0
-1.71421 1.17656 0 0 -0
0.254936 0.12468 0.371168 0 -0
-1.63318 0.988223 0.187848 0.321139 -0
1 0.2 0.04 0.008 0.0016
1 0.201596 0.0806386 0.0401916 0.0192511
1 1.13766 0.104439 -0.143042 -0.116744
1 0.0492161 -0.00712016 0.0206996 0.0488346
1 1.13598 0.123456 -0.0978265 -0.0207155
-1.63318 0.988223 0.187848 0.321139 0
0 0 0 0 1
-0.749657 -0.120571 -1.31185 -0.390745 1.81176
0.347197 0.0452513 -1.17491 -1.44482 1.78275
0.817348 -0.194357 -0.191474 -0.570419 0.490196
1 1.13598 0.123456 -0.0978265 -0.0207155
0 -0 -0 -0 -0
0 0.761055 -0.0278022 -0.120437 0.741477
0 0.444535 -0.018969 -0.368666 0.820434
0 -0.351294 -0.00940577 -0.173381 0.396468

Maximum coefficient scaled is 1.811765 .

irks(0,[0.2,0.4,0.6,0.8,1], [0.085,0.308,0.9092],[ 1 0 0; 0.150324 0.539995 1; 0.547644 1 0; ])

[alpha,beta,gamma,stageerror]

0.000333333 -0.00119969 0.00258858 0.00217395
-0 0 0 0.00629171
0.085 -0.0895058 0.0122238 -0.00781957
0.308 -0.219578 0.0245447 0.00318636
0.9092 -0.503946 -0.0146989 0.000333333

Error estimation

phi=[625.000000 -2500.000000 3750.000000 -2500.000000 625.000000 ]
phi1=[-1138.952734 1027.905467 222.094533 -736.047266 272.209453 ] phi0=352.790547 gives y^{p+1}(t_n-theta*h) at theta=0.512893