File html/i2a.html, c=[
0.333333 , 0.666667, 1]
AUBV (scaled)
| 0.444444 | -0 | -0 |
| 1.03751 | 0.444444 | -0 |
| 0.767025 | -0.38114 | 0.444444 |
|
| 1 | -0.111111 | -0.185185 |
| 1 | -0.815284 | -0.839819 |
| 1 | 0.169671 | 0.107948 |
|
| 0.767025 | -0.38114 | 0.444444 |
| 0 | 0 | 1 |
| -1.03752 | 0.310864 | 0.638599 |
|
| 1 | 0.169671 | 0.107948 |
| 0 | 0 | 0 |
| 0 | 0.0880614 | 0 |
|
Maximum coefficient scaled is 1.037524 .
irks(0.444444,[0.333333,0.666667,1], [0.33852],[ 1; ])
[alpha,beta,gamma,stageerror]
| 0.0048011 | -0.0220597 | 0.0564104 | -0.049863 |
| 0 | 0 | 0 | -0.24917 |
| 0.33852 | -0.0692799 | -0.0453159 | 0.0048011 |
Error estimation
phi=[9.000000 -18.000000 9.000000 ]
phi1=[-11.305393 2.305393 2.231536 ] phi0=6.768464 gives y^{p+1}(t_n-theta*h) at theta=0.584017