# International Journal of Optic and Photonic Engineering

(ISSN: 2631-5092)

Volume 5, Issue 2

Research Article

#### DOI: 10.35840/2631-5092/4526

# Applications of OHAM and MOHAM for Time-Fractional Klein-Fock-Gordon Equation

Jafar Biazar^{1*} and Saghi Safaei^{2}

### Table of Content

### Figures

**Figure 1:** The results of OHAM, MOHAM and Exact....

The results of OHAM, MOHAM and Exact solution for *α* = 2, at *t* = 0.6.

**Figure 2:** The Absolut Error of OHAM and.....

The Absolut Error of OHAM and MOHAM for *α* = 2, at *t* = 0.6.

**Figure 3:** The PLOT 3D of a) OHAM; b) MOHAM....

The PLOT 3D of a) OHAM; b) MOHAM and c) Exact solution for *α* = 2, at *t* = 0.6.

**Figure 4:** The results of OHAM and MOHAM....

The results of OHAM and MOHAM for different values of *α* = 2, at *t* = 0.6.

### Tables

** Table 1:** The value of control parameters *c _{i}*, for different values of

*α*.

** Table 2:** Values of control parameters *c _{ij}*.

** Table 3:** The results of OHAM, MOHAM and Exact solution for different value of *x* at *t* = 0.6.

### References

- I Podlubny (1999) Fractional Differential Equations. Academic Press, New York.
- MA Abdelkawy, Antonio M Lopes, Mohhamad M Babtin (2020) Shifted fractional Jacobi collocation method for solving fractional functional differential equations of variable order. Chaos, Solitons & Fractals 134: 109721.
- M Jleli, S Kumar, R Kumar, B Samet (2019) Analytical approach for time fractional wave equations in the sense of Yang-Abdel-Aty-Cattani via the homotopy perturbation transform method. Alexandria Engineering Journal 59: 2859-2863.
- KM Saad, Eman HF, Shareef AL, Alomari AK, Dumitru Baleanu (2020) On exact solutions for time-fractional Korteweg-de Vries and Korteweg-de Vries-Burger's equations using homotopy analysis transform method. Chinese Journal of Physics 63: 149-162.
- Mehmet Giyas Sakar, Fevzi Erdogan (2013) The homotopy analysis method for solving the time-fractional Fornberg-Whitham equation and comparison with Adomian's decomposition method. Applied Mathematical Modelling 37: 8876-8885.
- Mehmet Giyas Sakar, Hilmi Ergo¨ren (2015) Alternative variational iteration method for solving the time-fractional Fornberg-Whitham equation. Applied Mathematical Modelling 39: 3972-3979.
- Zaid Odibat (2019) On the optimal selection of the linear operator and the initial approximation in the application of the homotopy analysis method to nonlinear fractional differential equations. Applied Numerical Mathematics 137: 203-212.
- V Marinca, N Herisanu, C Bota, B Marinca (2009) An optimal homotopy asymptotic method applied to the steady flow of a fourth-grade fluid past a porous plate. Applied Mathematics Letters 22: 245-251.
- V Marinca, N Herisanu, I Nemes (2008) Optimal homotopy asymptotic method with application to thin film flow. Central European Journal of Physics 6: 648-653.
- V Marinca, N Herisanu (2008) Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer. International Communications in Heat and Mass Transfer 35: 710-715.
- N Herisanu, V Marinca, T Dordea, G Madescu (2008) A new analytical approach to nonlinear vibration of an electrical machine. Proceedings of the Romanian Academy, Series A 9: 229-236.
- N Herisanu, V Marinca (2010) Explicit analytical approximation to large-amplitude non-linear oscillations of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia. Meccanica 45: 847-855.
- S Iqbal, M Idrees, AM Siddiqui, A Ansari (2010) Some solutions of the linear and nonlinear Klein-Gordon equations using the optimal homotopy asymptotic method. Applied Mathematics and Computation 216: 2898-2909.
- S Iqbal, A Javed (2011) Application of optimal homotopy asymptotic method for the analytic solution of singular Lane-Emden type equation. Applied Mathematics and Computation 217: 7753-7761.
- MS Hashmi, N Khan, S Iqbal (2012) Numerical solutions of weakly singular Volterra integral equations using the optimal homotopy asymptotic method. Computers & Mathematics with Applications 64: 1567-1574.
- A Zeb, S Iqbal, AM Siddiqui, T Haroon (2013) Application of the optimal homotopy asymptotic method to flow with heat transfer of a pseudoplastic fluid inside a circular pipe. Journal of the Chinese Institute of Engineers 36: 797-805.
- S Iqbal, AR Ansari, AM Siddiqui, A Javed (2011) Use of optimal homotopy asymptotic method and galerkin's finite element formulation in the study of heat transfer flow of a third-grade fluid between parallel plates. Journal of Heat Transfer 133: 091702.
- NR Anakira, AK Alomari, AF Jameel, I Hashim (2016) Multistage Optimal homotopy asymptotic method for solving initial-value problems. J Nonlinear Sci Appl 9: 1826-1843.
- K Diethelm (2010) The analysis of fractional differential equations. Lecture Notes in Mathematics, Springer-Verlag, Berlin Heidelberg.

