删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Analytical and Approximate Solutions for Complex Nonlinear Schr【-逻*辑*与-】#xF6;dinger Equation via Gen

本站小编 Free考研考试/2022-01-02

Mostafa M. A. Khater,1,*, Dian-Chen Lu,2,?, Raghda A. M. Attia1, Mustafa Inç,2,§ Department of Mathematics, Faculty of Science, Jiangsu University, China
Science Faculty, Firat University, 23119 Elazig, Turkey

Corresponding authors: * E-mail:mostafa.khater2024$@$yahoo.com? E-mail:dclu$@$ujs.edu.cn? E-mail:rony_695$@$yahoo.com§ E-mail:minc$@$firat.edu.tr

Received:2019-03-13Online:2019-11-1


Abstract
This article studies the performance of analytical, semi-analytical and numerical scheme on the complex nonlinear Schrödinger (NLS) equation. The generalized auxiliary equation method is surveyed to get the explicit wave solutions that are used to examine the semi-analytical and numerical solutions that are obtained by the Adomian decomposition method, and B-spline schemes (cubic, quantic, and septic). The complex NLS equation relates to many physical phenomena in different branches of science like a quantum state, fiber optics, and water waves. It describes the evolution of slowly varying packets of quasi-monochromatic waves, wave propagation, and the envelope of modulated wave groups, respectively. Moreover, it relates to Bose-Einstein condensates which is a state of matter of a dilute gas of bosons cooled to temperatures very close to absolute zero. Some of the obtained solutions are studied under specific conditions on the parameters to constitute and study the dynamical behavior of this model in two and three-dimensional.
Keywords: complex NLS equation;generalized auxiliary equation method;adomian decomposition method;B-splines schemes (cubic & quintic & septic)


PDF (2921KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Mostafa M. A. Khater, Dian-Chen Lu, Raghda A. M. Attia, Mustafa Inç. Analytical and Approximate Solutions for Complex Nonlinear Schrödinger Equation via Generalized Auxiliary Equation and Numerical Schemes. [J], 2019, (11): 1267-1274 doi:10.1088/0253-6102/71/11/1267

1 Introduction

Recently, the nonlinear partial differential equations have been used for modelling many complex natural phenomena such as giant permafrost explosions, desert roses, blood falls, bismuth crystals, synchronized hordes of cicadas, Salar de Uyuni, Pele's hair lava, penitentes, aka fire tornadoes, fire whirls, halos, fire rainbows, volcanic lightning, and so on. Many analytical methods have been derived to study the explicit wave solutions of these complex phenomena that help in testing the accuracy of numerical schemes. It is also considered as a motivation for discovering more accurate numerical schemes that study the approximate solutions of these phenomena. For achieving this goal, the exact analytical wave solutions methods have experienced significant fast growth of development. Many methods have been derived to investigate the analytical and numerical solutions of many phenomena modeling such as modified exp $( -\Phi(\xi)) $ expansion method, generalized auxiliary equation method, $( {G'}/{G}) $-expansion method, Khater method, modified auxiliary equation method (modified Khater method), generalized Kudryashov method, generalized Riccati expansion method, generalized tanh-function method, improved Bernoulli subequation function, improved F-expansion method, sinh-cosh method, and so on.[1-15]

In 1925, Erwin Schrödinger derived the famous nonlinear evolution equation and called it by a complex nonlinear Schrödinger (NLS) equation.[16-25] The NLS equation relates with a quantum state, fiber optics, and water waves. The NLS illustrates the evolution of slowly varying packets of quasi-monochromatic waves, wave propagation, and the envelope of modulated wave groups, respectively for these branches. In the following lines, we mention some simple steps to show how this equation derived.

In studying of a particle in an infinite potential well, the particle of fixed energy $E$ has a wave function that can be represented as a linear next combination

$$ S(x,t)=\Lambda \exp(i (\kappa x+\varpi t))\,, $$
where $\Lambda$, $(\kappa, \varpi)$, $x$, $t$, and $i $ are the amplitude of the wave, arbitrary constants, the positive direction, time, and $\sqrt{-1}$ respectively. The momentum and energy of free particle can be expressed in the next formulas

