File html/e2b.html, c=[
0.5 , 1, 1]
AUBV (scaled)
| 0 | 0 | 0 |
| 0.446158 | 0 | 0 |
| -0.139694 | 0.460582 | 0 |
|
| 1 | 0.5 | 0.25 |
| 1 | 0.553842 | 0.553842 |
| 1 | 0.679112 | 0.21853 |
|
| -0.139694 | 0.460582 | 0 |
| 0 | 0 | 1 |
| 0.561552 | -0.561553 | 0.780777 |
|
| 1 | 0.679112 | 0.21853 |
| 0 | 0 | 0 |
| 0 | -0.780776 | 0 |
|
Maximum coefficient scaled is 0.780777 .
irks(0,[0.5,1,1], [0.640388],[ 1; ])
[alpha,beta,gamma,stageerror]
| 0.0238095 | -0.0580673 | 0.099734 | 0.100882 |
| 0 | 0 | -0 | 0.288234 |
| 0.640388 | -0.236861 | -0.173236 | 0.0238095 |
Error estimation
phi=[0.000000 0.000000 0.000000 ]
phi1=[-8.000000 2.691950 1.308050 ] phi0=4.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000