**Author Details**

Jafar Biazar^{1*} and Saghi Safaei^{2}

^{1}Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, Rasht, Iran

^{2}Department of Applied Mathematics, University of Guilan, Campus 2, Rasht, Iran

**Corresponding author**

Jafar Biazar, Department of Applied Mathematics, Faculty of Mathematical sciences University of Guilan, P.O. Box. 41635-19141, P.C.41938336997, Rasht, Iran.

Accepted: November 16, 2020 | Published Online: November 18, 2020

Citation: Biazar J, Safaei S (2020) Applications of OHAM and MOHAM for Time-Fractional Klein-Fock-Gordon Equation. Int J Opt Photonic Eng 5:026.

Copyright: © 2020 Biazar J, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

## Abstract

In this paper the optimal homotopy asymptotic method (OHAM) and multistage optimal homotopy asymptotic method (MOHAM) are applied to obtain an analytic approximate solution to a time-fractional Klein-Fock-Gordon (FKFG) equation. The FKFG equation plays an important role in characterizing the relativistic electrons. The MOHAM relies on OHAM to obtain analytic approximate solutions, it actually applies OHAM in each subinterval and we show that it achieves better results than OHAM over the large intervals; this is one of the advantages of this method which can be used for large intervals and to obtain good results. The convergence of the method is also addressed.

## Keywords

Optimal homotopy asymptotic method, Multistage Optimal homotopy asymptotic method, Convergence, Time-fractional Klein-Fock-Gordon (FKFG) equation

## Introduction

Fractional calculus (FC) is a part of mathematical analysis that studies the derivation of integrals and derivatives of rational orders [1]. The concept of fractional calculus (integrals and derivatives of any rational order) is established over 300 years ago, and nowadays is a very important subject. Gradually, researchers in different fields of sciences have discovered that fractional differential models have much better descriptors for different phenomena. Fractional calculus has widespread applications of in physics, chemistry, economics, dynamic systems, medical engineering, biological sciences, imaging, etc. On the other hand, physicists Klein, Fock, and Gordon have developed an equation called Klein-Fock-Gordon, which describes the relativity of electrons. This equation is a kind of the wave equation, also relates to the Schrödinger equation, and is a quantitative version of the relation of relative energy and motion. This equation is theoretically similar to the Dirac equation.

Consider the time-fractional Klein-Fock-Gordon equation as follows

$${D}^{\alpha}u\left(x,t\right)\text{}=\text{}{u}_{xx}\left(x,t\right)+au\left(x,t\right)+b{u}^{n}\left(x,t\right),\text{}1\alpha \le 2,\text{(1)}$$

Subject to the following initial conditions

$$u(x,0)\text{}=\text{}f(x)\text{}and\text{}{u}_{x}(x,0)\text{}=\text{}g(x),\text{(2)}$$

Where *a* and *b* are real constants and *n* is a positive integer.

The initial guess to the solution is as follows

$${u}_{0}(x,t)\text{}=\text{}f(x)+g(x)t.$$

Various methods have been developed to obtain approximate solutions of the fractional time differential equations such as, fractional Jacobi collocation method [2], homotopy perturbation transform method [3], homotopy analysis transform method [4], Adomian decomposition method [5], variational iteration [6], and homotopy analysis method [7]. Recently, Optimal homotopy analysis method (OHAM) was proposed by Marinca, et al. [8-12], which was used to obtain analytic approximate solutions for some nonlinear problems [13-17].

OHAM results in to satisfactory solutions on short domains, but when the interval becomes longer, the accuracy of the method decreases, so a new approach was proposed by Anakira, et al. which is called multistage optimal homotopy asymptotic method (MOHAM) that suitable for analytic approximate solutions for large intervals [18].

Finally, the approximate solutions obtained from both methods are compared with the exact solution.

## Basic Definitions of Fractional Calculus

In this section, some basic definitions of fractional calculus are explained briefly [1,19].

**Definition 4.1:** A real-valued function $f(t)$ with $t>0$ can be defined on the space ${C}_{\mu}$, $\mu \in \mathbb{R}$, if there is a real number such $\rho >\mu $

**Definition 4.2:** The Riemann-Liouville's integral of fractional order $\alpha >0$ for a continuous function $f\in {C}_{\mu}$ with $\mu >-1$ is defined as follows

${J}^{\alpha}f(t)\text{}=\text{}\frac{1}{\Gamma (\alpha )}{\displaystyle \underset{0}{\overset{t}{\int}}{(t-s)}^{\alpha -1}f(s)ds,\text{}\alpha 0,}\text{(3)}$

${J}^{0}f(t)\text{}=\text{}f(t).\text{(4)}$

By considering $f\in {C}_{\mu}$, $\mu >-1$, $\alpha ,\beta >0$ and $\gamma >-1$, the main properties of the operator ${J}^{\alpha}$ are listed as follows

$(1)\text{}{J}^{\alpha}{J}^{\beta}f(t)\text{}=\text{}{J}^{\alpha +\beta}f(t),\text{(5)}$

