Many nonlinear differential equations arising from practical problems may permit nontrivial multiple solutions relevant to applications, and these multiple solutions are helpful to deeply understand these practical problems and to improve some applications. Developing an efficient numerical method for finding multiple solutions is very necessary due to the nonlinearity and multiple solutions of these equations. Moreover, providing an efficient iteration plays an important role in successfully obtaining multiple solutions with fast and stable convergence. In the current paper, an efficient algorithm for finding multiple solutions of semilinear elliptic systems is proposed, where the trust region Levenberg-Marquardt method is firstly used to iterate the resulted nonlinear algebraic system. When the nonlinear term in these equations has only the first derivative, our algorithm can efficiently find multiple solutions as well. Several numerical experiments are tested to show the efficiency of our algorithm, and some solutions which have not been shown in the literature are also found and shown.
Li et al. (Sat,) studied this question.