$$ p=\kappa h\,, \quad E=\varpi h\,. $$

So that, Eq. (1) can be rewritten as follows

$$ S_{x x}=-\kappa^{2} S\,. $$
Using the relation following between $\kappa$, $h$, $p$, and $m$:

$$ \frac{p^{2}}{2 m}=\frac{h^{2} \kappa^{2}}{2 m} \Rightarrow \kappa^{2}=\frac{p^{2}}{h^{2}}\,. $$

This relation leads to the next formula

$$ S_{xx}=\frac{-p^{2}}{h^{2}} S\,. $$
Similarly, we can write

$$ S_{t}=-i \varpi S \Rightarrow S_{t}=E S\,. $$
For both of a kinetic energy and a potential energy, we find

$$ E=\frac{p^{2}}{2 m}+\varphi(x) \Rightarrow E S =\frac{p^{2}}{2 m} S+\varphi(x) S\,, $$
where $S$, and $\varphi$ are wave, and potential functions, respectively. According to Eqs. (2), we have

$$ 2 m \varphi(x) S-h^{2} S_{x\,x}-2 i m h S_{t}=0\,. $$
Equation (6) is the nonlinear complex Schrödinger equation. This equation takes many various formulas, we will mention some of them in the following equations

(i) Schrödinger equation:

$$ i \psi_{t}=-\frac{\nabla^{2} \psi}{2M}+V(x,y)\psi\,, $$
where

$$ \nabla^{2}=\frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}}\,.$$

(ii) Nonlinear cubic Schrödinger equation:

$$ i h \psi_{t}=-\frac{h^{2}}{2 m} \psi_{x x}+V\psi-\gamma |\psi|^{2} \psi\,. $$
(iii) Nonlinear Schrödinger equation:

$$ i u_{t}+u_{x x}+|u|^{2}u=0\,. $$
(iv) Klein-Gordon-Schrödinger equations:

$$ i \psi_{t}+\frac{1}{2} \Delta \psi=-\phi \psi\,, \quad \phi_{t t}-\Delta \phi+ M^{2} \phi=|\psi|^{2}\,. $$
(v) Cubic-quintic Schrödinger equation:

$$ i u_{t}+u_{x\,x}=a_{1} u |u|^{2}+a_{2} u |u|^{4}\,. $$
(vi) Perturbed nonlinear Schrödinger equation:

$$ i u_{t}+u_{x x}+\alpha u |u|^{2}+i [\gamma_{1} u_{x x x} \quad +\gamma_{2} |u|^{2} u_{x}+\gamma_{3} u (|u|^{2})_{x}]=0\,. $$
(vii) Higher-order nonlinear Schrödinger equation:

$$ \psi_{z}=i\alpha_{1}\psi_{t\,t}+i\alpha_{2}\psi|\psi|^{2}+\alpha_{3} \psi_{t t t} \hphantom{\psi_{z}=}+\alpha_{4} (\psi|\psi|^{2})_{t}+\alpha_{5}(|\psi|^{2})_{t}\,. $$
Many various analytical solutions methods have been applying to these different forms of the famous NLS equation. In our paper, we select one of these forms and implement an analytical technique (generalized auxiliary equation method)[26-28] and four different numerical schemes (Adomian decomposition method, and B-spline schemes)[29-30] that never applied before to this form of the complex NLS equation to estimate the explicit wave solutions and the approximate solutions. This goal of this study is showing which one of these numerical schemes is more accurate than other.

Rest of the article is organized as follows. In Sec. 2, we apply the generalized auxiliary equation method on the complex NLS equation. In Sec. 3, we test our obtained solutions, by comparison, them with approximate solutions that obtained three numerical scheme. In Sec. 5, we evince the conclusion of our paper.

2 Analytical Solutions

The complex NLS equation has the next formula

$$ i\nu_{t}+\nu_{x x}+2\delta |\nu|^{2}\nu=0\,, $$
where $i=\sqrt{-1}$ and $\nu(x,t)$ represents the wave function. Using the next wave transformation $\nu(x,t)=e^{i\eta} \nu(\theta)$ where $\{\theta=x+ct$ and $\eta=rx+st\}$ on Eq. (14), yields

