with boundary conditions u(0,t)=u(1,t)=A and v(0,t)=v(1,t)=B/A. All parameters are given fixed values for which a stable periodic solution is known to exist.
The continuation parameter is the independent time variable,
namely PAR(14).
The AUTO-constants DS, DSMIN, and DSMAX
then control the step size
in space-time, here consisting of PAR(14) and (u(x),v(x)).
Initial data at time zero are
and
.
Note that in the subroutine STPNT the space derivatives of u and v
must also be provided;
see the equations-file bru.f.
Euler time integration is only first order accurate, so that the time step must be sufficiently small to ensure correct results. This option has been added only as a convenience, and should generally be used only to locate stationary states. Indeed, in the case of the asymptotic periodic state of this demo, the number of required steps is very large and use of a better time integrator is advisable.
| COMMAND | ACTION |
| mkdir bru | create an empty work directory |
| cd bru | change directory |
| @dm bru | copy the demo files to the work directory |
| cp r.bru.1 r.bru | get the constants-file |
| @r bru | time integration |
| @sv bru | save output-files as p.bru, q.bru, d.bru |