D02EJF Example Program Results
Case 1: calculating Jacobian internally,
intermediate output, root-finding
Calculation with TOL = 0.1E-02
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94163 0.00003 0.05834
4.00 0.90551 0.00002 0.09447
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Calculation with TOL = 0.1E-03
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94161 0.00003 0.05837
4.00 0.90551 0.00002 0.09446
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Case 2: calculating Jacobian by PEDERV,
intermediate output, root-finding
Calculation with TOL = 0.1E-02
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94163 0.00003 0.05834
4.00 0.90551 0.00002 0.09447
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Calculation with TOL = 0.1E-03
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94161 0.00003 0.05837
4.00 0.90551 0.00002 0.09446
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Case 3: calculating Jacobian internally,
no intermediate output, root-finding
Calculation with TOL = 0.1E-02
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Calculation with TOL = 0.1E-03
Root of Y(1)-0.9 at 4.377
Solution is 0.90000 0.00002 0.09998
Case 4: calculating Jacobian internally,
intermediate output, no root-finding
Calculation with TOL = 0.1E-02
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94163 0.00003 0.05834
4.00 0.90551 0.00002 0.09447
6.00 0.87930 0.00002 0.12068
8.00 0.85858 0.00002 0.14140
10.00 0.84136 0.00002 0.15862
Calculation with TOL = 0.1E-03
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
2.00 0.94161 0.00003 0.05837
4.00 0.90551 0.00002 0.09446
6.00 0.87926 0.00002 0.12072
8.00 0.85854 0.00002 0.14145
10.00 0.84136 0.00002 0.15863
Case 5: calculating Jacobian internally,
no intermediate output, no root-finding (integrate to XEND)
Calculation with TOL = 0.1E-02
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
10.00 0.84136 0.00002 0.15862
Calculation with TOL = 0.1E-03
X Y(1) Y(2) Y(3)
0.00 1.00000 0.00000 0.00000
10.00 0.84136 0.00002 0.15863