$$ \nu^{\prime\prime}-(r^{2}+s) \nu+2 \delta \nu^{3}=0\,, $$
where $c=-2 r$. According to the generalized auxiliary equation method, the general solutions of Eq. (15) are in the next form

$$ \nu(\theta)=\sum_{i=0}^{n} a_{i}\phi^{i}(\theta)\,, $$
where $\phi(\theta)$ is the solution of the next ordinary equation

$$ \phi'(\theta)=\sqrt{\sum_{i=0}^{m} h_{i}\phi^{i}(\theta)}\,. $$
Balancing the terms in Eq. (15), we get a relation between $n$ and $m$ in the following form

$$ 2 n+2=m\,. $$

Taking $n=2$, leads to $m=6$. Thus, the general solution of Eq. (15) is in the next form

$$ \nu(\theta)=a_{0}+a_{1}\phi(\theta)+a_{2}\phi^{2}(\theta)\,. $$
Substituting Eqs. (17) and (18) into Eq. (15) and collecting all coefficients of same power of $\phi^{i}(\theta)\, \{i=0,1,2,\ldots,6\}$, leads to a system of algebraic equations. Solving this system, we obtain

$$a_0= \frac{i h_4}{2 \sqrt{\delta } \sqrt{h_6}}\,, \quad a_1= 0\,,\quad a_2= \frac{2 i \sqrt{h_6}}{\sqrt{\delta }}\,, \\ s= \frac{-2 h_6 r^2-h_4^2}{2 h_6}\,,\quad h_2= \frac{h_4^2}{4 h_6}\,, $$

where $\delta<0$, $h_{6}>0$. So that, the solitary wave solutions of Eq. (14) are in the next forms.

it Family 1 When, $h_{0}=h_{1}=h_{3}=h_{5}=0$:

it Case 1 When, $h_{2}>0$

