BR PT TY LAB PAR(1) ... PAR(2) PAR(29) PAR(30) 1 10 11 7.114523E+00 ... 7.081751E-02 -4.649861E-01 3.183429E-03 1 20 12 9.176810E+00 ... 7.678731E-02 -4.684912E-01 1.609294E-02 1 30 13 1.210834E+01 ... 8.543468E-02 -4.718871E-01 3.069638E-02 1 40 EP 14 1.503788E+01 ... 9.428036E-02 -4.743794E-01 4.144558E-02The fact that PAR(29) and PAR(30) do not change sign indicates that there are no further non-hyperbolic equilibria along this family. Note that restarting in the opposite direction with IRS=11, DS=-0.02
BR PT TY LAB PAR(1) ... PAR(2) PAR(29) PAR(30) 1 10 UZ 15 6.610459E+00 ... 6.932482E-02 -4.636603E-01 1.725013E-09Note that the values of PAR(1) and PAR(2) differ from that at label 4 only in the sixth significant figure.
Actually, the program runs further and eventually computes the point and the whole lower family of emanating from it, however, the solutions between and should be considered as spurious22.1, therefore we do not save these data. The reliable way to compute the lower family of is to restart computation of saddle homoclinic orbits in the other direction from the point
BR PT TY LAB PAR(1) ... PAR(2) PAR(29) PAR(30) 1 10 15 4.966429E+00 ... 6.298418E-02 -4.382426E-01 4.946824E-03 1 20 16 4.925379E+00 ... 7.961214E-02 -3.399102E-01 3.288447E-02 1 30 17 7.092267E+00 ... 1.587114E-01 -1.692842E-01 3.876291E-02 1 40 EP 18 1.101819E+01 ... 2.809825E-01 -3.482651E-02 2.104384E-02The data are appended to the stored results in b.1, s.1 and d.1. One could now display all data using the AUTO command @p 1 to reproduce the curve shown in Figure 22.1.