```
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
```