File html/e5a.html, c=[
0.166667 , 0.333333, 0.5, 0.666667, 0.833333, 1]
AUBV (scaled)
| -0 | 0 | -0 | 0 | -0 | 0 |
| 0.33754 | 0 | -0 | 0 | -0 | 0 |
| 0.536092 | 0.242934 | 0 | -0 | 0 | -0 |
| 1.13918 | -0.258709 | 0.394478 | -0 | 0 | -0 |
| 0.83731 | 0.364167 | -0.407896 | 0.513482 | 0 | -0 |
| 1.66114 | -0.824245 | 0.192754 | 0.356627 | 0.2483 | -0 |
|
| 1 | 0.166667 | 0.0277778 | 0.00462963 | 0.000771605 | 0.000128601 |
| 1 | -0.004207 | -0.00140233 | 0.00890868 | 0.00609493 | 0.00281299 |
| 1 | -0.279026 | -0.0906532 | -0.000652275 | 0.0165822 | 0.0141858 |
| 1 | -0.608284 | -0.157288 | -0.00825744 | 0.0175233 | 0.0199877 |
| 1 | -0.473731 | -0.104185 | 0.00881732 | 0.00817274 | -0.00350797 |
| 1 | -0.634577 | -0.0863068 | -0.00104005 | -0.0024664 | 0.0332931 |
|
| 1.66114 | -0.824245 | 0.192754 | 0.356627 | 0.2483 | 0 |
| -0 | -0 | -0 | 0 | 0 | 1 |
| -4.00226 | 0.0686735 | 5.93957 | -3.74506 | -4.66044 | 4.55405 |
| -12.1518 | 9.47631 | 4.73845 | -1.3149 | -12.2351 | 7.56982 |
| -13.3586 | 14.8117 | -3.10503 | 3.03249 | -9.72373 | 4.91779 |
| -5.86424 | 6.84466 | -2.93117 | 1.73188 | -2.54679 | 1.12613 |
|
| 1 | -0.634577 | -0.0863068 | -0.00104005 | -0.0024664 | 0.0332931 |
| 0 | 0 | 0 | 0 | 0 | 0 |
| 0 | 1.84546 | 0.00143678 | -0.103549 | 0.104581 | 0.321193 |
| 0 | 3.91728 | 6.89854e-05 | -0.16635 | 0.0531044 | 0.932886 |
| 0 | 3.42539 | 0.0107406 | -0.0341797 | -0.15104 | 0.970023 |
| 0 | 1.63952 | 0.00602414 | 0.0237623 | -0.101606 | 0.315953 |
|
Maximum coefficient scaled is 14.811692 .
irks(-1.278e-14,[0.166667,0.333333,0.5,0.666667,0.833333,1], [0.0074,0.0674,0.3361,0.8734],[ -0.22483 1 0 0; 0.342782 0.115651 0.591259 1; 0.7828 0.863139 1 0; 1 0 0 0; ])
[alpha,beta,gamma,stageerror]
| 1.90259e-05 | 0.000346447 | -0.000148034 | 0.000166555 |
| 0 | -0 | 0 | 0.000202257 |
| 0.0074 | -0.00937364 | 0.00291048 | 5.52048e-06 |
| 0.0674 | -0.0693327 | 0.0178656 | -0.000256992 |
| 0.3361 | -0.277229 | 0.0510796 | -0.000203034 |
| 0.8734 | -0.487808 | 0.0610808 | 1.90259e-05 |
Error estimation
phi=[-7776.000000 38880.000000 -77760.000000 77760.000000 -38880.000000 7776.000000 ]
phi1=[42423.069333 -86617.673334 89570.231111 -47737.673334 11319.069333 -590.511556 ] phi0=-8366.511556 gives y^{p+1}(t_n-theta*h) at theta=0.595990