File html/e5b.html, c=[
0.2 , 0.4, 0.6, 0.8, 1, 1]
AUBV (scaled)
| -0 | -0 | 0 | 0 | -0 | 0 |
| -0.0616368 | -0 | 0 | 0 | -0 | 0 |
| -0.524361 | 0.241996 | 0 | 0 | -0 | 0 |
| -0.0339762 | -1.21336 | 0.987761 | 0 | -0 | 0 |
| -0.130498 | -1.88414 | 1.2998 | 0.213985 | -0 | 0 |
| -1.39939 | -5.11795 | 1.57531 | 5.25887 | -2.731 | 0 |
|
| 1 | 0.2 | 0.04 | 0.008 | 0.0016 | 0.00032 |
| 1 | 0.461637 | 0.184655 | 0.0713964 | 0.0275724 | 0.0107331 |
| 1 | 0.882365 | 0.376148 | 0.162765 | 0.0844285 | 0.0509794 |
| 1 | 1.05957 | 0.438964 | 0.0317072 | -0.132119 | -0.156807 |
| 1 | 1.50085 | 0.657376 | 0.105413 | -0.0747516 | -0.0382968 |
| 1 | 3.41417 | 0.811561 | 0.0191855 | 1.14775 | 3.53033 |
|
| -1.39939 | -5.11795 | 1.57531 | 5.25887 | -2.731 | 0 |
| 0 | 0 | -0 | 0 | 0 | 1 |
| 5.86211 | -4.82629 | 1.75058 | -6.03705 | 5.81961 | -0.997842 |
| 2.83945 | 4.93339 | -4.62701 | -4.85666 | 4.18109 | -0.0465228 |
| 0.576809 | 4.10256 | -1.41636 | -5.33145 | 3.59682 | -0.359532 |
| 0.0531365 | 0.453402 | 0.00701614 | -0.959763 | 0.653503 | -0.119904 |
|
| 1 | 3.41417 | 0.811561 | 0.0191855 | 1.14775 | 3.53033 |
| 0 | 0 | 0 | 0 | 0 | 0 |
| 0 | -1.57113 | 0.431223 | -0.151643 | -1.38777 | -2.30848 |
| 0 | -2.42374 | -0.0285611 | 0.209494 | 0.0520968 | 1.61772 |
| 0 | -1.16885 | 0.242608 | 0.0157495 | -0.875312 | 0.12044 |
| 0 | -0.0873902 | 0.0760278 | 0.010362 | -0.292633 | 0.234594 |
|
Maximum coefficient scaled is 6.037046 .
irks(-7.6316e-15,[0.2,0.4,0.6,0.8,1,1], [-0.0695,0.1387,0.0194,0.5997],[ 0.299354 -0.057819 -0.220706 1; -0.152033 -0.291561 1 0; 1 0 0 0; 0.550099 1 0 0; ])
[alpha,beta,gamma,stageerror]
| 1.90259e-05 | 0.0495535 | -0.049355 | -0.00120209 |
| 0 | 0 | -0 | -0.00468452 |
| -0.0695 | 0.0974205 | -0.0557413 | -0.00893998 |
| 0.1387 | -0.160549 | 0.00787109 | -0.0155838 |
| 0.0194 | 0.192261 | -0.0651948 | -0.0205354 |
| 0.5997 | -0.0225329 | -0.317707 | 1.90259e-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 4465.761144 -1340.761144 ] phi0=-3125.000000 gives y^{p+1}(t_n-theta*h) at theta=0.500000