An $R_2$-Reversible Homoclinic Solution.

make third

Copies the files rev.c.3 and rev.dat.3 to rev.f and rev.dat, and runs them with the constants stored in c.rev.3 and h.rev.3. The orbit contained in the data file is a ``multi-pulse'' homoclinic solution for $P=1.6$, with truncation (half-)interval PAR(11) = 47.4464189. which is reversible under $R_2$. This reversibility is specified in h.rev.1 via NREV=1, (IREV(I), I=1,NDIM) = 1 0 1 0. The output

  BR    PT  TY LAB    PAR(1)        L2-NORM       MAX U(1)   ...
   1    15  UZ   2  1.700000E+00  3.836401E-01  4.890015E-01  
   1    16  LP   3  1.711574E+00  3.922135E-01  5.442385E-01  
   1    19  UZ   4  1.600000E+00  4.329404E-01  7.769491E-01  
   1    31  UZ   5  1.000000E+00  4.808488E-01  1.083298E+00  
   1    86  UZ   6 -9.664802E-10  5.158463E-01  1.258650E+00

contains the label of a limit point (ILP was set to 1 in c.rev.3, which corresponds to a ``coalescence'' of two reversible homoclinic orbits. The two solutions on either side of this limit point are displayed in Figure 26.3. The computation ends in a no-convergence point. The solution here is depicted in Figure 26.4. The lack of convergence is due to the large peak and trough of the solution rapidly moving to the left as $P \to -2$ (cf. ChSp:93).

Figure 26.3: Two $R_2$-reversible homoclinic orbits at $P=1.6$ corresponding to labels 1 (smaller amplitude) and 5 (larger amplitude)

Figure 26.4: An $R_2$-reversible homoclinic orbit at label 8

Continuing from the initial solution in the other parameter direction

make fourth

we obtain the output

  BR    PT  TY LAB    PAR(1)        L2-NORM       MAX U(1)   ...
   1     7  UZ   8  1.600000E+00  3.701709E-01  3.836833E-01  
   1    33  UZ   9  9.999980E-01  3.614405E-01  1.775035E-01  
   1    93  UZ  10 -7.819855E-06  3.713007E-01  4.698309E-02

which again ends at a no convergence error for similar reasons.