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