Numerical solution for the time-fractional Fokker–Planck equation via shifted Chebyshev polynomials of the fourth kind

*Correspondence: dlsuthar@gmail.com 1Department of Mathematics, Wollo University, Dessie, Ethiopia Abstract This paper provides a numerical approach for solving the time-fractional Fokker–Planck equation (FFPE). The authors use the shifted Chebyshev collocation method and the finite difference method (FDM) to present the fractional Fokker–Planck equation into systems of nonlinear equations; the Newton–Raphson method is used to produce approximate results for the nonlinear systems. The results obtained from the FFPE demonstrate the simplicity and efficiency of the proposed method.


Introduction
Applications of fractional differential equation (FDE) in science and engineering are becoming vibrant, particularly in the fields of physics, finance, viscoelasticity, chemistry, and fluid mechanics. For this reason many authors are attracted to knowing the properties of FDEs (for details, see [2, 13, 15, 16, 18-20, 24, 25, 28]) and vast applications in modeling and engineering fields [1,6,11,12]. Also see the literature [3,5,21,23] for further applications of FDE in different disciplines.
The Fokker-Planck equation (FPE) was studied for the first time by Risken [22]. In addition, several researchers [4,17] have worked on fractional Fokker-Planck equations among mathematical models developed in physical and biological sciences. It results from a diffusion estimation of certain stochastic processes that have been re-enacted as Markovian and continuous. It is a generalized diffusion equation that governs the evolution of the probability density in time. In most cases, fractional differential equations (FDEs) cannot be solved using the exact methods [8][9][10], that is why the recent research has used the properties of shifted Chebyshev polynomials [14,27] and the finite difference method (FDM) to simplify the fractional initial value problem (IVP) to a set of nonlinear equations. In this paper, we are limited to the computational solution for the two-variable time-fractional Fokker-Planck equation [17] as follows: (1.1) where u = u(x, t), α is a parameter that defines the fractional derivative order (0 < α ≤ 1). A i (x, t, u) and B i (x, t, u) are arbitrary constants and n is an integer such that n = 0. In recent decades, Chebyshev polynomials have been among the most useful approximations because of their suitability for solving complex problems in science and engineering that can be expressed in integral equations, ordinary differential equations, and partial differential equations in integer and fractional orders. In this article, we use the fourthkind Chebyshev polynomials with some important properties and their analytical form. Here, the approximate solution of FFPE is computed with the help of the shifted Chebyshev collocation method and the finite difference scheme. The shifted Chebyshev collocation method is used to reduce FFPE into systems of nonlinear differential equations, while FDM can be used for rewriting these systems into systems of nonlinear equations, and hence the Newton-Raphson method is used to arrive at the required approximate solution.
The organization of the paper is as follows. In Sect. 2, we provide the basic properties of Chebyshev polynomials of the fourth kind. In Sect. 3, the shifted Chebyshev collocation method and the finite difference method (FDM) are implemented to solve the fractional Fokker-Planck equation. In Sect. 4, conclusions are provided.

Preliminaries
This section discusses mathematical description, fractional derivative notation, and some essential properties for the fourth-kind Chebyshev polynomials.

Definition 2.1
In the Caputo sense the fractional derivative of f (z) is defined as in [19]: where z > 0, m ∈ C, and m -1 < α < m.
The sequential property of the Caputo fractional derivative emerges analogous to the integer order differentiation For the Caputo fractional derivative, we have Γ ( -α+1) z -α , for ∈ N 0 and ≥ α , 0, for ∈ N 0 and ≥ α , the notation α is used in the sense of a ceiling function which gives an integer smaller than or equal to α.
We recall from Mason and Handscomb [14] that the fourth-kind Chebyshev polynomials W n (t) are orthogonal degree polynomials of n in t based on [-1, 1]. Hence, for t = cos θ and θ ∈ [0, π], W n (t) = sin(n + 1 2 )θ sin θ 2 , (2.5) and it can be extracted from the special case β = -α = 1 2 of the Jacobi polynomial P (α,β) n (t) explicitly: and W n (t) is the orthogonal polynomial on [-1, 1] w.r.t. the inner product as follows: π, n = m, where W n (t) has a weight function 1-t 1+t . The Chebyshev polynomials of the fourth kind W n (t) can be generated using the recurrence relation W n+1 (t) = 2tW n (t) -W n-1 (t), n = 1, 2, . . . (2.9) starting with the values W 0 (t) = 1, W 1 (t) = 2t + 1, W 2 (t) = 4t 2 + 2t -1. The analytical form of the fourth-kind Chebyshev polynomial W n (t) of degree n can be expressed using (2.6) and the properties of Jacobi polynomial in the following form: . However, in this paper, we are only concerned with the fourth-kind shifted Chebyshev polynomial of degree n symbolized by W * n (t): These shifted polynomials of the fourth kind are orthogonal at the support interval [0, 1] under the inner product shown as follows: (2.12) On the right-hand side of (2.12), 1-t t is the weight function of W * n (t) and is normalized by W * n (1) = 1. The fourth-kind shifted Chebyshev polynomial is obtained by means of a recurrence connection: starting with the values W * 0 (t) = 1, W * 1 (t) = 4t -1, W * 2 (t) = 16t 2 -12t + 1. The possible expression of the analytic form for the fourth-kind shifted Chebyshev polynomials W * n (t) of degree n is shown as follows: In a spectral method, it is possible to expand g(t), the square integrable function, in [0, 1] and define it by an infinite W * n (t) expansion as follows: where c i are constants. Now, we can estimate a number of the coefficients c i , and then g(t) can be approximated by a finite sum of terms (m + 1) such as where the coefficients c i (i = 0, 1, 2, . . . , m) are given by

