“…where we used the subsolution-supersolution method (see previous studies [7][8][9][10][11][12][13][14][15][16][17][18] ), Galerkin approach, and a variational method to prove the existence of positive solutions in different cases depending on the parameters p and . We chose the logarithmic nonlinearity of source terms, because it appears in several branches of physics such as inflationary cosmology, nuclear physics, optics, and geophysics (see previous studies 2,4,12,[19][20][21][22][23] ).…”