File html/i4a.html, c=[
0.2 , 0.4, 0.6, 0.8, 1]
AUBV (scaled)
| 0.272727 | 0 | 0 | -0 | 0 |
| -0.10398 | 0.272727 | 0 | -0 | 0 |
| -1.58615 | 0.744169 | 0.272727 | -0 | 0 |
| -0.873658 | 0.537885 | -0.163299 | 0.272727 | 0 |
| 0.295489 | -1.18482 | -0.66803 | 1.00717 | 0.272727 |
|
| 1 | -0.0727273 | -0.0690909 | -0.0247273 | -0.00712727 |
| 1 | 0.231253 | -0.0165898 | -0.0544315 | -0.0408908 |
| 1 | 1.16926 | 0.0718537 | -0.245408 | -0.245787 |
| 1 | 1.02634 | 0.31875 | 0.0113804 | -0.117597 |
| 1 | 1.27746 | 0.474372 | -0.497217 | -1.28255 |
|
| 0.295489 | -1.18482 | -0.66803 | 1.00717 | 0.272727 |
| 0 | 0 | 0 | 0 | 1 |
| -2.02941 | -2.00462 | -0.694653 | 2.22612 | 0.316166 |
| 1.39282 | -0.377735 | 1.14441 | -0.870992 | 0.751488 |
| 0.531901 | 0.117034 | 0.371322 | -0.449735 | 0.201056 |
|
| 1 | 1.27746 | 0.474372 | -0.497217 | -1.28255 |
| 0 | -0 | -0 | -0 | -0 |
| 0 | 2.1864 | 0.0549222 | -0.266674 | 1.35455 |
| 0 | -2.03999 | -1.73762 | -0.803952 | 1.84119 |
| 0 | -0.771578 | -0.43451 | -0.260709 | 0.74903 |
|
Maximum coefficient scaled is 2.226120 .
irks(0.272727,[0.2,0.4,0.6,0.8,1], [0.1715,-0.1238,0.6617],[ -0.0811981 0.407917 1; 1 0 0; 0.827005 1 0; ])
[alpha,beta,gamma,stageerror]
| 0.000207802 | -0.127948 | 0.127762 | -0.00562636 |
| -0 | 0 | 0 | -0.00162551 |
| 0.1715 | -0.20204 | 0.0885587 | 0.00293573 |
| -0.1238 | 0.7859 | -0.532481 | 0.0222982 |
| 0.6617 | 0.220794 | -0.782441 | 0.000207802 |
Error estimation
phi=[625.000000 -2500.000000 3750.000000 -2500.000000 625.000000 ]
phi1=[-1723.513482 2197.026964 -947.026964 -151.486518 155.297304 ] phi0=469.702696 gives y^{p+1}(t_n-theta*h) at theta=0.550305