$(2)\text{}{J}^{\alpha}{J}^{\beta}f(t)\text{}=\text{}{J}^{\beta}{J}^{\alpha}f(t),\text{(6)}$

$$(3)\text{}{J}^{\alpha}{t}^{\gamma}\text{}=\text{}\frac{\Gamma (\gamma +1)}{\Gamma (\alpha +\gamma +1)}{t}^{\alpha +\gamma}.\text{(7)}$$

**Definition 4.3:** The fractional-order derivative of $f(t),$ in Caputo sense, is defined as follows

${D}^{\alpha}f(t)\text{}=\text{}\frac{1}{\Gamma (m-\alpha )}{\displaystyle \underset{0}{\overset{t}{\int}}{(t-s)}^{m-\alpha -1}{f}^{m}(s)ds,}\text{(8)}$

$for\text{}m-1\alpha \le m,\text{}m\in \mathbb{N},\text{}t0,\text{}f\in {C}_{-1}^{m}.$

**Lemma 4.1:** By assuming $m-1<\alpha \le m,\text{}m\in \mathbb{N},\text{}f\in {C}_{-1}^{m},$ and $\mu \ge -1$, the following properties will be valid.

$(1)\text{}{D}^{\alpha}{J}^{\alpha}f(t)\text{}=\text{}f(t),\text{(9)}$

$(2)\text{}({J}^{\alpha}{D}^{\alpha})f(t)\text{}=\text{}f(t)-{\displaystyle \sum _{k=0}^{m-1}{f}^{(k)}({0}^{+})\frac{{t}^{k}}{k!}}\text{(10)}$

## Basic Principles of the Proposed Techniques

A short introduction to the methods that will be used in this research.

### OHAM

Let's consider the following fractional equation with the boundary conditions

$L(u(\xi ))+N(u(\xi ))+f(\xi )\text{}=\text{}0,\text{}B(u(\xi ),\frac{du(\xi )}{d\xi})\text{}=\text{}0,\text{(11)}$

Where $\xi $ is an independent variable, $L,\text{}f,\text{}N$ , and $u(\xi )$ represent, respectively a linear operator, a nonlinear operator, a known function, and an unknown function. Also, $B(u(\xi ))$ shows a boundary operator.

An optimal convex homotopy ${\rm H}(\phi (\xi ;p,{c}_{i})):\Omega \times [0,1]\to \mathbb{R}$ is constructed as follows,

$$\begin{array}{l}{\rm H}(\phi (\xi ;p,{c}_{i}))\text{}=\text{}(1-p)[L(\phi (\xi ;p,{c}_{i}))+f(\xi )]\text{}=\text{(12)}\\ \text{}H(\xi ;p,{c}_{i})[L(\phi (\xi ;p,{c}_{i}))+f(\xi )+N(\phi (\xi ;p,{c}_{i}))],\text{}B(\phi (\xi ;p,{c}_{i}))\text{}=\text{}0.\end{array}$$

A deformation equation of zero-order as the following

$(1-p)[L(\phi (\xi ,p))+f(\xi )]\text{}=\text{}H(p)[L(\phi (\xi ,p))+f(\xi )+N(\phi (\xi ,p))],\text{}B(\phi (\xi ,p),\frac{d\phi (\xi ,p)}{d\xi})\text{}=\text{}0,\text{(13)}$

Where *p* is an embedding parameter in the interval $[0,1],\text{}H(p,\xi )$, is an auxiliary function with non-zero and zero outputs for $p\ne 0$ and $p\text{}=\text{}0$, respectively and ${u}_{0}(\xi )$ represents the initial condition of $u(\xi )$, and $\phi (\xi ,p)$ is an unknown function. By inserting $p\text{}=\text{}0$ and 1 into Eq. (13), the following functions are obtained

$\phi (\xi ,0)\text{}=\text{}{u}_{0}(\xi ),\text{(14)}$

$$\phi (\xi ,1)\text{}=\text{}u(\xi ).\text{(15)}$$

Therefore, $\phi (\xi ,p)$ will change continuously from the initial guess to the solution, ${u}_{0}(\xi )$ , to $u(\xi )$ when *p* increase from 0 up to 1.

By putting $p\text{=}0$ into Eq. (13), the initial solution ${u}_{0}(\xi )$ is determined as a solution for the problem

$$L({u}_{0}(\xi ))+f(\xi )\text{}=\text{}0,\text{}B({u}_{0},\frac{d{u}_{0}}{d\xi})\text{}=\text{}0.\text{(16)}$$

Next, choose an auxiliary function $H(p)$ in the following form

$$H(p)\text{}=\text{}p\text{}{c}_{1}+{p}^{2}{c}_{2}+{p}^{3}{c}_{3}+\mathrm{...}\text{},\text{(17)}$$

Where ${c}_{1},{c}_{2},{c}_{3},\mathrm{...}$ are called the convergence control parameters, that will be determined later. The auxiliary function $H(p)$ can be expressed in other forms see the paper by Herisanu and Marinca [8].