Main results
The approximate formula for the function g m (t) given in (2.16) is presented in the following theorem.

Theorem 3.1 Let g m (t) be an approximate function in terms of the fourth-kind shifted
Chebyshev polynomials given by (2.16). Suppose that α > 0, we get Proof Using the definition of approximated function g m (t) given in (2.16) and the Caputo fractional differentiation properties given in (2.2), we obtain Applying equations (2.3) and (2.4), we get

Numerical examples
In this section, we present numerical examples to show the efficacy and validity of the proposed method and to compare it with the existing method.
Example 4.1 Let us consider the time-fractional Fokker-Planck equation with the initial condition The exact solution to equation (4.1) for the non-fractional case at α = 1 is u( In order to use the shifted Chebyshev polynomials of the fourth kind, we approximate u(x, t) with m = 3: From [26], we have Substituting (4.4), (4.5), and (4.6) into (4.1), we get For suitable collocation points t p , we use the roots of shifted Chebyshev polynomials of the From (4.7) and the Chebyshev collocation method, we have By the property of shifted Chebyshev polynomials of the fourth kind with W * 0 (t p ) = 1, we obtain a nonlinear system of differential equations: where Using the finite difference method to solve the systems of ordinary differentials (4.10)-(4.12), we use the notations T = T final , 0 < t j ≤ T, t j = j t t = T N , j = 1, 2, . . . , N . System (4.10)-(4.12) is discretized and takes the following form: Rearranging (4.13)-(4.16) and applying the Newton-Raphson method, we arrive at the required solution for the problem on (4.1). Graphs and Table: In Fig. 1, the numerical solutions attained for Example 4.1 using the shifted Chebyshev collocation method are presented and are compared with the previous method [7] and with the exact solution as well. And it is easy to observe, the proposed method has an excellent agreement with the exact solution. Table 1 presents the numerical results of Example 4.1 at various levels of value of α. The results obtained by our method with α = 0.2, α = 0.4, α = 0.6, and α = 0.8 are given in the 3rd, 5th, 7th, and 9th columns, respectively. To validate the accuracy of the proposed method, the numerical solution is compared with the previous findings provided  in Habenom et al. [7] and the exact solutions, which shows that the present method is in good agreement with the exact solution. The findings in Table 1 demonstrate that the accuracy of the method presented is better than that of other methods. In comparison, the computing expense (CPU time) of the suggested approach tends to be lower than in other approaches, since we require just a few terms to find an answer with high accuracy. We provide a comparative study of two methods (TSM = Taylor series method (see [7]); CPFK = Chebyshev polynomial of the fourth kind (present method)) in the numerical form as follows. To use the shifted fourth-type Chebyshev polynomials, we approximate u(x, t) with m = 5 on (2.16) as follows: From the above theorem, we have For suitable collocation points t p , we use the roots of shifted Chebyshev polynomials of the fourth kind W * m+1α (t), p = 0, 1, 2, . . . , mα . For fractional α with m = 5, we have Of the property of the fourth-kind shifted Chebyshev polynomials, W * 0 (t p ) = 1. Equation (4.26) will lead us to a nonlinear system of ODE given below.
From the initial condition (IC) given at (4.18), we have where

Conclusion
In this paper, to solve a Fokker-Planck time-fractional equation, the Chebyshev collocation method and the finite difference method are applied. Using the properties of the shifted Chebyshev fourth-kind polynomials, we reduce the time-fractional Fokker-Planck equation to the system of differential equations. Also, using FDM we reduce such systems of DEs into nonlinear equations where the Newton-Raphson method can be used to solve them. As it can be observed from the tabulated values and graphical solutions in Sect. 3, the behavior of numerical solution for different values of α is presented, which shows that the proposed method is dominant over the previous study in [7] taking the exact solutions as reference values. From this we conclude that the shifted Chebyshev polynomial of the fourth kind gives us better approximation, and we recommend other scholars in this area to increase the accuracy of such methods even more than this by taking a large value of the series (2.16). The advantage of the methods is that the approximation requires only few terms of fourth-kind Chebyshev polynomials. Also, we compared the present method with the method in [7] to show that the present method demonstrated the effectiveness and high accuracy. All the numerical results were obtained with the help of MATLAB 2018a.