Abstract:Recently, fractional differential equations (FDEs) have attracted much more attention in modeling real-life problems. Since most FDEs do not have exact solutions, numerical solution methods are used commonly. Therefore, in this study, we have demonstrated a novel approximate-analytical solution method, which is called the Laplace homotopy analysis method (LHAM) using the Caputo-Fabrizio (CF) fractional derivative operator. The recommended method is obtained by combining Laplace transform (LT) and the homotopy analysis method (HAM). We have used the fractional operator suggested by Caputo and Fabrizio in 2015 based on the exponential kernel. We have considered the LHAM with this derivative in order to obtain the solutions of the fractional Black-Scholes equations (FBSEs) with the initial conditions. In addition to this, the convergence and stability analysis of the model have been constructed. According to the results of this study, it can be concluded that the LHAM in the sense of the CF fractional derivative is an effective and accurate method, which is computable in the series easily in a short time.