File html/e3a.html, c=[
0.25 , 0.5, 0.75, 1]
AUBV (scaled)
| 0 | 0 | 0 | 0 |
| 0.845232 | -0 | 0 | -0 |
| -0.0449691 | 0.882754 | 0 | 0 |
| 1.03061 | 0.483638 | 0.157846 | 0 |
|
| 1 | 0.25 | 0.0625 | 0.015625 |
| 1 | -0.345232 | -0.172616 | -0.033481 |
| 1 | -0.087785 | -0.29777 | -0.231759 |
| 1 | -0.672094 | -0.235712 | 0.177668 |
|
| 1.03061 | 0.483638 | 0.157846 | 0 |
| 0 | 0 | 0 | 1 |
| -0.712022 | -0.63459 | -0.0665772 | 1.06962 |
| -0.707945 | -0.57053 | -0.210071 | 0.622123 |
|
| 1 | -0.672094 | -0.235712 | 0.177668 |
| 0 | 0 | 0 | 0 |
| 0 | 0.343573 | -0.0487644 | 0.512949 |
| 0 | 0.866424 | -0.00463586 | 0.0487644 |
|
Maximum coefficient scaled is 1.069616 .
irks(0,[0.25,0.5,0.75,1], [0.2679,0.5731],[ 1 0; 0.2852 1; ])
[alpha,beta,gamma,stageerror]
| 0.00320513 | 0.0234429 | -0.0151096 | 0.0100271 |
| 0 | -0 | -0 | -0.0259168 |
| 0.2679 | -0.210659 | 0.0571253 | -0.0671131 |
| 0.5731 | -0.209704 | 0.149161 | 0.00320513 |
Error estimation
phi=[-64.000000 192.000000 -192.000000 64.000000 ]
phi1=[329.071228 -397.606843 201.071228 -34.267807 ] phi0=-98.267807 gives y^{p+1}(t_n-theta*h) at theta=0.758859