Predictor Laplace Fractional Power Series Method for Finding Multiple Solutions of Fractional Boundary Value Problems

This research focuses on finding multiple solutions (MSs) to nonlinear fractional boundary value problems (BVPs) through a new development, namely the predictor Laplace fractional power series method. This method predicts the missing initial values by applying boundary or force conditions. This research provides a set of theorems necessary for deriving the recurrence relations to find the series terms. Several examples demonstrate the efficacy, convergence, and accuracy of the algorithm. Under Caputo’s definition of the fractional derivative with symmetric order, the obtained results are visualized numerically and graphically. The behavior of the generated solutions indicates that altering the fractional derivative parameters within their domain symmetrically changes these solutions, ultimately aligning them with the standard derivative. The results are compared with the homotopy analysis method and are presented in various figures and tables.