File html/e6a.html, c=[
0.142857 , 0.285714, 0.428571, 0.571429, 0.714286, 0.857143, 1]
AUBV (scaled)
| 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| -0.00124917 | 0 | -0 | 0 | 0 | -0 | 0 |
| -2.2183 | 1.33945 | -0 | 0 | 0 | 0 | 0 |
| -16.5919 | 10.8788 | -0.720646 | 0 | -0 | 0 | -0 |
| -96.8 | 68.4419 | -8.8095 | 0.468064 | -0 | 0 | -0 |
| -270.543 | 196.026 | -29.6931 | 2.51207 | -0.103536 | 0 | 0 |
| -125.503 | 246.354 | -173.91 | -25.281 | 84.0171 | -27.144 | 0 |
|
| 1 | 0.142857 | 0.0204082 | 0.00291545 | 0.000416493 | 5.9499e-05 | 8.49986e-06 |
| 1 | 0.286963 | 0.0819896 | 0.0234001 | 0.00667846 | 0.00190657 | 0.000544437 |
| 1 | 1.30743 | 0.0520771 | -0.113496 | -0.0653575 | -0.0255518 | -0.00831325 |
| 1 | 7.00524 | -0.531665 | -1.06468 | -0.487907 | -0.145438 | -0.0210228 |
| 1 | 37.4138 | -3.92624 | -6.07461 | -2.57158 | -0.656475 | -0.0213862 |
| 1 | 102.658 | -11.2541 | -16.7536 | -6.96789 | -1.70125 | 0.0268387 |
| 1 | 22.4666 | 0.550622 | 0.175214 | -0.992009 | -0.227857 | 4.18015 |
|
| -125.503 | 246.354 | -173.91 | -25.281 | 84.0171 | -27.144 | 0 |
| 0 | 0 | 0 | 0 | 0 | 0 | 1 |
| 231.098 | -203.472 | -98.8397 | 325.94 | -253.758 | 71.7422 | -0.853333 |
| -94.2675 | 126.213 | -214.838 | 307.657 | -218.829 | 61.2421 | -0.395556 |
| 79.2598 | -67.4244 | -83.6532 | 234.618 | -195.264 | 56.4865 | -0.395278 |
| 81.5029 | -134.323 | 50.9727 | 84.2501 | -93.8516 | 28.1104 | -0.281278 |
| 13.5374 | -18.7307 | 3.47291 | 15.3739 | -14.6523 | 4.28942 | -0.0925926 |
|
| 1 | 22.4666 | 0.550622 | 0.175214 | -0.992009 | -0.227857 | 4.18015 |
| 0 | -0 | -0 | -0 | -0 | -0 | -0 |
| 0 | -71.8567 | 4.69066 | 6.69563 | 2.7507 | 0.128751 | -4.27577 |
| 0 | 33.2177 | -4.23209 | -5.98986 | -2.34254 | -0.277927 | 1.25425 |
| 0 | -23.6282 | 2.35455 | 3.48162 | 1.53223 | -0.215688 | 0.634521 |
| 0 | -16.3804 | -0.0603297 | -0.175634 | -0.222779 | -0.51521 | 0.64 |
| 0 | -3.19796 | 0.0523013 | 0.0351006 | -0.0536804 | -0.228935 | 0.282179 |
|
Maximum coefficient scaled is 325.939978 .
irks(0,[0.142857,0.285714,0.428571,0.571429,0.714286,0.857143,1], [-0.015,0.0256,0.0356,0.1423,0.5063],[ 0.16553 1 0 0 0; -0.447468 0.272903 -0.0456143 1 0; 1 0 0 0 0; -0.149893 0.451252 0.342036 0.707133 1; 0.604852 0.068652 1 0 0; ])
[alpha,beta,gamma,stageerror]
| 1.5623e-06 | 0.005359 | -0.0053342 | -0.000139927 |
| -0 | 0 | 0 | -0.0005025 |
| -0.015 | 0.00868082 | -0.00108632 | -0.00100838 |
| 0.0256 | -0.0143301 | -0.0136312 | -0.00146399 |
| 0.0356 | -0.0562064 | 0.0637822 | -0.00105814 |
| 0.1423 | 0.0442923 | -0.0807388 | 0.000664819 |
| 0.5063 | 0.118929 | -0.232968 | 1.5623e-06 |
Error estimation
phi=[117649.000000 -705894.000000 1764735.000000 -2352980.000000 1764735.000000 -705894.000000 117649.000000 ]
phi1=[-1604199.163201 4459650.489603 -6844505.816005 6256260.816005 -3400809.489603 1015954.163201 -128329.309029 ] phi0=245978.309029 gives y^{p+1}(t_n-theta*h) at theta=0.727254