nag_ode_ivp_bdf_zero_simple (d02ejc) Example Program Results Case 1: calculating Jacobian internally intermediate output, root-finding Calculation with tol = 1.0e-03 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 (User-supplied callback g, first invocation.) (User-supplied callback fcn, first invocation.) 2.00 0.9416 0.00003 0.0583 4.00 0.9055 0.00002 0.0945 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Calculation with tol = 1.0e-04 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 2.00 0.9416 0.00003 0.0584 4.00 0.9055 0.00002 0.0945 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Case 2: calculating Jacobian by pederv intermediate output, root-finding Calculation with tol = 1.0e-03 X Y(1) Y(2) Y(3) (User-supplied callback pederv, first invocation.) 0.00 1.0000 0.00000 0.0000 2.00 0.9416 0.00003 0.0583 4.00 0.9055 0.00002 0.0945 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Calculation with tol = 1.0e-04 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 2.00 0.9416 0.00003 0.0584 4.00 0.9055 0.00002 0.0945 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Case 3: calculating Jacobian internally no intermediate output, root-finding Calculation with tol = 1.0e-03 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Calculation with tol = 1.0e-04 Root of Y(1)-0.9 at 4.377 Solution is 0.9000 0.00002 0.1000 Case 4: calculating Jacobian internally intermediate output, no root-finding Calculation with tol = 1.0e-03 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 2.00 0.9416 0.00003 0.0583 4.00 0.9055 0.00002 0.0945 6.00 0.8793 0.00002 0.1207 8.00 0.8586 0.00002 0.1414 10.00 0.8414 0.00002 0.1586 10.00 0.8414 0.00002 0.1586 Calculation with tol = 1.0e-04 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 2.00 0.9416 0.00003 0.0584 4.00 0.9055 0.00002 0.0945 6.00 0.8793 0.00002 0.1207 8.00 0.8585 0.00002 0.1414 10.00 0.8414 0.00002 0.1586 10.00 0.8414 0.00002 0.1586 Case 5: calculating Jacobian internally no intermediate output, no root-finding (integrate to xend) Calculation with tol = 1.0e-03 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 10.00 0.8414 0.00002 0.1586 Calculation with tol = 1.0e-04 X Y(1) Y(2) Y(3) 0.00 1.0000 0.00000 0.0000 10.00 0.8414 0.00002 0.1586