File html/i4b.html, c=[
0.25 , 0.5, 0.75, 1, 1]
AUBV (scaled)
| 0.272727 | 0 | 0 | -0 | 0 |
| -1.37006 | 0.272727 | 0 | 0 | 0 |
| 1.15639 | 0.363734 | 0.272727 | -0 | 0 |
| -1.40391 | 1.40166 | 0.253643 | 0.272727 | 0 |
| 1.35248 | -0.0571837 | 0.276579 | -0.171017 | 0.272727 |
|
| 1 | -0.0227273 | -0.0738636 | -0.0355114 | -0.0131392 |
| 1 | 1.59734 | 0.662305 | 0.177342 | 0.0117654 |
| 1 | -1.04285 | -0.78852 | -0.527976 | -0.397963 |
| 1 | 0.475885 | -0.62562 | -1.03421 | -1.13202 |
| 1 | -0.673583 | -0.237343 | 0.017442 | 0.0704959 |
|
| 1.35248 | -0.0571837 | 0.276579 | -0.171017 | 0.272727 |
| 0 | 0 | 0 | 0 | 1 |
| -1.4554 | -0.775621 | -1.53558 | 0.696452 | 1.54112 |
| -0.619158 | -1.18882 | -1.97884 | 0.829075 | 1.10769 |
| 0.0121999 | -0.420177 | -0.655617 | 0.262169 | 0.283811 |
|
| 1 | -0.673583 | -0.237343 | 0.017442 | 0.0704959 |
| 0 | 0 | 0 | 0 | 0 |
| 0 | 1.52902 | 0.33154 | -0.266829 | 0.119768 |
| 0 | 1.85005 | 0.593129 | -0.463297 | 0.225337 |
| 0 | 0.517613 | 0.305541 | -0.218743 | 0.131757 |
|
Maximum coefficient scaled is 1.978838 .
irks(0.272727,[0.25,0.5,0.75,1,1], [0.0234,0.1864,0.703],[ -0.0829033 0.254099 1; -0.196813 1 0; 1 0 0; ])
[alpha,beta,gamma,stageerror]
| 0.000207802 | 0.000268293 | -0.000454921 | -0.00238854 |
| 0 | -0 | -0 | 0.0133762 |
| 0.0234 | -0.0113625 | -0.00508768 | -0.0400385 |
| 0.1864 | -0.0869652 | -0.020674 | -0.0824038 |
| 0.703 | -0.281056 | -0.0267534 | 0.000207802 |
Error estimation
phi=[0.000000 0.000000 0.000000 0.000000 0.000000 ]
phi1=[-1024.000000 1536.000000 -1024.000000 775.246414 -519.246414 ] phi0=256.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000