$$ \nu_{1}(x,t)=\frac{i h_4 \exp\big[{i (r (x-r t)-\frac{h_4^2 t}{2 h_6})}\big] \big(1-\frac{4 h_2 h_6 \text{sech}^2(\sqrt{h_2} (-2 r t+x+\vartheta ))}{h_4^2-h_2 h_6 (\epsilon \tanh (\sqrt{h_2} (-2 r t+x+\vartheta ))+1){}^2}\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,, $$
$$ \nu_{2}(x,t)=\frac{i h_4 \exp\big[{i (r (x-r t)-\frac{h_4^2 t}{2 h_6})}\big] \big(\frac{4 h_2 h_6 \text{csch}^2 ((\sqrt{h_2} (-2 r t+x+\vartheta ))}{h_4^2-h_2 h_6 (\epsilon \coth (\sqrt{h_2} (-2 r t+x+\vartheta ))+1){}^2}+1\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,. $$
When, $h_{2}>0$, $h_{6}>0$

$$ \nu_{3}(x,t)=\frac{i \exp\big[{i (r (x-r t)-\frac{h_4^2 t}{2 h_6})}\big] \big(h_4-\frac{4 h_2 h_6 \text{sech}^2(\sqrt{h_2} (-2 r t+x+\vartheta ))}{2 \sqrt{h_2 h_6} \epsilon \tanh (\sqrt{h_2} (-2 r t+x+\vartheta ))+h_4}\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,, $$
$$ \nu_{4}(x,t)=\frac{i \exp\big[i (r (x-r t)-\frac{h_4^2 t}{2 h_6})\big] \big(\frac{4 h_2 h_6 \text{csch}^2(\sqrt{h_2} (-2 r t+x+\vartheta ))}{2 \sqrt{h_2 h_6} \epsilon \coth (\sqrt{h_2} (-2 r t+x+\vartheta ))+h_4}+h_4\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,. $$
When, $h_{2}<0$, $h_{6}>0$


$$\nu_{5}(x,t)=\frac{i \exp\big[{i (r (x-r t)-\frac{h_4^2 t}{2 h_6}\big)} \big(h_4-\frac{4 h_2 h_6 \text{sech}^2(\sqrt{h_2} (-2 r t+x+\vartheta ))}{2 \sqrt{-h_2 h_6} \epsilon \tan (\sqrt{-h_2} (-2 r t+x+\vartheta ))+h_4}\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,, $$
$$ \nu_{6}(x,t)=\frac{i \exp\big[{i (r (x-r t)-\frac{h_4^2 t}{2 h_6}\big)} \big(\frac{4 h_2 h_6 \text{csch}^2(\sqrt{h_2} (-2 r t+x+\vartheta ))}{2 \sqrt{-h_2 h_6} \epsilon \cot (\sqrt{-h_2} (-2 r t+x+\vartheta ))+h_4}+h_4\big)}{2 \sqrt{\delta } \sqrt{h_6}}\,. $$
it Family 2 When, $h_{0}=h_{1}=h_{3}=h_{5}=0$:

it Case 1 When, $h_{4} ^{2}=4 h_{2} h_{6}$

$$ \nu_{7}(x,t)=\frac{i \exp\big[{-\frac{1}{2} i (\frac{h_4^2 t}{h_6}+2 r (r t-x))}\big] \big(h_4^2-4 h_2 h_6 (\epsilon \tanh (\sqrt{h_2} \epsilon (-2 r t+x+\vartheta ))+1)\big)}{2 \sqrt{\delta } h_4 \sqrt{h_6}}\,, $$
$$ \nu_{8}(x,t)=\frac{i \exp\big[{-\frac{1}{2} i (\frac{h_4^2 t}{h_6}+2 r (r t-x))}\big] \big(h_4^2-4 h_2 h_6 (\epsilon \coth (\sqrt{h_2} \epsilon (-2 r t+x+\vartheta ))+1)\big)}{2 \sqrt{\delta } h_4 \sqrt{h_6}}\,, $$
where $\epsilon=\pm1$.

3 Approximate Solutions

3.1 Adomian Decomposition Method

Applying the Adomian decomposition method to Eq. (25) allows obtaining the following solutions:

$$ u_{0}=-\frac{\theta}{2}\,,$$
$$u_{1} =\frac{\theta^3}{6}-\frac{\theta^5}{20}\,, $$
$$ u_{2} =\frac{\theta^{10}}{600}-\frac{\theta^8}{112}+\frac{\theta^7}{420}-\frac{\theta^5}{60}\,, $$
$$ u_{3} =-\frac{\theta^{13}}{5200}-\frac{\theta^{12}}{39600}+\frac{\theta^{11}}{550} \\ \hphantom{ u_{3} =} +\frac{\theta^{10}}{5040}-\frac{67 \theta^9}{7560}+\frac{31 \theta^7}{1260}\,. $$
According to Eqs. (27)--(30), we get the approximate solution of Eq. (25) in the next formula

$$u_{\rm approximate}=-\frac{\theta^{13}}{5200}-\frac{\theta^{12}}{39600}+\frac{\theta^{11}}{550}+\frac{47 \theta^{10}}{25200}-\frac{67 \theta^9}{7560} \\ -\frac{\theta^8}{112}+\frac{17 \theta^7}{630}-\frac{\theta^5}{15}+\frac{\theta^3}{6}-\frac{\theta}{2}+\cdots$$

Now, we give a comparison between exact and approximate solutions at different points.

3.2 Cubic Splines

Using this kind of techniques allows putting the approximate solution of Eq. (25) in the next formula

$$ u(\theta)=\sum_{i=-1} ^{n+2} c_{i}B_{i}\,, $$
where $c_{i}$, $B_{i}$ satisfy the collection condition, respectively:

$$ Lu(\theta)=f(\theta_{i},u(\theta_{i}))\,, \quad \text{where}\ \ (i=0,1,\ldots,n)\,. $$

$$ B_{i}(\theta)=\frac{1}{h^{3}}\left\{ \begin{aligned} (\theta-\theta_{i-2})^{7}\,, & \quad \text{where} \ \ \theta \in [\theta_{i-2},\theta_{i-1}]\,, \\ \vdots \\ (\theta_{i+2}-\theta)^{3}\,, & \quad \text{where} \ \ \theta\in [\theta_{i+3},\theta_{i+4}]\,, \\ 0\,, & \quad \text{otherwise} \,, \end{aligned} \right. $$
where

$$ i \in [-1,n+2]\,, $$

we obtain

$$ u_{i}(\theta)=c_{i-3}+120c_{i-2}+1191c_{i-1}+2416c_{i} \\ +1191c_{i+1}+120 c_{i+2}+c_{i+3}\,. $$
Substituting Eq. (33), its derivatives and boundary conditions of Eq. (25) into Eq. (15), gets $(n+2)$ of equations. Solving this system of equations to get the value of $c_{i}$, we obtain the following Table 2.


Table 1
Table 1Exact and approximate solutions of Eq. (15) at different value of $\theta$ to show the absolute error between both of them.

New window|CSV


Table 2
Table 2Applying the cubic spline technique to determine exact and approximate solutions of Eq. (15) at different values of $\theta$ to study the absolute error of them.

New window|CSV

3.3 Quintic Splines

Using this kind of techniques allows putting the approximate solution of Eq. (25) in the next formula

$$ u(\theta)=\sum_{i=-2} ^{n+1} c_{i}B_{i}\,, $$
where $c_{i}$, $B_{i}$ satisfy the collection condition, respectively:

$$Lu(\theta)=f(\theta_{i},u(\theta_{i})) \ \ \text{where} \ \ (i=0,1,\ldots,n)\,. $$

$$ B_{i}(\theta)=\frac{1}{h^{5}}\left\{ \begin{aligned} (\theta-\theta_{i-3})^{5}\,, & \quad \text{where} \ \ \theta \in [\theta_{i-3},\theta_{i-2}]\,, \\ \vdots \\ (\theta_{i+3}-\theta)^{5}\,, & \quad \text{where} \ \ \theta \in [\theta_{i+2},\theta_{i+3}]\,, \\ 0\,, & \quad \text{otherwise}\,. \end{aligned} \right. $$
where

$ i \in [-2,n+2]$,

we obtain

$$ u_{i}(\theta)=c_{i-2}+26c_{i-1}+66c_{i}+26c_{i+1}+c_{i+2}\,. $$
Substituting Eq. (36), its derivatives and boundary conditions of Eq. (25) into Eq. (15), leads to $(n+5)$ of equations. Solving this system of equations to get the value of $c_{i}$, we obtain the following Table 3.


Table 3
Table 3Implementation of the quintic spline technique to determine exact and approximate solutions of Eq. (15) at various values of $\theta$ to represent the absolute error of them.

New window|CSV

3.4 Septic Splines

Using this kind of techniques allows putting the approximate solution of Eq. (25) in the next formula

$$ u(\theta)=\sum_{i=-3} ^{n+3} c_{i}B_{i}\,, $$
where $c_{i}$, $B_{i}$ satisfy the collection condition, respectively:

$$L u(\theta)=f(\theta_{i},u(\theta_{i})) \quad \text{where} \ \ (i=0,1,\ldots,n)\,. $$

$$ B_{i}(\theta)=\frac{1}{h^{7}}\left\{ \begin{aligned} (\theta-\theta_{i-4})^{7}\,, & \quad \text{where} \ \ \theta\in [\theta_{i-4},\theta_{i-3}]\,, \\ \qquad \quad \vdots \\ (\theta_{i+4}-\theta)^{7}\,, & \quad \text{where} \ \ \theta \in [\theta_{i+3},\theta_{i+4}]\,, \\ 0\,, & \quad \text{otherwise}\,, \end{aligned} \right. $$
where $i \in [-3,n+3]$, we obtain

$$ u_{i}(\theta)=c_{i-3}+120 c_{i-2}+1191 c_{i-1}+2416 c_{i} \\ +1191 c_{i+1}+120 c_{i+2}+c_{i+3}\,. $$
Substituting Eq. (39), its derivatives and boundary conditions of Eq. (25) into Eq. (15), gets $(n+7)$ of equations. Solving this system of equations to get the value of $c_{i}$, we obtain the following Table 4.


Table 4
Table 4Using septic spline technique to evaluate exact and approximate solutions of Eq. (15) at different values of $\theta$ to show the absolute error of them.

New window|CSV

4 Physical Explanation of Tables and Figures

In this section, we interest to clear of physical properties of our represented figures and tables.

4.1 Interpretation of Shown Tables

This part shows the mathematical and physical meaning of the shown tables as follows:

(i) Table 1 represents the convergence between the values of exact and approximate solutions. Moreover, it shows the fact of this method that is be more accurate when the value of $\theta$ near to zero.

(ii) Table 2 gives the values of exact and approximate solutions at different values of $\theta$ that show the values of the absolute value of error at the parties to the period are equal to zero and become bigger when taking away from the parties to the period until to reach to the middle of this period and come back again to decrease.

(iii) Table 3 aims to show the values of exact and approximate solutions for making a discussion about absolute error that is remarked that is increasing when the value of $\theta$ increase.

(iv) Table 4 aims to show the values of exact and approximate solutions for making a discussion about the absolute error that is remarked that is increasing when the value of $\theta$ increase.

4.2 Interpretation of Represented Figures

This subsection gives an interpretation of the shown figures in our paper.

(i) Figures 1--4 represent three and two-dimensional plots of Eqs. (25) and (26), when $( \delta =-4$, $h_2=1$, $h_6=4$, $h_4=4$, $r=2$, $\epsilon =1$, $\vartheta =3 )$.

Fig. 1

New window|Download| PPT slide
Fig. 1(Color online) Three-dimensional sketches for absolute, real, and imaginary solutions of Eq. (25), respectively.



Fig. 2

New window|Download| PPT slide
Fig. 2(Color online) Two-dimensional sketches for absolute, real, and imaginary solutions of Eq. (25), respectively.



Fig. 3

New window|Download| PPT slide
Fig. 3(Color online) Three-dimensional sketches for absolute, real, and imaginary solutions of Eq. (26), respectively.



Fig. 4

New window|Download| PPT slide
Fig. 4(Color online) Two-dimensional sketches for absolute, real, and imaginary solutions of Eq. (26), respectively.



(ii) Figure 5 shows the analytical and approximate solutions of Eqs. (25) and (31) via two and three-dimensional plots. We also represent the contour plot of approximate solutions (31).

Fig. 5

New window|Download| PPT slide
Fig. 5(Color online) List line plot in two and three-dimensional of Eqs. (25) and (31) and list contour plot of approximate solution (31).



(iii) Figure 6 represents the numerical that obtained by cubic spline and solitary solution of Eq. (14) via specialized radar, waterfall and verticals panel curves.

Fig. 6

New window|Download| PPT slide
Fig. 6(Color online) Specialized radar, waterfall and verticals panel curves between exact and approximate solutions (25), (33).



(iv) Figure 7 represents the line spline connected and specialized Radar plots for exact and numerical solutions that obtained by quintic spline.

Fig. 7

New window|Download| PPT slide
Fig. 7(Color online) Line spline connected and specialized Radar plots between exact Eq. (25) which represent by line (B), approximate Eq. (36) which represent by line (C) solutions and Error (B).



(v) Figure 8 represents the list line plot, Pareto chart-Binned data, and scatter matrix curves for exact and numerical solutions that obtained by septic spline.

Fig. 8

New window|Download| PPT slide
Fig. 8(Color online) List line plot, Pareto chart-Binned data and scatter matrix curves between exact and approximate solutions (25), (39).



5 Conclusion

This part of our paper gives a summary of the basic points in our paper. Analytical and numerical schemes have been successfully utilized to establish exact and approximate solutions of the complex NLS equation. Based on these five techniques, we got plenty of exact solutions that compared with four different approximate solutions to calculate the value of absolute error. Based on the values of the absolute errors that shown in Figs. 9 and 10, we see that the Adomian decomposition method, quintic, and septic schemes obtained more accurate approximate solutions than cubic spline scheme. We also plot some figures to show how the difference between exact and approximate solutions are very tiny such that in some figures, we cannot recognize the difference between them. The obtained results will be helpful in different fields of science.

Fig. 9

New window|Download| PPT slide
Fig. 9(Color online) Representation of absolute error for Tables 1--2.



Fig. 10

New window|Download| PPT slide
Fig. 10(Color online) Representation of absolute error for Tables 3--4.




Reference By original order
By published year
By cited within times
By Impact factor

F. Z. Geng, S. P. Qian, S. Li , J. Journal of Computational and Applied Mathematics 255(2014) 97.
[Cited within: 1]

P. K. Sahu and S. Saha Ray , Mathematical Methods in the Applied Sciences 38 ( 2015) 274.


R. Pakrouh, M. J. Hosseini, A. A. Ranjbar, R. Bahrampoury , Energy Conversion and Management 103 ( 2015) 542.


A. N. Hirani, K. B. Nakshatrala, J. H. Chaudhry, International Journal for Computational Methods in Engineering Science and Mechanics 16 ( 2015) 151.


M. M. Bhatti, M. A. Abbas, M. M. Rashidi , Applied Mathematics and Computation 316 ( 2018) 381.


A. Amir, J. Selvaraj, N. A. Rahim, Renewable and Sustainable Energy Reviews 62 ( 2016) 350.


K. Vasilic, W. Schmidt, H. C. Kühne , et al., Cement and Concrete Research 88 ( 2016) 1.


B. Yu, X. Jiang, H. Xu , Numerical Algorithms 68 ( 2015) 923.


R. A. Attia, D. Lu, M. M. A. Khater, Mathematical and Computational Applications 24 ( 2019) 10.


M. M. A. Khater, D. Lu, R. A. M. Attia , AIP Advances 9 ( 2019) 025003.


M. M. A. Khater, R. A. M. Attia, D. Lu, Mathematical and Computational Applications 24 ( 2019) 1.


R. A. M. Attia, D. Lu, M. M. A. Khater , Phys. J. ( 2018) 1.


H. Bulut, T. A. Sulaiman, H. M. Baskonus, Superlattices and Microstructures 123 ( 2018) 12.


M. S. Osman and A. M. Wazwaz , Applied Mathematics and Computation 321 ( 2018) 282.


A. M. Wazwaz and M. S. Osman , Computers & Mathematics with Applications 76 ( 2018) 276.
[Cited within: 1]

L. Demkowicz, J. Gopalakrishnan, S. Nagaraj, P. Sepulveda , SIAM Journal on Numerical Analysis 55 ( 2017) 1740.
[Cited within: 1]

E. K. Sklyanin and L. D. Faddeev , In Fifty Years of Mathematical Physics: Selected Works of Ludwig Faddeev
World Scientific Publishing Company (2016) pp. 290-292.



X. Zhao, Z. Z. Sun, Z. P. Hao , SIAM Journal on Scientific Computing 36 ( 2014) A2865.


C. Q. Dai, J. Liu, Y. Fan, D. G. Yu , Nonlinear Dynamics 88 ( 2017) 1373.


H. Triki and A. M. Wazwaz , Waves in Random and Complex Media 27 ( 2017) 587.


M. Eslami , Optik-International Journal for Light and Electron Optics 126 ( 2015) 3987.


H. Triki and A. M. Wazwaz , Rom. J. Phys. 61 ( 2016) 360.


P. Felmer and C. Torres , Commun. Pure Appl. 1 ( 2014) 2.


A. Maireche , The African Rev. Phys. 10 ( 2015) 373.


A. Maireche , The African Rev. Phys. 9 ( 2014) 185.
[Cited within: 1]

Sheng Zhang , Computers & Mathematics with Applications 54 ( 2007) 1028.
[Cited within: 1]

S. Jiong , Phys. Lett. A 309 ( 2003) 387.


C. P. Liu and X. P. Liu , Phys. Lett. A 348 ( 2006) 222.
[Cited within: 1]

J. Rashidinia and S. Sharifi , International Journal of Mathematical Modelling & Computations 5 ( 2015) 111.
[Cited within: 1]

M. Donatelli, C. Garoni, C. Manni , et al., SIAM Journal on Numerical Analysis 55 ( 2017) 31.
[Cited within: 1]

相关话题/Analytical Approximate Solutions