By expanding $\phi (\xi ;p,{c}_{i})$ in Taylor's series about *p*, the following expansion will be obtained

$$\phi (\xi ;p,{c}_{i})\text{}=\text{}{u}_{0}(\xi )+{\displaystyle \sum _{m=1}^{\infty}{u}_{m}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{m}){p}^{m}}\text{(18)}$$

By putting Eqs. (14-18) into (13), and equating the coefficients of the terms with identical powers of *p*, one will obtain the governing equation of the initial approximation ${u}_{0}(\xi )$, given by Eq.(16), and then the governing equation of the first order problem is defined as

$L({u}_{1}(\xi ))+f(\xi )\text{}=\text{}{c}_{1}{N}_{0}({u}_{0}(\xi )),\text{}B({u}_{1},\frac{d{u}_{1}}{d\xi})\text{}=\text{}0.\text{(19)}$

And the governing equation of the *m ^{th}* order is defined as

$L({u}_{m}(\xi ))-L({u}_{m-1}(\xi ))\text{}=\text{}{c}_{m}{N}_{0}({u}_{0}(\xi ))+{\displaystyle \sum _{j=1}^{m-1}{c}_{j}[L({u}_{m-j}(\xi ))+}$

${N}_{m-j}({u}_{0}(\xi ),{u}_{1}(\xi ),\mathrm{...},{u}_{m-1}(\xi ))],\text{}B({u}_{m},\frac{d{u}_{m}}{d\xi})\text{}=\text{}0,\text{}m\text{}=\text{}2,3,\mathrm{...}\text{(20)}$

Where ${N}_{m}({u}_{0}(\xi ),{u}_{1}(\xi ),\mathrm{...},{u}_{m}(\xi ))$ is the coefficient of ${p}^{m}$ in the expansion of $N(\phi (\xi ,p))$ about the embedding parameter *p*

$N(\phi (\xi ,p))\text{}=\text{}{N}_{0}({u}_{0}(\xi ))+{\displaystyle \sum _{m=1}^{\infty}{N}_{m}({u}_{0}(\xi ),{u}_{1}(\xi ),\mathrm{...},{u}_{m}(\xi )){p}^{m}.}\text{(21)}$

Solving Eq. (20) gives various approximates solutions ${u}_{m}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{m})$, but there exist still *m* unknowns, auxiliary parameters, ${c}_{1},{c}_{2},{c}_{3},\mathrm{..},{c}_{m}$ in the obtained solutions. It is assumed that the auxiliary parameters $({c}_{1},{c}_{2},{c}_{3},\mathrm{..},{c}_{m})$, the linear operator *L*, and the deformation equation of the zero order (16) are appropriately determined in order to ensure the convergence of the series (18) at $p\text{}=\text{}1$. Hence, putting Eqs. (14) and (15) into Eq. (18) for $p\text{}=\text{}1$ gives the solution $u(\xi )$ as

$u(\xi ,{c}_{1},{c}_{2},\mathrm{...})\text{}=\text{}{u}_{0}(\xi )+{\displaystyle \sum _{m=1}^{\infty}{u}_{m}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{m})}\text{(22)}$

The approximate solution of Eq. (11) can be calculated as

$u(\xi ,{c}_{i})\text{}=\text{}{u}_{0}(\xi )+{\displaystyle \sum _{k=1}^{m}{u}_{i}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{k}).}\text{(23)}$

Substitution of Eq. (23) in Eq. (11), results in the following residual

$R(\xi ,{c}_{i})\text{}=\text{}L(u(\xi ,{c}_{i}))+f(\xi ),\text{}i\text{}=\text{}1,2,3,\mathrm{...}\text{(24)}$

By supposing that $R(\xi ,{c}_{i})\text{}=\text{}0$, the exact solution will be $u(\xi ,{c}_{i})$. However, such a case could not be true for a nonlinear equation. By least squares technique the functional $J({c}_{i})$ should be minimized

$$J({c}_{i})\text{}=\text{}{\displaystyle \underset{a}{\overset{b}{\int}}{R}^{2}(\xi ,{c}_{i})d\xi ,}\text{(25)}$$

Where *a* and *b* are two values relating to the problem in hand. The optimal values of the unknown coefficients ${c}_{i}\text{}(i\text{}=\text{}1,2,\mathrm{...},m)$ can be determined based on the following conditions

$\frac{\partial J}{\partial {c}_{1}}\text{}=\text{}\frac{\partial J}{\partial {c}_{2}}\text{}=\text{}\frac{\partial J}{\partial {c}_{3}}\text{}=\text{}\mathrm{...}\text{}=\text{}\frac{\partial J}{\partial {c}_{m}}\text{(26)}$

In order to get an analytic approximate solution at the level *m*, the obtained optimal coefficients will be substituted in Equation (23).

### MOHAM

