File html/i3a.html, c=[
0.25 , 0.5, 0.75, 1]
AUBV (scaled)
| 0.225 | -0 | 0 | -0 |
| 0.211287 | 0.225 | 0 | -0 |
| 0.946338 | -0.342943 | 0.225 | -0 |
| 0.52149 | -0.662474 | 0.490476 | 0.225 |
|
| 1 | 0.025 | -0.05 | -0.0265625 |
| 1 | 0.063713 | -0.0806435 | -0.0833663 |
| 1 | -0.0783954 | 0.0947737 | 0.121956 |
| 1 | 0.425507 | 0.216014 | -0.103603 |
|
| 0.52149 | -0.662474 | 0.490476 | 0.225 |
| 0 | 0 | 0 | 1 |
| -0.0423385 | 0.695379 | -0.784079 | 1.0116 |
| 0.077564 | 0.246379 | -0.321806 | 0.274145 |
|
| 1 | 0.425507 | 0.216014 | -0.103603 |
| 0 | 0 | 0 | 0 |
| 0 | -0.880558 | -0.521284 | 0.774748 |
| 0 | -0.276282 | -0.350743 | 0.521284 |
|
Maximum coefficient scaled is 1.011596 .
irks(0.225,[0.25,0.5,0.75,1], [0.3312,1.005],[ 1 0; -0.49541 1; ])
[alpha,beta,gamma,stageerror]
| 0.000542057 | -0.0229588 | 0.0246866 | -0.0131524 |
| 0 | -0 | 0 | -0.029952 |
| 0.3312 | -0.253506 | -0.0793501 | 0.0381657 |
| 1.005 | -0.595887 | -0.0829382 | 0.000542057 |
Error estimation
phi=[-64.000000 192.000000 -192.000000 64.000000 ]
phi1=[110.801474 -70.202212 -17.198526 20.299631 ] phi0=-43.700369 gives y^{p+1}(t_n-theta*h) at theta=0.545705