File html/e2a.html, c=[
0.333333 , 0.666667, 1]
AUBV (scaled)
| 0 | 0 | 0 |
| 0.428571 | 0 | 0 |
| -0.107143 | 0.75 | 0 |
|
| 1 | 0.333333 | 0.111111 |
| 1 | 0.238095 | 0.15873 |
| 1 | 0.357143 | 0.0714286 |
|
| -0.107143 | 0.75 | 0 |
| 0 | 0 | 1 |
| 1.5 | -1.5 | 1 |
|
| 1 | 0.357143 | 0.0714286 |
| 0 | 0 | 0 |
| 0 | -1 | -0 |
|
Maximum coefficient scaled is 1.500000 .
irks(0,[0.333333,0.666667,1], [0.5],[ 1; ])
[alpha,beta,gamma,stageerror]
| 0.0238095 | -0.00198413 | 0.0436508 | 0.0339506 |
| 0 | 0 | 0 | 0.0652557 |
| 0.5 | -0.203704 | -0.0462963 | 0.0238095 |
Error estimation
phi=[9.000000 -18.000000 9.000000 ]
phi1=[-19.018868 10.018868 -0.339623 ] phi0=9.339623 gives y^{p+1}(t_n-theta*h) at theta=0.679245