File html/i5a.html, c=[
0.166667 , 0.333333, 0.5, 0.666667, 0.833333, 1]
AUBV (scaled)
| 0.333333 | -0 | 0 | 0 | -0 | 0 |
| -0.0864424 | 0.333333 | -0 | -0 | 0 | -0 |
| -2.1685 | -2.23619 | 0.333333 | -0 | 0 | -0 |
| -4.73161 | -3.89265 | -0.276319 | 0.333333 | 0 | -0 |
| -6.75188 | -7.90757 | 0.790245 | -0.448352 | 0.333333 | -0 |
| -4.26488 | -11.932 | 3.38925 | -2.2397 | 0.662808 | 0.333333 |
|
| 1 | -0.166667 | -0.0833333 | -0.0231481 | -0.00540123 | -0.00115741 |
| 1 | 0.0864424 | -0.082297 | -0.0668705 | -0.0354363 | -0.0161274 |
| 1 | 4.57136 | 2.13029 | 0.801105 | 0.267278 | 0.0734859 |
| 1 | 9.23392 | 4.44862 | 1.75094 | 0.60494 | 0.147361 |
| 1 | 14.8175 | 7.46878 | 3.08789 | 1.14343 | 0.308158 |
| 1 | 15.0512 | 8.20199 | 3.39623 | 0.938913 | -0.358964 |
|
| -4.26488 | -11.932 | 3.38925 | -2.2397 | 0.662808 | 0.333333 |
| 0 | 0 | -0 | 0 | 0 | 1 |
| -14.389 | -5.67602 | 13.1001 | 12.8472 | -15.3351 | 3.34034 |
| 9.13285 | 11.3394 | -10.8492 | -13.8107 | 13.6236 | -1.96365 |
| -9.72217 | 2.55393 | -2.04938 | 7.6132 | -5.40784 | 1.29153 |
| -1.43374 | 0.0716829 | 0.00665129 | 1.25309 | -0.987629 | 0.209082 |
|
| 1 | 15.0512 | 8.20199 | 3.39623 | 0.938913 | -0.358964 |
| 0 | 0 | 0 | 0 | 0 | 0 |
| 0 | 6.11251 | -1.77151 | 1.06355 | 7.46759 | 13.8988 |
| 0 | -7.47232 | -0.119121 | 0.51875 | 0.262291 | 3.26351 |
| 0 | 5.72072 | -0.133454 | -1.2633 | 0.155329 | 4.58315 |
| 0 | 0.88086 | -0.019435 | -0.149884 | -0.0226945 | 1.09743 |
|
Maximum coefficient scaled is 15.335113 .
irks(0.333333,[0.166667,0.333333,0.5,0.666667,0.833333,1], [-0.0489,0.4228,-0.8814,0.9021],[ 0.42033 1 0 0; -0.00527116 0.489107 0.337476 1; 1 0 0 0; 0.53546 0.743626 1 0; ])
[alpha,beta,gamma,stageerror]
| -1.71468e-05 | 3.77113 | -3.77082 | 0.000595659 |
| 0 | 0 | 0 | -0.00152926 |
| -0.0489 | 0.968295 | -0.924182 | -0.00488393 |
| 0.4228 | -0.640452 | 0.210144 | -0.00652775 |
| -0.8814 | 4.23804 | -3.02013 | -0.00478093 |
| 0.9021 | 0.927348 | -2.62253 | -1.71468e-05 |
Error estimation
phi=[-7776.000000 38880.000000 -77760.000000 77760.000000 -38880.000000 7776.000000 ]
phi1=[-27058.262446 87085.656114 -142034.208153 125965.656114 -58162.262446 10989.710408 ] phi0=3213.710408 gives y^{p+1}(t_n-theta*h) at theta=0.347786