nag_ode_ivp_bdf_gen (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.00000      0.00000      0.00000
(User-supplied callback g, first invocation.)
(User-supplied callback fcn, first invocation.)
    2.00      0.94163      0.00003      0.05834
    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 

 Calculation with tol = 1.0e-04
     X         Y(1)         Y(2)         Y(3)
    0.00      1.00000      0.00000      0.00000
    2.00      0.94161      0.00003      0.05836
    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 = 1.0e-03
     X         Y(1)         Y(2)         Y(3)
(User-supplied callback pederv, first invocation.)
    0.00      1.00000      0.00000      0.00000
    2.00      0.94163      0.00003      0.05834
    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 

 Calculation with tol = 1.0e-04
     X         Y(1)         Y(2)         Y(3)
    0.00      1.00000      0.00000      0.00000
    2.00      0.94161      0.00003      0.05836
    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 = 1.0e-03
  Root of Y(1)-0.9 at 4.377
  Solution is 0.90000 0.00002 0.09998 

 Calculation with tol = 1.0e-04
  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 = 1.0e-03
     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.09446
    6.00      0.87928      0.00002      0.12070
    8.00      0.85859      0.00002      0.14139
   10.00      0.84143      0.00002      0.15855
   10.00      0.84143      0.00002      0.15855

 Calculation with tol = 1.0e-04
     X         Y(1)         Y(2)         Y(3)
    0.00      1.00000      0.00000      0.00000
    2.00      0.94161      0.00003      0.05836
    4.00      0.90551      0.00002      0.09446
    6.00      0.87926      0.00002      0.12072
    8.00      0.85854      0.00002      0.14144
   10.00      0.84136      0.00002      0.15862
   10.00      0.84136      0.00002      0.15862

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.00000      0.00000      0.00000
   10.00      0.84143      0.00002      0.15855

 Calculation with tol = 1.0e-04
     X         Y(1)         Y(2)         Y(3)
    0.00      1.00000      0.00000      0.00000
   10.00      0.84136      0.00002      0.15862