File html/e3b.html, c=[
0.333333 , 0.666667, 1, 1]
AUBV (scaled)
| 0 | 0 | 0 | 0 |
| 0.49229 | -0 | 0 | -0 |
| 1.15715 | -0.406894 | 0 | -0 |
| 1.15718 | 1.15681 | -0.446374 | 0 |
|
| 1 | 0.333333 | 0.111111 | 0.037037 |
| 1 | 0.174376 | 0.116251 | 0.132199 |
| 1 | 0.24974 | 0.77109 | 1.15681 |
| 1 | -0.867611 | -0.421116 | 0.410983 |
|
| 1.15718 | 1.15681 | -0.446374 | 0 |
| 0 | 0 | 0 | 1 |
| 0.128185 | -0.307297 | 0.00256208 | 0.747775 |
| -0.531814 | -0.640847 | 0.24758 | 0.370865 |
|
| 1 | -0.867611 | -0.421116 | 0.410983 |
| 0 | -0 | -0 | -0 |
| 0 | -0.571225 | -0.176401 | 1.11599 |
| 0 | 0.554217 | -0.0278831 | 0.176401 |
|
Maximum coefficient scaled is 1.157175 .
irks(0,[0.333333,0.666667,1,1], [0.4494,0.6721],[ 1 0; 0.4742 1; ])
[alpha,beta,gamma,stageerror]
| 0.00320513 | 0.077124 | -0.0687907 | 0.0296298 |
| -0 | 0 | -0 | 0.0461217 |
| 0.4494 | -0.389452 | 0.0874102 | 0.357463 |
| 0.6721 | -0.316714 | 0.314396 | 0.00320513 |
Error estimation
phi=[0.000000 0.000000 0.000000 0.000000 ]
phi1=[81.000000 -81.000000 3.526543 23.473457 ] phi0=-27.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000