File html/i3b.html, c=[
0.333333 , 0.666667, 1, 1]
AUBV (scaled)
| 0.225 | -0 | 0 | 0 |
| 0.37915 | 0.225 | 0 | 0 |
| 0.576331 | 0.157279 | 0.225 | 0 |
| 0.740628 | 0.864647 | -0.404278 | 0.225 |
|
| 1 | 0.108333 | -0.0388889 | -0.037963 |
| 1 | 0.0625165 | -0.108322 | -0.130087 |
| 1 | 0.0413899 | -0.0439261 | -0.0768157 |
| 1 | -0.425997 | -0.288058 | 0.138096 |
|
| 0.740628 | 0.864647 | -0.404278 | 0.225 |
| 0 | 0 | 0 | 1 |
| -0.715262 | -0.740829 | 0.350709 | 0.892767 |
| -0.570461 | -0.893763 | 0.381671 | 0.404828 |
|
| 1 | -0.425997 | -0.288058 | 0.138096 |
| 0 | 0 | 0 | 0 |
| 0 | 0.212614 | -0.0223407 | 0.495762 |
| 0 | 0.677725 | -0.00100675 | 0.0223407 |
|
Maximum coefficient scaled is 0.893763 .
irks(0.225,[0.333333,0.666667,1,1], [0.2463,0.5101],[ 1 0; 0.13519 1; ])
[alpha,beta,gamma,stageerror]
| 0.000542057 | 0.0323451 | -0.0306173 | -0.00889114 |
| -0 | -0 | -0 | -0.0296206 |
| 0.2463 | -0.193098 | 0.0674604 | -0.0190979 |
| 0.5101 | -0.163906 | 0.150004 | 0.000542057 |
Error estimation
phi=[0.000000 0.000000 0.000000 0.000000 ]
phi1=[81.000000 -81.000000 86.238349 -59.238349 ] phi0=-27.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000