This paper brings together methods to solve and/or characterize solutions of some problems of mathematical physics equations involving p-Laplacian and p-pseudo-Laplacian. Using surjectivity or variational approaches, one may obtain or characterize weak solutions for Dirichlet or Newmann problems for these important operators. This article details three ways to use surjectivity results for a special type of operator involving the duality mapping and a Nemytskii operator, three methods starting from Ekeland’s variational principle and, lastly, one with a generalized variational principle to solve or describe the above-mentioned solutions. The relevance of these operators and the possibility of their involvement in the modeling of an important class of real phenomena determined the author to group these seven procedures together, presented in detail, followed by many applications, accompanied by a general overview of specialty domains. The use of certain variational methods facilitates the complete solution of the problem via appropriate numerical methods and computational algorithms. The exposure of the sequence of theoretical results, together with their demonstration in as much detail as possible has been fulfilled as an opportunity for the complete development of these topics.