Although the OHAM is used to obtain approximate solutions of nonlinear problems. It has some disadvantage in nonlinear problems with large domain. To control this drawback, we introduce in this section a multistage OHAM to obtain the nonlinear problem with long of the domain. A simple way to confirm the validity of the approximate solutions of large *T* is by dividing the interval $[0,\text{}T]$ in to subinterval as $[{t}_{0},\text{}{t}_{1}),\mathrm{...},[{t}_{j-1},\text{}{t}_{j}],$ where ${t}_{j}\text{}=\text{}T$ and utilizing the MOHAM solution on each subinterval. The initial approximation in each interval is given from the solution in prior interval. First, by considering the following initial condition

$${y}_{i}({t}_{j})\text{}=\text{}{\alpha}_{i}.\text{(27)}$$

then we can obtain the initial approximation ${y}_{i,0}(t)\text{}=\text{}\alpha $ and the following zero-order equation

$(1-p)[L({\phi}_{i}(t,p))-{y}_{i,0}(t)]\text{}=\text{}{H}_{i}(p)[L({\phi}_{i}(t,p))+{f}_{i}(t)+{N}_{i}({\phi}_{i}(t,p))].\text{(28)}$

Next, we pick out the auxiliary function ${H}_{i}(p)$ in the form

$${H}_{i}(p)\text{}=\text{}{c}_{1,j}p+{c}_{2,j}{p}^{2}+{c}_{3,j}{p}^{3}+\mathrm{...},\text{(29)}$$

Or

${H}_{i}(p)\text{}=\text{}({c}_{1,j}+{c}_{2,j}+{c}_{3,j}+\mathrm{...})p.\text{(30)}$

Then, the first, second and *m ^{th}* order-approximate solutions can be generated subject to initial condition ${y}_{i,1}({t}_{j})\text{}=\text{}{y}_{i,2}({t}_{j})\text{}=\text{}\mathrm{...}\text{}=\text{}{y}_{i,m}({t}_{j})\text{}=\text{}0,$ and the approximate solution as follows,

${y}_{i}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j})\text{}=\text{}{y}_{i,0}(t)+{\displaystyle \sum _{k=1}^{m}{y}_{i,k}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{k,j}).}\text{(31)}$

Substituting Eq. (31) into Eq. (11) yields the following residual

${R}_{i}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j})\text{}=\text{}L({y}_{i}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j}))+{f}_{i}(t)+{N}_{i}({y}_{i}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j})).\text{(32)}$

If ${R}_{i}=0$, then ${y}_{i}$ will be the exact solution. Generally, such a case will not arise for nonlinear problems, but we can minimize the function

$${J}_{i}({c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j})\text{}=\text{}{\displaystyle \underset{{t}_{j}}{\overset{{t}_{j}+h}{\int}}{R}_{i}^{2}(t,{c}_{1,j},{c}_{2,j},\mathrm{...},{c}_{m,j})dt,}\text{(33)}$$

Where *h* is the length of the subinterval $\left[{t}_{j},{t}_{j+1}\right]$ and $N\text{}=\text{}[T/h]$ the number of subinterval. Now, we can solve Eq. (33) for $j\text{}=\text{}0,1,\mathrm{...},N$ with changing the initial approximation ${\alpha}_{i}$ in each subinterval from the solution in the last point of the prior interval. The unknown convergence control parameters ${c}_{i,j}(i\text{}=\text{}1,2,3,\mathrm{...},m,\text{}j\text{}=\text{}1,2,3,\mathrm{...},N)$ can be obtained by solving the following system of equations

$\frac{\partial {J}_{i}}{\partial {c}_{1,j}}\text{}=\text{}\frac{\partial {J}_{i}}{\partial {c}_{2,j}}\text{}=\text{}\mathrm{...}\text{}=\text{}\frac{\partial {J}_{i}}{\partial {c}_{m,j}}\text{}=\text{}0.\text{(34)}$

Thus, the analytic approximate solution will be as follows

$y(t)\text{}=\text{}\left\{\begin{array}{l}{y}_{1}(t),\text{}{t}_{0}\le t{t}_{1},\\ {y}_{2}(t),\text{}{t}_{1}\le t{t}_{2},\\ \vdots \\ {y}_{N}(t),\text{}{t}_{N-1}\le t\le T.\end{array}\right.\text{(35)}$

By this way, we successfully gain the solution of the initial value problem for a large interval *T*. It should be noted that if $J\text{}=\text{}0$ the MOHAM expresses the OHAM. One of the benefits of MOHAM is that it provides a simple way to control convergence and regulate convergence region and adjust the convergence region though the auxiliary function ${H}_{i}(p)$ involving several convergent control parameters ${c}_{i,j}{}^{,}$s. In general, this method eliminates the difficulty of finding approximate solutions in large ranges.

## Convergence Theorem

If the series (23) converges to $u(\xi )$, where ${u}_{n}(\xi )\in L({R}^{+})$ is produced by Eqs. (16)-(19), and the satisfies the ${n}^{th}$ order deformation equation (20), then $u(\xi )$ is the exact solution of (11) [18].

**Proof: **

Since the series $\sum _{n=1}^{\infty}{u}_{n}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{n})$ is convergent, it can be written as follows

$$S(\xi )\text{}=\text{}{\displaystyle \sum _{n=1}^{\infty}{u}_{n}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{n})}\text{(36)}$$

