“…To discretize (1.1) we have used two methods. To compute the small positive solutions bifurcating from u = 0 we implemented a pseudo-spectral method combining a trigonometric spectral method with collocation at equidistant points, as in R. Gómez-Reñasco and J. López-Gómez [28,31], J. López-Gómez, J. C. Eilbeck, K. Duncan and M. Molina-Meyer [42], J. López-Gómez and M. Molina-Meyer [44,45,46], J. López-Gómez, M. Molina-Meyer and A. Tellini [47], J. López-Gómez, M. Molina-Meyer and P. H. Rabinowitz [48], and M. Fencl and J. López-Gómez [25]. This gives high accuracy at a rather reasonable computational cost (see, e.g., Canuto, Hussaini, Quarteroni and Zang [13]).…”