nag_pde_parab_1d_cd_ode (d03plc) Example Program Results


Example 1

 Method parameters:
  Number of mesh points used =  201
  Relative tolerance used    =    2.500e-04
  Absolute tolerance used    =    1.000e-05


 Integration Results:
  Global error is less than  100 times the local error tolerance.

 Integration Statistics:
  Number of time steps           (nearest  50) =    150
  Number of function evaluations (nearest 100) =   1400
  Number of Jacobian evaluations (nearest  20) =     20
  Number of iterations           (nearest 100) =    400


Example 2

 Problem parameters and initial conditions:
  gamma          = 1.400
      e(x<0.5,0) = 2.500      e(x>0.5,0) = 0.250
    rho(x<0.5,0) = 1.000    rho(x>0.5,0) = 0.125

 Method parameters:
  Number of mesh points used =  141
  Relative tolerance used    =    5.000e-04
  Absolute tolerance used    =    5.000e-03

 Solution
   t        x        d        v        p
 0.100   0.0000   1.0000   0.0000   1.0000
         0.1000   1.0000  -0.0000   1.0000
         0.2000   1.0000  -0.0000   1.0000
         0.3000   1.0000  -0.0000   1.0000
         0.4000   0.8668   0.1665   0.8188
         0.5000   0.4299   0.9182   0.3071
         0.6000   0.2969   0.9274   0.3028
         0.7000   0.1250   0.0000   0.1000
         0.8000   0.1250  -0.0000   0.1000
         0.9000   0.1250  -0.0000   0.1000
         1.0000   0.1250   0.0000   0.1000

 0.200   0.0000   1.0000   0.0000   1.0000
         0.1000   1.0000  -0.0000   1.0000
         0.2000   1.0000  -0.0000   1.0000
         0.3000   0.8718   0.1601   0.8253
         0.4000   0.6113   0.5543   0.5022
         0.5000   0.4245   0.9314   0.3014
         0.6000   0.4259   0.9277   0.3030
         0.7000   0.2772   0.9272   0.3031
         0.8000   0.2657   0.9276   0.3032
         0.9000   0.1250  -0.0000   0.1000
         1.0000   0.1250   0.0000   0.1000


 Integration Statistics:
  Number of time steps           (nearest  50) =    150
  Number of function evaluations (nearest  50) =    400
  Number of Jacobian evaluations (nearest   1) =      1
  Number of iterations           (nearest   1) =      2