File html/i5b.html, c=[
0.2 , 0.4, 0.6, 0.8, 1, 1]
AUBV (scaled)
| 0.333333 | 0 | -0 | -0 | 0 | 0 |
| -0.0499118 | 0.333333 | -0 | -0 | 0 | 0 |
| 0.565215 | -0.349134 | 0.333333 | -0 | 0 | 0 |
| 0.535747 | -1.71315 | 0.883819 | 0.333333 | 0 | 0 |
| 0.973119 | -3.2072 | 1.57123 | 0.104471 | 0.333333 | 0 |
| 6.56081 | -9.45014 | 1.08954 | 6.83006 | -3.46669 | 0.333333 |
|
| 1 | -0.133333 | -0.0933333 | -0.032 | -0.00906667 | -0.00234667 |
| 1 | 0.116578 | -0.0867019 | -0.0900106 | -0.0581362 | -0.0320274 |
| 1 | 0.0505857 | 0.0132212 | -0.0442414 | -0.0871085 | -0.0980726 |
| 1 | 0.760254 | 0.202308 | -0.3245 | -0.615263 | -0.712703 |
| 1 | 1.22504 | 0.457212 | -0.474835 | -1.11493 | -1.49605 |
| 1 | -0.896918 | -0.0330515 | -0.141593 | 0.813371 | 3.12991 |
|
| 6.56081 | -9.45014 | 1.08954 | 6.83006 | -3.46669 | 0.333333 |
| 0 | 0 | -0 | 0 | 0 | 1 |
| 3.15637 | -7.8322 | 2.55608 | -3.12908 | 3.53191 | -0.174677 |
| -2.60131 | 7.66182 | -1.04218 | -8.97901 | 5.01739 | 0.283553 |
| 0.132881 | 0.0450398 | 1.75874 | -5.78241 | 3.30753 | -0.157463 |
| 1.01517 | -2.26648 | 0.998732 | -0.541188 | 0.516769 | -0.161805 |
|
| 1 | -0.896918 | -0.0330515 | -0.141593 | 0.813371 | 3.12991 |
| 0 | 0 | 0 | 0 | 0 | 0 |
| 0 | 1.8916 | 1.22798 | -0.443736 | -1.32498 | -1.05687 |
| 0 | -0.340269 | -0.0737892 | 0.0969065 | 0.207498 | 1.59972 |
| 0 | 0.695683 | 0.752049 | -0.284981 | -1.29323 | -0.0544567 |
| 0 | 0.43881 | 0.364615 | -0.13835 | -0.626674 | -0.0316581 |
|
Maximum coefficient scaled is 9.450143 .
irks(0.333333,[0.2,0.4,0.6,0.8,1,1], [-0.0575,0.0472,-0.1525,-0.1387],[ 0.0819494 0.057375 1 0; 1 0 0 0; -0.302922 0.412735 -0.0146192 1; -0.78609 1 0 0; ])
[alpha,beta,gamma,stageerror]
| -1.71468e-05 | 0.0435988 | -0.0432975 | 0.00249112 |
| -0 | 0 | 0 | 0.0021684 |
| -0.0575 | 0.100943 | -0.108918 | -0.0001842 |
| 0.0472 | 0.0669665 | -0.11792 | -0.00460993 |
| -0.1525 | 0.472942 | -0.446894 | -0.0104878 |
| -0.1387 | 1.10315 | -1.27489 | -1.71468e-05 |
Error estimation
phi=[0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 ]
phi1=[15625.000000 -31250.000000 31250.000000 -15625.000000 3570.131940 -445.131940 ] phi0=-3125.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000