nag_ode_bvp_coll_nlin_solve (d02tlc) Example Program Results
lrcomm = 3012, licomm = 127
Tolerance = 1.0e-03 r = 1.000e+06
(User-supplied callback guess, first invocation.)
(User-supplied callback gafun, first invocation.)
(User-supplied callback gajac, first invocation.)
(User-supplied callback gbfun, first invocation.)
(User-supplied callback gbjac, first invocation.)
(User-supplied callback ffun, first invocation.)
(User-supplied callback fjac, first invocation.)
Used a mesh of 21 points
Maximum error = 8.05e-08 in interval 5 for component 3
Mesh points:
1(1) 0.0000e+00 2(3) 5.0000e-02 3(2) 1.0000e-01 4(3) 1.5000e-01
5(2) 2.0000e-01 6(3) 2.5000e-01 7(2) 3.0000e-01 8(3) 3.5000e-01
9(2) 4.0000e-01 10(3) 4.5000e-01 11(2) 5.0000e-01 12(3) 5.5000e-01
13(2) 6.0000e-01 14(3) 6.5000e-01 15(2) 7.0000e-01 16(3) 7.5000e-01
17(2) 8.0000e-01 18(3) 8.5000e-01 19(2) 9.0000e-01 20(3) 9.5000e-01
21(1) 1.0000e+00
x f f' g
0.000 0.0000 0.0000 1.0000
0.050 0.0070 0.1805 0.4416
0.100 0.0141 0.0977 0.1886
0.150 0.0171 0.0252 0.0952
0.200 0.0172 -0.0165 0.0595
0.250 0.0157 -0.0400 0.0427
0.300 0.0133 -0.0540 0.0322
0.350 0.0104 -0.0628 0.0236
0.400 0.0071 -0.0683 0.0156
0.450 0.0036 -0.0714 0.0078
0.500 -0.0000 -0.0724 -0.0000
0.550 -0.0036 -0.0714 -0.0078
0.600 -0.0071 -0.0683 -0.0156
0.650 -0.0104 -0.0628 -0.0236
0.700 -0.0133 -0.0540 -0.0322
0.750 -0.0157 -0.0400 -0.0427
0.800 -0.0172 -0.0165 -0.0595
0.850 -0.0171 0.0252 -0.0952
0.900 -0.0141 0.0977 -0.1886
0.950 -0.0070 0.1805 -0.4416
1.000 -0.0000 0.0000 -1.0000
Tolerance = 1.0e-03 r = 1.000e+08
Used a mesh of 21 points
Maximum error = 1.97e-06 in interval 16 for component 3
Mesh points:
1(1) 0.0000e+00 2(3) 1.7442e-02 3(2) 3.4884e-02 4(3) 5.1845e-02
5(2) 6.8806e-02 6(3) 8.6127e-02 7(2) 1.0345e-01 8(3) 1.3538e-01
9(2) 1.6731e-01 10(3) 2.3097e-01 11(2) 2.9464e-01 12(3) 4.6995e-01
13(2) 6.4526e-01 14(3) 7.4273e-01 15(2) 8.4020e-01 16(3) 8.8093e-01
17(2) 9.2167e-01 18(3) 9.4424e-01 19(2) 9.6681e-01 20(3) 9.8341e-01
21(1) 1.0000e+00
x f f' g
0.000 0.0000 0.0000 1.0000
0.017 0.0025 0.1718 0.3954
0.035 0.0047 0.0836 0.1403
0.052 0.0056 0.0271 0.0527
0.069 0.0058 0.0026 0.0214
0.086 0.0057 -0.0074 0.0097
0.103 0.0056 -0.0113 0.0052
0.135 0.0052 -0.0135 0.0027
0.167 0.0047 -0.0140 0.0021
0.231 0.0038 -0.0142 0.0016
0.295 0.0029 -0.0142 0.0012
0.470 0.0004 -0.0143 0.0002
0.645 -0.0021 -0.0143 -0.0008
0.743 -0.0035 -0.0142 -0.0014
0.840 -0.0048 -0.0139 -0.0022
0.881 -0.0054 -0.0128 -0.0035
0.922 -0.0058 -0.0040 -0.0136
0.944 -0.0057 0.0194 -0.0424
0.967 -0.0046 0.0916 -0.1551
0.983 -0.0023 0.1747 -0.4155
1.000 0.0000 0.0000 -1.0000