so, we have

$\underset{n\to \infty}{\mathrm{lim}}{u}_{n}(\xi ,{c}_{1},{c}_{2},\mathrm{...},{c}_{n})\text{}=\text{}0.\text{(37)}$

The left hand-side of (20) satisfies

${u}_{1}(\xi ,{c}_{1})+{\displaystyle \sum _{j=2}^{n}{u}_{j}(\xi ,{\overrightarrow{c}}_{j})-{\displaystyle \sum _{j=2}^{n}{u}_{j-1}(\xi ,{\overrightarrow{c}}_{j-1})\text{}=\text{}{u}_{n}(\xi ,{\overrightarrow{c}}_{n}).}}\text{(38)}$

According to (38) and the limit, we have

${u}_{1}(\xi ,{c}_{1})+{\displaystyle \sum _{j=2}^{\infty}{u}_{j}(\xi ,{\overrightarrow{c}}_{j})-}{\displaystyle \sum _{j=2}^{\infty}{u}_{j-1}(\xi ,{\overrightarrow{c}}_{j-1})\text{}=\text{}\underset{n\to \infty}{\mathrm{lim}}}\text{}{u}_{n}(\xi ,{\overrightarrow{c}}_{n})\text{}=\text{}0.\text{(39)}$

Applying the linear operator

$L({u}_{1}(\xi ,{c}_{1}))+L({\displaystyle \sum _{j=2}^{\infty}{u}_{j}(\xi ,{\overrightarrow{c}}_{j}))-L(}{\displaystyle \sum _{j=2}^{\infty}{u}_{j-1}(\xi ,{\overrightarrow{c}}_{j-1}))\text{}=\text{}0.}\text{(40)}$

Eq.(40) can be written as following

$L({u}_{1}(\xi ,{c}_{1}))+{\displaystyle \sum _{j=2}^{\infty}(L({u}_{j}(\xi ,{\overrightarrow{c}}_{j}))-}L({u}_{j-1}(\xi ,{\overrightarrow{c}}_{j-1})))\text{}=\text{}$

${c}_{1}{N}_{0}({u}_{0})+f(\xi )+{\displaystyle \sum _{j=2}^{\infty}({c}_{j}{N}_{0}({u}_{0})+{\displaystyle \sum _{k=1}^{j-1}{c}_{k}[L({u}_{j-k})+{N}_{j-k}({\overrightarrow{u}}_{j-1})]})}\text{}=\text{}$

${c}_{1}{N}_{0}({u}_{0})+f(\xi )+{\displaystyle \sum _{j=2}^{\infty}{c}_{j}{N}_{0}({u}_{0})+{\displaystyle \sum _{j=2}^{\infty}{\displaystyle \sum _{k=1}^{j-1}({c}_{k}[L({u}_{j-k})+{N}_{j-k}({\overline{u}}_{j-1})])}\text{}=}}\text{}$

