Electronic International Standard Serial Number (EISSN)
It has been recently shown in  that Problem (1), for the special choice (2), admits an arbitrarily large number of positive solutions, provided that lambda is sufficiently negative. Moreover, using b as the main bifurcation parameter, some fundamental qualitative properties of the associated global bifurcation diagrams have been established. Based on them, the authors computed such bifurcation diagrams by coupling some adaptation of the classical path-following solvers with spectral methods and collocation (see ). In this paper, we complete our original program by computing the global bifurcation diagrams for a wider relevant class of weight functions a(x)'s. The numerics suggests that the analytical results of  should be true for general symmetric weight functions, whereas some of them can fail if a(x) becomes asymmetric around 0.5. In any of these circumstances, the more negative lambda, the larger the number of positive solutions of Problem (1). As an astonishing ecological consequence, facilitation in competitive environments within polluted habitat patches causes complex dynamics.
global bifurcation diagrams; isolas; path following; pseudo-spectral methods; large solutions; superlinear indefinite problems