File html/e4a.html, c=[
0.2 , 0.4, 0.6, 0.8, 1]
AUBV (scaled)
| 0 | 0 | 0 | 0 | 0 |
| 0.198404 | 0 | 0 | 0 | 0 |
| -1.71421 | 1.17656 | 0 | 0 | -0 |
| 0.254936 | 0.12468 | 0.371168 | 0 | -0 |
| -1.63318 | 0.988223 | 0.187848 | 0.321139 | -0 |
|
| 1 | 0.2 | 0.04 | 0.008 | 0.0016 |
| 1 | 0.201596 | 0.0806386 | 0.0401916 | 0.0192511 |
| 1 | 1.13766 | 0.104439 | -0.143042 | -0.116744 |
| 1 | 0.0492161 | -0.00712016 | 0.0206996 | 0.0488346 |
| 1 | 1.13598 | 0.123456 | -0.0978265 | -0.0207155 |
|
| -1.63318 | 0.988223 | 0.187848 | 0.321139 | 0 |
| 0 | 0 | 0 | 0 | 1 |
| -0.749657 | -0.120571 | -1.31185 | -0.390745 | 1.81176 |
| 0.347197 | 0.0452513 | -1.17491 | -1.44482 | 1.78275 |
| 0.817348 | -0.194357 | -0.191474 | -0.570419 | 0.490196 |
|
| 1 | 1.13598 | 0.123456 | -0.0978265 | -0.0207155 |
| 0 | -0 | -0 | -0 | -0 |
| 0 | 0.761055 | -0.0278022 | -0.120437 | 0.741477 |
| 0 | 0.444535 | -0.018969 | -0.368666 | 0.820434 |
| 0 | -0.351294 | -0.00940577 | -0.173381 | 0.396468 |
|
Maximum coefficient scaled is 1.811765 .
irks(0,[0.2,0.4,0.6,0.8,1], [0.085,0.308,0.9092],[ 1 0 0; 0.150324 0.539995 1; 0.547644 1 0; ])
[alpha,beta,gamma,stageerror]
| 0.000333333 | -0.00119969 | 0.00258858 | 0.00217395 |
| -0 | 0 | 0 | 0.00629171 |
| 0.085 | -0.0895058 | 0.0122238 | -0.00781957 |
| 0.308 | -0.219578 | 0.0245447 | 0.00318636 |
| 0.9092 | -0.503946 | -0.0146989 | 0.000333333 |
Error estimation
phi=[625.000000 -2500.000000 3750.000000 -2500.000000 625.000000 ]
phi1=[-1138.952734 1027.905467 222.094533 -736.047266 272.209453 ] phi0=352.790547 gives y^{p+1}(t_n-theta*h) at theta=0.512893