In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solution...In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solutions are only valid for small values of the independent variable. The DTM solutions diverge for some differential equations that extremely have nonlinear behaviors or have boundary-conditions at infinity. For this reason the governing boundary-layer equations are solved by the Multi-step Differential Transform Method (MDTM). The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, or perturbation. It is a semi analytical-numerical technique that formulizes Taylor series in a very different manner. By applying the MDTM the interval of convergence for the series solution is increased. The MDTM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions for systems of differential equations. It is predicted that the MDTM can be applied to a wide range of engineering applications.展开更多
This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)metho...This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.展开更多
This paper aims to investigate a new efficient method for solving time fractional partial differential equations.In this orientation,a reliable formable transform decomposition method has been designed and developed,w...This paper aims to investigate a new efficient method for solving time fractional partial differential equations.In this orientation,a reliable formable transform decomposition method has been designed and developed,which is a novel combination of the formable integral transform and the decomposition method.Basically,certain accurate solutions for time-fractional partial differential equations have been presented.Themethod under concern demandsmore simple calculations and fewer efforts compared to the existingmethods.Besides,the posed formable transformdecompositionmethod has been utilized to yield a series solution for given fractional partial differential equations.Moreover,several interesting formulas relevant to the formable integral transform are applied to fractional operators which are performed as an excellent application to the existing theory.Furthermore,the formable transform decomposition method has been employed for finding a series solution to a time-fractional Klein-Gordon equation.Over and above,some numerical simulations are also provided to ensure reliability and accuracy of the new approach.展开更多
Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to sol...Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.展开更多
In this paper,a reliable algorithm based on an adaptation of the standard differential transform method is presented,which is the multi-step differential transform method(MSDTM).The solutions of non-linear oscillators...In this paper,a reliable algorithm based on an adaptation of the standard differential transform method is presented,which is the multi-step differential transform method(MSDTM).The solutions of non-linear oscillators were obtained by MSDTM.Figurative comparisons between the MSDTM and the classical fourthorder Runge-Kutta method(RK4)reveal that the proposed technique is a promising tool to solve non-linear oscillators.展开更多
This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence ...This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations (PDEs) are transformed into a system of coupled non-linear ordinary differential equations (ODEs) with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method (DTM). The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.展开更多
In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.T...In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.The reduced differential transforms method(RDTM)is one of the interesting methods for finding the approximate solutions for nonlinear problems.We apply the RDTM to discuss the analytic approximate solutions to the SI1I2R model for the spread of virus HCV-subtype and SIR childhood disease model.We discuss the numerical results at some special values of parameters in the approximate solutions.We use the computer software package such as Mathematical to find more iteration when calculating the approximate solutions.Graphical results and discussed quantitatively are presented to illustrate behavior of the obtained approximate solutions.展开更多
In this paper, the genera]ised two-dimensiona] differentia] transform method (DTM) of solving the time-fractiona] coupled KdV equations is proposed. The fractional derivative is described in the Caputo sense. The pr...In this paper, the genera]ised two-dimensiona] differentia] transform method (DTM) of solving the time-fractiona] coupled KdV equations is proposed. The fractional derivative is described in the Caputo sense. The presented method is a numerical method based on the generalised Taylor series expansion which constructs an analytical solution in the form of a polynomial. An illustrative example shows that the genera]ised two-dimensional DTM is effective for the coupled equations.展开更多
In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomi...In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].展开更多
In this work, approximate analytical solutions to the lid-driven square cavity flow problem, which satisfied two-dimensional unsteady incompressible Navier-Stokes equations, are presented using the kinetically reduced...In this work, approximate analytical solutions to the lid-driven square cavity flow problem, which satisfied two-dimensional unsteady incompressible Navier-Stokes equations, are presented using the kinetically reduced local Navier-Stokes equations. Reduced differential transform method and perturbation-iteration algorithm are applied to solve this problem. The convergence analysis was discussed for both methods. The numerical results of both methods are given at some Reynolds numbers and low Mach numbers, and compared with results of earlier studies in the review of the literatures. These two methods are easy and fast to implement, and the results are close to each other and other numerical results, so it can be said that these methods are useful in finding approximate analytical solutions to the unsteady incompressible flow problems at low Mach numbers.展开更多
In this paper a new method for solving Goursat problem is introduced using Reduced Differential Transform Method (RDTM). The approximate analytical solution of the problem is calculated in the form of series with easi...In this paper a new method for solving Goursat problem is introduced using Reduced Differential Transform Method (RDTM). The approximate analytical solution of the problem is calculated in the form of series with easily computable components. The comparison of the methodology presented in this paper with some other well known techniques demonstrates the effectiveness and power of the newly proposed methodology.展开更多
In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Trans...In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.展开更多
One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transfor...One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transform method (ETM)) for computing delay differential equations (DDEs). Here, a reconstructed Elzaki transform method (RETM) is proposed for the solution of DDEs where Mamadu-Njoseh polynomials are applied as basis functions in the approximation of the analytic solution. Using this strategy, a numerical illustration as in Ref.[1] is provided to the RETM as a basis for comparison to guarantee accuracy and consistency of the method. All numerical computations were performed with MAPLE 18 software.展开更多
In this paper, the Combined Laplace Transform-Adomian Decomposition Method is used to solve nth-order integro-differential equations. The results show that the method is very simple and effective.
An approximate solution to Richards' equation is presented, mathematically describing a sort of unsaturated single phase fluid flow in porous media. The approach is a differential transform method (DTM) with interm...An approximate solution to Richards' equation is presented, mathematically describing a sort of unsaturated single phase fluid flow in porous media. The approach is a differential transform method (DTM) with intermediate variables. Two examples are given to demonstrate the accuracy of the presented solution.展开更多
An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of...An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of the shock wave. It accelerates the convergence of the series solution and provides an exact Dower series solution.展开更多
This paper concentrates on the differential transform method (DTM) to solve some delay differential equations (DDEs). Based on the method of steps for DDEs and using the computer algebra system Mathematica, we success...This paper concentrates on the differential transform method (DTM) to solve some delay differential equations (DDEs). Based on the method of steps for DDEs and using the computer algebra system Mathematica, we successfully apply DTM to find the analytic solution to some DDEs, including a neural delay differential equation. The results confirm the feasibility and efficiency of DTM.展开更多
In this paper, Differential Transform Method (DTM) is proposed for the closed form solution of linear and non-linear stiff systems. First, we apply DTM to find the series solution which can be easily converted into ex...In this paper, Differential Transform Method (DTM) is proposed for the closed form solution of linear and non-linear stiff systems. First, we apply DTM to find the series solution which can be easily converted into exact solution. The method is described and illustrated with different examples and figures are plotted accordingly. The obtained result confirm that DTM is very easy, effective and convenient.展开更多
In this paper, the Adomian Decomposition Method (ADM) and the Differential Transform Method (DTM) are applied to solve the multi-pantograph delay equations. The sufficient conditions are given to assure the convergenc...In this paper, the Adomian Decomposition Method (ADM) and the Differential Transform Method (DTM) are applied to solve the multi-pantograph delay equations. The sufficient conditions are given to assure the convergence of these methods. Several examples are presented to demonstrate the efficiency and reliability of the ADM and the DTM;numerical results are discussed, compared with exact solution. The results of the ADM and the DTM show its better performance than others. These methods give the desired accurate results only in a few terms and in a series form of the solution. The approach is simple and effective. These methods are used to solve many linear and nonlinear problems and reduce the size of computational work.展开更多
In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential e...In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential equations with initial boundary value (IBVP). According to boundary condition, the initial condition is expanded into a Fourier series. After that, the IBVP is transformed to an iterative relation in K-domain. The series solution or exact solution can be obtained. The rationality and practicability of the algorithm FDTM are verified by comparisons of the results obtained by FDTM and the existing analytical solutions.展开更多
文摘In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solutions are only valid for small values of the independent variable. The DTM solutions diverge for some differential equations that extremely have nonlinear behaviors or have boundary-conditions at infinity. For this reason the governing boundary-layer equations are solved by the Multi-step Differential Transform Method (MDTM). The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, or perturbation. It is a semi analytical-numerical technique that formulizes Taylor series in a very different manner. By applying the MDTM the interval of convergence for the series solution is increased. The MDTM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions for systems of differential equations. It is predicted that the MDTM can be applied to a wide range of engineering applications.
基金the National Natural Science Foundation of China(Grant Nos.71961022,11902163,12265020,and 12262024)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(Grant Nos.2019BS01011 and 2022MS01003)+5 种基金2022 Inner Mongolia Autonomous Region Grassland Talents Project-Young Innovative and Entrepreneurial Talents(Mingjing Du)2022 Talent Development Foundation of Inner Mongolia Autonomous Region of China(Ming-Jing Du)the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region Program(Grant No.NJYT-20-B18)the Key Project of High-quality Economic Development Research Base of Yellow River Basin in 2022(Grant No.21HZD03)2022 Inner Mongolia Autonomous Region International Science and Technology Cooperation High-end Foreign Experts Introduction Project(Ge Kai)MOE(Ministry of Education in China)Humanities and Social Sciences Foundation(Grants No.20YJC860005).
文摘This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.
基金funded by the Deanship of Research in Zarqa University,Jordan。
文摘This paper aims to investigate a new efficient method for solving time fractional partial differential equations.In this orientation,a reliable formable transform decomposition method has been designed and developed,which is a novel combination of the formable integral transform and the decomposition method.Basically,certain accurate solutions for time-fractional partial differential equations have been presented.Themethod under concern demandsmore simple calculations and fewer efforts compared to the existingmethods.Besides,the posed formable transformdecompositionmethod has been utilized to yield a series solution for given fractional partial differential equations.Moreover,several interesting formulas relevant to the formable integral transform are applied to fractional operators which are performed as an excellent application to the existing theory.Furthermore,the formable transform decomposition method has been employed for finding a series solution to a time-fractional Klein-Gordon equation.Over and above,some numerical simulations are also provided to ensure reliability and accuracy of the new approach.
文摘Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.
文摘In this paper,a reliable algorithm based on an adaptation of the standard differential transform method is presented,which is the multi-step differential transform method(MSDTM).The solutions of non-linear oscillators were obtained by MSDTM.Figurative comparisons between the MSDTM and the classical fourthorder Runge-Kutta method(RK4)reveal that the proposed technique is a promising tool to solve non-linear oscillators.
文摘This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations (PDEs) are transformed into a system of coupled non-linear ordinary differential equations (ODEs) with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method (DTM). The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.
文摘In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.The reduced differential transforms method(RDTM)is one of the interesting methods for finding the approximate solutions for nonlinear problems.We apply the RDTM to discuss the analytic approximate solutions to the SI1I2R model for the spread of virus HCV-subtype and SIR childhood disease model.We discuss the numerical results at some special values of parameters in the approximate solutions.We use the computer software package such as Mathematical to find more iteration when calculating the approximate solutions.Graphical results and discussed quantitatively are presented to illustrate behavior of the obtained approximate solutions.
基金Project supported by the Natural Science Foundation of Inner Mongolia of China (Grant No. 20080404MS0104)the Young Scientists Fund of Inner Mongolia University of China (Grant No. ND0811)
文摘In this paper, the genera]ised two-dimensiona] differentia] transform method (DTM) of solving the time-fractiona] coupled KdV equations is proposed. The fractional derivative is described in the Caputo sense. The presented method is a numerical method based on the generalised Taylor series expansion which constructs an analytical solution in the form of a polynomial. An illustrative example shows that the genera]ised two-dimensional DTM is effective for the coupled equations.
文摘In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].
文摘In this work, approximate analytical solutions to the lid-driven square cavity flow problem, which satisfied two-dimensional unsteady incompressible Navier-Stokes equations, are presented using the kinetically reduced local Navier-Stokes equations. Reduced differential transform method and perturbation-iteration algorithm are applied to solve this problem. The convergence analysis was discussed for both methods. The numerical results of both methods are given at some Reynolds numbers and low Mach numbers, and compared with results of earlier studies in the review of the literatures. These two methods are easy and fast to implement, and the results are close to each other and other numerical results, so it can be said that these methods are useful in finding approximate analytical solutions to the unsteady incompressible flow problems at low Mach numbers.
文摘In this paper a new method for solving Goursat problem is introduced using Reduced Differential Transform Method (RDTM). The approximate analytical solution of the problem is calculated in the form of series with easily computable components. The comparison of the methodology presented in this paper with some other well known techniques demonstrates the effectiveness and power of the newly proposed methodology.
文摘In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.
文摘One of the solution techniques used for ordinary differential equations, partial and integral equations is the Elzaki Transform. This paper is an extension of Mamadu and Njoseh [1] numerical procedure (Elzaki transform method (ETM)) for computing delay differential equations (DDEs). Here, a reconstructed Elzaki transform method (RETM) is proposed for the solution of DDEs where Mamadu-Njoseh polynomials are applied as basis functions in the approximation of the analytic solution. Using this strategy, a numerical illustration as in Ref.[1] is provided to the RETM as a basis for comparison to guarantee accuracy and consistency of the method. All numerical computations were performed with MAPLE 18 software.
文摘In this paper, the Combined Laplace Transform-Adomian Decomposition Method is used to solve nth-order integro-differential equations. The results show that the method is very simple and effective.
基金Project supported by the National Basic Research Program of China(973 Program)(No.2011CB013800)
文摘An approximate solution to Richards' equation is presented, mathematically describing a sort of unsaturated single phase fluid flow in porous media. The approach is a differential transform method (DTM) with intermediate variables. Two examples are given to demonstrate the accuracy of the presented solution.
基金Project supported by the National Natural Science Foundation of China(Nos.50909017,51109031, 50921001,11072053,and 51009022)the Doctoral Foundation of Ministry of Education of China(No.20100041120037)+1 种基金the Fundamental Research Funds for the Central Universities (Nos.DUT12LK52 and DUT12LK34)the Major State Basic Research Development Program of China(973 Program)(Nos.2010CB832704 and 2013CB036101)
文摘An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of the shock wave. It accelerates the convergence of the series solution and provides an exact Dower series solution.
文摘This paper concentrates on the differential transform method (DTM) to solve some delay differential equations (DDEs). Based on the method of steps for DDEs and using the computer algebra system Mathematica, we successfully apply DTM to find the analytic solution to some DDEs, including a neural delay differential equation. The results confirm the feasibility and efficiency of DTM.
文摘In this paper, Differential Transform Method (DTM) is proposed for the closed form solution of linear and non-linear stiff systems. First, we apply DTM to find the series solution which can be easily converted into exact solution. The method is described and illustrated with different examples and figures are plotted accordingly. The obtained result confirm that DTM is very easy, effective and convenient.
文摘In this paper, the Adomian Decomposition Method (ADM) and the Differential Transform Method (DTM) are applied to solve the multi-pantograph delay equations. The sufficient conditions are given to assure the convergence of these methods. Several examples are presented to demonstrate the efficiency and reliability of the ADM and the DTM;numerical results are discussed, compared with exact solution. The results of the ADM and the DTM show its better performance than others. These methods give the desired accurate results only in a few terms and in a series form of the solution. The approach is simple and effective. These methods are used to solve many linear and nonlinear problems and reduce the size of computational work.
文摘In this paper, a new Fourier-differential transform method (FDTM) based on differential transformation method (DTM) is proposed. The method can effectively and quickly solve linear and nonlinear partial differential equations with initial boundary value (IBVP). According to boundary condition, the initial condition is expanded into a Fourier series. After that, the IBVP is transformed to an iterative relation in K-domain. The series solution or exact solution can be obtained. The rationality and practicability of the algorithm FDTM are verified by comparisons of the results obtained by FDTM and the existing analytical solutions.