Reaction-diffusion equations have several applications in the field of population dynamics and some of them are characterized by the presence of a factor which describes different types of food sources in a heterogeneous habitat. In this context, to study persistence or extinction of populations it is relevant the search of nontrivial steady states. Our paper focuses on a one-dimensional model given by a parameter-dependent equation of the form u + λa + (t) − µa − (t) g(u) = 0, where g : [0, 1] → R is a continuous function such that g(0) = g(1) = 0, g(s) > 0 for every 0 < s < 1 and lim s→0 + g(s)/s = 0, and the weight a(t) has two positive humps separated by a negative one. In this manner, we consider bounded habitats which include two favorable food sources and an unfavorable one. We deal with various boundary conditions, including the Dirichlet and Neumann ones, and we prove the existence of eight positive solutions when λ and µ are positive and sufficiently large. Throughout the paper, numerical simulations are exploited to discuss the results and to explore some open problems. * Work partially supported by the Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). Progetto di Ricerca 2017: "Problemi differenziali con peso indefinito: tra metodi topologici e aspetti dinamici". Guglielmo Feltrin is supported by the Belgian F.R.S.-FNRS -Fonds de la Recherche Scientifique, Chargé de recherches project: "Quantitative and qualitative properties of positive solutions to indefinite problems arising from population genetic models: topological methods and numerical analysis", and partially by the project "Existence and asymptotic behavior of solutions to systems of semilinear elliptic partial differential equations" (T.1110.14).