$\sum _{j=2}^{\infty}{\displaystyle \sum _{k=1}^{j-1}{c}_{k}[L({u}_{j-k})]+{\displaystyle \sum _{j=2}^{\infty}{\displaystyle \sum _{k=1}^{j-1}({c}_{j}{N}_{0}({u}_{0})+{c}_{k}[{N}_{j-k}({\overline{u}}_{j-1})])+f(\xi )\text{}=\text{}0.}}}}\text{(41)$

So by choosing the optimal ${c}_{k},\text{}k\text{}=\text{}1,2,3,\mathrm{...}$, Eq.(41) is converted to the following

$L(u(\xi ))+N(u(\xi ))+f(\xi )\text{}=\text{}0.\text{(42)}$

Which is the exact solution of the problem.

## Solution of the Fractional KFG Equation

We consider one of the nonlinear cases of FKFG equation, for $a\text{}=\text{}-\frac{3}{4},b\text{}=\text{}\frac{3}{2},$ and $n\text{}=\text{}2,$ therefore Eq. (1) becomes

${D}_{t}^{\alpha}u(x,t)\text{}=\text{}{u}_{xx}(x,t)-\frac{3}{4}u(x,t)+\frac{3}{2}{u}^{3}(x,t),\text{}1\alpha \le 2.\text{(43)}$

With the initial condition

$u(x,0)\text{}=\text{}\mathrm{sech}(x)\text{}and\text{}{u}_{t}(x,0)\text{}=\text{}\frac{1}{2}\mathrm{sech}(x)\mathrm{tanh}(x).\text{(44)}$

and the exact solution

$$u(x,t)\text{}=\text{}-\mathrm{sech}(x+0.5t).$$

Now, we solve Eq. (43) by OHAM and MOHAM.

### Solution of FKFG equation by OHAM

Having the linear operator $L\text{}=\text{}{D}_{t}^{\alpha}$ and nonlinear operator

$$N\text{}=\text{}-{u}_{xx}(x,t)+\frac{3}{4}u(x,t)-\frac{3}{2}{u}^{3}(x,t)$$

and using Eq. (16), we have

${u}_{0}(x,t)\text{}=\text{}\frac{1}{2}\mathrm{sech}(x)\mathrm{tanh}(x)t+\mathrm{sech}(x),\text{(45)}$

The following problems are resulted from Eqs. (19) - (20)

*The First-order problem:*

${D}_{t}^{\alpha}({u}_{1}(x,t))\text{}=\text{}{c}_{1}[-{D}_{xx}({u}_{0}(x,t))+\frac{3}{4}{u}_{0}(x,t)-\frac{3}{2}{u}_{0}^{3}(x,t)],\text{(46)}$

*The Second-order problem:*

$${D}_{t}^{\alpha}({u}_{2}(x,t))-{D}_{t}^{\alpha}({u}_{1}(x,t))\text{}=\text{}{c}_{1}[{D}_{t}^{\alpha}{u}_{1}(x,t)-{D}_{xx}{u}_{1}(x,t)+\frac{3}{4}{u}_{1}(x,t)-$$

$\frac{3}{2}(3{u}_{0}^{2}(x,t){u}_{1}(x,t))]+{c}_{2}[-{D}_{xx}{u}_{0}(x,t)+\frac{3}{4}{u}_{0}(x,t)-\frac{3}{2}{u}_{0}^{3}(x,t)],\text{(47)}$

*The Third-order problem:*

${D}_{t}^{\alpha}{u}_{3}(x,t)-{D}_{t}^{\alpha}{u}_{2}(x,t)\text{}=\text{}{c}_{1}[{D}_{t}^{\alpha}{u}_{2}(x,t)-{D}_{xx}{u}_{2}(x,t)+\frac{3}{4}{u}_{2}(x,t)-\frac{3}{2}(3{u}_{0}^{2}(x,t){u}_{2}(x,t)+$

$3{u}_{1}^{2}(x,t){u}_{0}(x,t))]+{c}_{2}[{D}_{t}^{\alpha}{u}_{1}(x,t)-{D}_{xx}{u}_{1}(x,t)+\frac{3}{4}{u}_{1}(x,t)-\frac{3}{2}(3{u}_{0}^{2}(x,t){u}_{1}(x,t))]+$

${c}_{3}[-{D}_{xx}{u}_{0}(x,t)+\frac{3}{4}{u}_{0}(x,t)-\frac{3}{2}{u}_{0}^{3}(x,t)],\text{(48)}$

By solving the Eqs. (46) - (48), analytic approximate solutions will be obtained.

In this research, a four terms approximation for $u(x,t)$ is considered. By substituting the solutions of the zero-order, first-order, second-order and third-order into Eq. (43) and using the least squares technique, the parameters ${c}_{1},{c}_{2}\text{}and\text{}{c}_{3}$ are determined for different values of *α*, as in Table 1.

### Solution of FKFG equation by MOHAM

In this section, we utilize MOHAM to Eq. (43). We will consider the auxiliary function $H(p,t)$ as the following

$H(p,t)\text{}=\text{}({c}_{1,j}+{c}_{2,j}t+{c}_{3,j}{t}^{2})p.\text{(49)}$

This approach leads to the following sequence of equations

$${p}^{0}:\text{}{D}_{t}^{\alpha}({u}_{0j}(x,t))\text{}=\text{}{\alpha}_{j},\text{(50)}$$

${p}^{1}:\text{}{D}_{t}^{\alpha}({u}_{1j}(x,t))\text{}=\text{}{c}_{1j}[-{D}_{xx}{u}_{0j}(x,t)+\frac{3}{4}{u}_{0j}(x,t)-\frac{3}{2}{u}_{0j}^{3}(x,t)],\text{(51)}$

$${p}^{2}:\text{}{D}_{t}^{\alpha}({u}_{2j}(x,t))\text{}=\text{}-{D}_{t}^{\alpha}({u}_{1j}(x,t))\text{}=\text{}{c}_{1j}[{D}_{t}^{\alpha}{u}_{1j}(x,t)-{D}_{xx}{u}_{1j}(x,t)+\frac{3}{4}{u}_{1j}(x,t)-,$$

$\frac{3}{2}(3{u}_{0j}^{2}(x,t){u}_{1j}(x,t))]+{c}_{2j}[-{D}_{xx}{u}_{0j}(x,t)+\frac{3}{4}{u}_{0j}(x,t)-\frac{3}{2}{u}_{0j}^{3}(x,t)].\text{(52)}$

${p}^{3}:\text{}{D}_{t}^{\alpha}{u}_{3j}(x,t)-{D}_{t}^{\alpha}{u}_{2j}(x,t)\text{}=\text{}{c}_{1j}[{D}_{t}^{\alpha}{u}_{2j}(x,t)-{D}_{xx}{u}_{2j}(x,t)+\frac{3}{4}{u}_{2j}(x,t)-$

$\frac{3}{2}(3{u}_{0j}^{2}(x,t){u}_{2j}(x,t)+3{u}_{1j}^{2}(x,t){u}_{0j}(x,t))]+{c}_{2j}[{D}_{t}^{\alpha}{u}_{1j}(x,t)-{D}_{xx}{u}_{1j}(x,t)+\frac{3}{4}{u}_{1j}(x,t)-$

$\frac{3}{2}(3{u}_{0j}^{2}(x,t){u}_{1j}(x,t))]+{c}_{3j}[-{D}_{xx}{u}_{0j}(x,t)+\frac{3}{4}{u}_{0j}(x,t)-\frac{3}{2}{u}_{0j}^{3}(x,t)],\text{(53)}$

As in OHAM, we obtain two approximate solutions to Equation (43)

$${u}_{j}(x,t)\text{}=\text{}{u}_{0j}(x,t)+{u}_{1j}(x,t).\text{(54)}$$

By applying the least squares technique, the parameters ${c}_{1j},{c}_{2j}\text{}and\text{}{c}_{3j}$ are determined for different values of *α* where $h\text{}=\text{}0.2$, and ${t}_{0}\text{}=\text{}0$ up to ${t}_{2}\text{}=\text{}T\text{}=\text{}1$ as in Table 2.

Approximate solution for $\alpha \text{}=\text{}2$ is as the following form

$$u\left(x,t\right)\text{}=\text{}\left\{\begin{array}{c}\begin{array}{l}0.5000000000\mathrm{sech}\left(x\right)\mathrm{tanh}\left(x\right)t-\left(\mathrm{sech}\left(x\right)-0.007922060576\right){\mathrm{sech}}^{3}\left(x\right){\mathrm{tanh}}^{3}\left(x\right){t}^{7}+\\ \left(0.003558142171{\mathrm{sech}}^{3}\left(x\right){\mathrm{tanh}}^{3}\left(x\right)+0.06654530884{\mathrm{sech}}^{3}\left(x\right){\mathrm{tanh}}^{2}\left(x\right)\right){t}^{6}+\mathrm{...},\text{}0\le x0.2\end{array}\\ \frac{1}{{\mathrm{cosh}}^{6}\left(x\right)}\left(\begin{array}{l}\text{}\left(-0.008882058897\mathrm{sinh}\left(x\right)\mathrm{cosh}{\left(x\right)}^{2}+0.008882058897\mathrm{sinh}\left(x\right)\right){t}^{7}+\\ \left(0.07461{\mathrm{cosh}}^{3}\left(x\right)+0.0048\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)-0.0746\mathrm{cosh}\left(x\right)-0.004833\mathrm{sinh}\left(x\right)\right){t}^{6}+\mathrm{...}\end{array}\right)\text{}0.2\le x0.4\text{}\\ \frac{1}{{\mathrm{cosh}}^{6}\left(x\right)}\left(\begin{array}{l}\text{}\left(0.00546237854\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)-0.00546237854\mathrm{sinh}\left(x\right)\right){t}^{7}+\\ \left(-0.045884{\mathrm{cosh}}^{3}\left(x\right)-0.12506588\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)+0.045884\mathrm{cosh}\left(x\right)+0.0125\mathrm{sinh}\left(x\right)\right){t}^{6}+\mathrm{...}\end{array}\right)\text{}0.4\le x0.6\text{}\\ \frac{1}{{\mathrm{cosh}}^{6}\left(x\right)}\left(\begin{array}{l}\text{}\left(0.002476348954\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)-0.002476348954\mathrm{sinh}\left(x\right)\right){t}^{7}+\\ \left(-0.02080{\mathrm{cosh}}^{3}\left(x\right)-0.006727\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)+0.0208\mathrm{cosh}\left(x\right)+0.006727\mathrm{sinh}\left(x\right)\right){t}^{6}+\mathrm{...}\end{array}\right)\text{}0.6\le x0.8\text{}\\ \frac{1}{{\mathrm{cosh}}^{6}\left(x\right)}\left(\begin{array}{l}\text{}\left(0.001751944671\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)-0.001751944671\mathrm{sinh}\left(x\right)\right){t}^{7}+\\ \left(-0.014716{\mathrm{cosh}}^{3}\left(x\right)-0.006178\mathrm{sinh}\left(x\right){\mathrm{cosh}}^{2}\left(x\right)+0.0147\mathrm{cosh}\left(x\right)+0.006178\mathrm{sinh}\left(x\right)\right){t}^{6}+\mathrm{...}\end{array}\right)\text{}0.8\le x\le 1\text{}\end{array}\right.$$

From Figure 1, Figure 2, Figure 3 and Figure 4 one can see that the solutions obtained by OHAM and MOHAM are nearly identical with the exact solution Table 3.

### Conclusion

In this study OHAM and MOHAM are used to derive an analytic approximate solution for the time-fractional Klein-Fock-Gordon (FKFG) equation. The results obtained from these methods show that MOHAM converges better than OHAM. One observes that the results agree very well with the exact solution. The MOHAM by dividing the interval $[o,\text{}T]$ can obtain better solution than OHAM. As far as the authors are aware, MOHAM has not been used to solve fractional partial differential equations, so far, the method has been tested on fractional-PDE and yields to satisfactory results. The Figures and Tables expose that good results are obtained by MOHAM and more accurate solution as compared to OHAM. The convergence of the method is proved.