File html/e4b.html, c=[
0.25 , 0.5, 0.75, 1, 1]
AUBV (scaled)
| 0 | 0 | 0 | 0 | 0 |
| 0.320112 | 0 | -0 | 0 | 0 |
| -0.475669 | 0.537635 | 0 | -0 | -0 |
| -0.456266 | 0.421881 | 0.247582 | -0 | 0 |
| -0.463912 | -0.125869 | -0.410284 | 0.701687 | -0 |
|
| 1 | 0.25 | 0.0625 | 0.015625 | 0.00390625 |
| 1 | 0.179888 | 0.0899439 | 0.064979 | 0.042493 |
| 1 | 0.688034 | 0.2627 | 0.107837 | 0.0773183 |
| 1 | 0.786803 | 0.434879 | 0.351344 | 0.399781 |
| 1 | 1.29838 | 0.569877 | -0.231322 | -1.02247 |
|
| -0.463912 | -0.125869 | -0.410284 | 0.701687 | 0 |
| 0 | 0 | 0 | 0 | 1 |
| -1.39531 | -1.09018 | -0.513546 | 1.53777 | 0.0820588 |
| 0.934241 | -1.06952 | -0.940365 | 0.603287 | 0.735 |
| 0.413923 | -0.177612 | -0.170126 | -0.0282253 | 0.245098 |
|
| 1 | 1.29838 | 0.569877 | -0.231322 | -1.02247 |
| 0 | -0 | -0 | -0 | -0 |
| 0 | 1.37921 | 0.318494 | 0.0863759 | 1.01959 |
| 0 | -0.262646 | -0.66363 | -0.801027 | 0.710088 |
| 0 | -0.283058 | -0.207906 | -0.307933 | 0.482532 |
|
Maximum coefficient scaled is 1.537770 .
irks(0,[0.25,0.5,0.75,1,1], [0.17,0.0279,0.7497],[ -0.324197 0.435914 1; 1 0 0; 0.876408 1 0; ])
[alpha,beta,gamma,stageerror]
| 0.000333333 | -0.0741312 | 0.0755201 | 0.00551532 |
| -0 | 0 | 0 | 0.00948308 |
| 0.17 | -0.228215 | 0.100766 | 0.025901 |
| 0.0279 | 0.286427 | -0.137343 | 0.0551316 |
| 0.7497 | -0.130257 | -0.403893 | 0.000333333 |
Error estimation
phi=[0.000000 0.000000 0.000000 0.000000 0.000000 ]
phi1=[-1024.000000 1536.000000 -1024.000000 319.700178 -63.700178 ] phi0=256.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000