In this paper, we compile the fractional power series method and the Laplace transform to design a new algorithm for solving the fractional Volterra integro-differential equation. For that, we assume the Laplace power series (LPS) solution in terms of power q=1m,m∈Z+, where the fractional derivative of order α=qγ, for which γ∈Z+. This assumption will help us to write the integral, the kernel, and the nonhomogeneous terms as a LPS with the same power. The recurrence relations for finding the series coefficients can be constructed using this form. To demonstrate the algorithm’s accuracy, the residual error is defined and calculated for several values of the fractional derivative. Two strongly nonlinear examples are discussed to provide the efficiency of the algorithm. The algorithm gains powerful results for this kind of fractional problem. Under Caputo meaning of the symmetry order, the obtained results are illustrated numerically and graphically. Geometrically, the behavior of the obtained solutions declares that the changing of the fractional derivative parameter values in their domain alters the style of these solutions in a symmetric meaning, as well as indicates harmony and symmetry, which leads them to fully coincide at the value of the ordinary derivative. From these simulations, the results report that the recommended novel algorithm is a straightforward, accurate, and superb tool to generate analytic-approximate solutions for integral and integro-differential equations of fractional order.