On the periodic solutions of the nonlinear oscillators
Abdelhalim Ebaid^{1} , Shaher Momani^{2} , ShihHsiang Chang^{3} , Mona D. Aljoufi^{4}
^{1, 4}Department of Mathematics, Faculty of Science, University of Tabuk, P. O. Box 741, Tabuk 71491, Saudi Arabia
^{2}Department of Mathematics, Faculty of Science, University of Jordan, Amman 11942, Jordan
^{3}Department of Mechanical Engineering, Far East University, Jhonghua Road, Tainan 74448, Taiwan, R. O. C.
^{1}Corresponding author
Journal of Vibroengineering, Vol. 16, Issue 1, 2014, p. 114.
Received 10 July 2013; received in revised form 29 October 2013; accepted 5 November 2013; published 15 February 2014
JVE Conferences
In this paper, a new approach is introduced to overcome the difficulty of applying the differential transformation method to the nonlinear oscillators described by $\ddot{x}\left(t\right)+f\left(x,\dot{x}\left(t\right),\leftx\left(t\right)\right\right)=0$. The obtained approximate periodic solutions are compared with those in open literatures and the results reveal that the present approach is very effective and convenient for a class of nonlinear oscillators with discontinuities.
Keywords: nonlinear oscillators, periodic solutions, differential transformation method, Cosineaftertreatment technique.
1. Introduction
Since the beginning of 1986, Zhou [1] and Pukhov [2] have developed a socalled differential transformation method (DTM) for electrical circuits problems. The method has been used extensively to solve effectively various linear and nonlinear ordinary, partial and integrodifferential equations [316]. However this method has a difficulty when applied to differential equations with discontinuous terms. So the main purpose of the present research is not only to solve this difficulty but also to show how to obtain the approximate periodic solutions for nonlinear oscillators with discontinuities. To our best knowledge, there is no paper reported in the literature on the application of the differential transformation method to differential equations with discontinuities.
It is well known that DTM constructs the solution in the form of a truncated series which is periodic only in a very small region [17]. In order to overcome this difficulty, a reliable aftertreatment (AT) technique has been developed recently in [18] to obtain the approximate periodic solutions in a wider range. As shown in [18], the proposed aftertreatment technique splits into two types, named as (SineAT technique, SAT) and (CosineAT technique, CAT) and for illustration we will reintroduce the basic idea of the CAT technique which will be used in this paper. It was shown in [18] that the approximate periodic solutions have been obtained without any need for Pade approximants or Laplace transform which may give the SAT and CAT techniques some power if compared with the modified differential transform method [17, 19]. Very recently Merdan and Gokdoan [19] showed that the CosineAT and SineAT give same results obtained by using the modified differential transform method and the classical fourthorder RungeKutta (RK4) method.
The rest of this paper is organized as follows: in the next section (Section 2) a brief description of the onedimensional DTM is provided; in Sections 3 the basic idea of the Cosineaftertreatment technique is discussed; in Section 4 the proposed techniques are implemented to obtain the approximate periodic solutions for three nonlinear oscillators with discontinuities; in Sections 5 and 6 some conclusions are given.
2. Onedimensional differential transform
The differential transform of a function $x\left(t\right)$ is defined as follows:
In Eq. (1) $x\left(t\right)$ is the original function and $X\left(k\right)$ is the transformed function. The differential inverse transform of $X\left(k\right)$ is defined as:
So:
Eq. (3) implies that the concept of differential transform is derived from the Taylor series expansion. In actual applications the function $x\left(t\right)$ is expressed by a truncated series and Eq. (2) can be written as:
Some of the fundamental mathematical operations performed by the differential transform method are listed in Table 1.
Table 1. The fundamental operations of onedimensional DTM
Original function $x\left(t\right)$

Transformed function $X\left(k\right)$

$\alpha u\left(t\right){}_{}{}^{+}\beta v\left(t\right)$

$\alpha U\left(t\right){}_{}{}^{+}\beta V\left(t\right)$

$\frac{{d}^{m}x\left(t\right)}{{dt}^{m}}$

$\frac{\left(k+m\right)!}{k!}X(k+m)$

$u\left(t\right)v\left(t\right)$

$\sum _{I=0}^{K}U\left(I\right)V(kI)$

$u\left(t\right)v\left(t\right)w\left(t\right)$

$\sum _{I=0}^{k}\sum _{m=0}^{kI}U\left(I\right)V\left(m\right)W(kIm)$

$u\left(t\right)\underset{0}{\overset{t}{\int}}v\left(t\right)dt$

$\sum _{I=1}^{k}U\left(kI\right)\frac{V(I1)}{I},k\ge 1$

${t}^{m}$

$\delta \left(km\right)=1\text{if}k=m,0\text{if}k\ne m$

${e}^{x}$

$\frac{1}{k!}$

$\mathrm{sin}(\lambda t+\omega )$

$\frac{{\lambda}^{k}}{k!}\mathrm{sin}\left(\frac{k\pi}{2}+\omega \right)$

$\mathrm{cos}(\lambda t+\omega )$

$\frac{{\lambda}^{k}}{k!}\mathrm{cos}\left(\frac{k\pi}{2}+\omega \right)$

3. Cosineaftertreatment technique (CATtechnique)
If the truncated series given by Eq. (4) is expressed only in evenpowers of the independent variable $t$, i. e.:
then the Cosineaftertreatment technique (CATtechnique) is based on the assumption that this truncated series can be expressed as another finite series in terms of the cosine trigonometric functions with different amplitudes and frequencies:
where $n$ is finite.
Our main question here is how to find the values of ${\lambda}_{j}$ and ${\Omega}_{j}$, $j=0,\dots ,n.$ For answer we begin by expanding both sides of Eq. (6) as power series of $t$ to obtain:
Now the equation of the coefficients of like powers yields:
${t}^{4}:\sum _{j=1}^{n}{\lambda}_{j}{\Omega}_{j}^{4}=4!X\left(4\right),{t}^{6}:\sum _{j=1}^{n}{\lambda}_{j}{\Omega}_{j}^{6}=6!X\left(6\right).$
In practical applications it is sufficient to express the truncated series ${\Phi}_{N}\left(t\right)$ in terms of two or three cosines with different amplitudes and frequencies. If we choose to express ${\Phi}_{N}\left(t\right)$ as an approximate periodic solution in terms of two cosines with two different amplitudes, ${\lambda}_{1}$and ${\lambda}_{2},$ and two different frequencies, ${\Omega}_{1}$ and${\Omega}_{2}$, we can rewrite Eq. (6) for $N=6$ and $n=2$ as:
In this case the four unknowns ${\lambda}_{1},{\lambda}_{2},{\Omega}_{1}$ and ${\Omega}_{2}$ can be determined by solving the following system of nonlinear algebraic equations analytically:
${\lambda}_{1}{\Omega}_{1}^{2}+{\lambda}_{2}{\Omega}_{2}^{2}=2!X\left(2\right),$
${\lambda}_{1}{\Omega}_{1}^{4}+{\lambda}_{2}{\Omega}_{2}^{4}=4!X\left(4\right),$
${\lambda}_{1}{\Omega}_{1}^{6}+{\lambda}_{2}{\Omega}_{2}^{6}=6!X\left(6\right).$
Moreover, if we choose to express ${\Phi}_{N}\left(t\right)$ as more accurate periodic solution in terms of three cosines, we can rewrite Eq. (6) for $N=10$ and $n=3$ as:
In this case the six unknowns ${\lambda}_{1},{\lambda}_{2},{\lambda}_{3},{\Omega}_{1},\mathrm{}{\Omega}_{2}$ and ${\Omega}_{3}$ can be determined by solving the following system of nonlinear algebraic equations numerically:
${\lambda}_{1}{\Omega}_{1}^{2}+{\lambda}_{2}{\Omega}_{2}^{2}+{\lambda}_{3}{\Omega}_{3}^{2}=2!X\left(2\right),$
${\lambda}_{1}{\Omega}_{1}^{4}+{\lambda}_{2}{\Omega}_{2}^{4}+{\lambda}_{3}{\Omega}_{3}^{4}=4!X\left(4\right),$
${\lambda}_{1}{\Omega}_{1}^{6}+{\lambda}_{2}{\Omega}_{2}^{6}+{\lambda}_{3}{\Omega}_{3}^{6}=6!X\left(6\right),$
${\lambda}_{1}{\Omega}_{1}^{8}+{\lambda}_{2}{\Omega}_{2}^{8}+{\lambda}_{3}{\Omega}_{3}^{8}=8!X\left(8\right),$
${\lambda}_{1}{\Omega}_{1}^{10}+{\lambda}_{2}{\Omega}_{2}^{10}+{\lambda}_{3}{\Omega}_{3}^{10}=10!X\left(10\right).$
4. Applications
In this section we discuss the application of the DTM with the proposed CosineAT technique to obtain the periodic solutions of nonlinear oscillators with discontinuities. As explained in the previous section, the use of the CosineAT technique is based on obtaining the series solution of the problem under consideration as a finite polynomial in even powers of $t$.
4.1. Example 1
Consider the nonlinear oscillator with a discontinuous term [2023]:
with the initial conditions:
In Eq. (13) the major source of the difficulty in using the differential transformation method is the existence of the discontinuous term $\leftx\right.$ In order to overcome this difficulty we introduce a new approach to deal with such term. Firstly we suppose that:
which can be written as:
We then take a differentiation with respect to $t$ to obtain:
By this Eq. (13) becomes:
Now by applying the differential transform to Eqs. (14), (17) and (18) we obtain the following recurrence scheme for $k\ge 0$:
$F\left(0\right)=\leftx\left(0\right)\right=a,$
$\left(k+1\right)\left(k+2\right)X\left(k+2\right)+\sum _{m=0}^{k}X\left(km\right)F\left(m\right)=0,$
$\sum _{m=0}^{k}\left(m+1\right)\left[F\left(m+1\right)F\left(km\right)X\left(m+1\right)X\left(km\right)\right]=0.$
Using Eqs. (19) by taking $N=6$, we obtain a system of algebraic equations for $k=0,\dots ,4$. By solving this system for the values of $X\left(2\right),\dots ,X\left(6\right)$ by using Mathematica, we get:
From the inverse transformation rule Eq. (4) we can construct the following truncated series solution up to ${t}^{6}$:
Since the truncated series ${\Phi}_{6}\left(t\right)$ given by the last equation is expressed as a polynomial in even powers of $t$, we can now deal with it by using the CosineAT technique. We assume an approximate periodic solution in the form:
In order to find ${\lambda}_{1},{\lambda}_{2},{\Omega}_{1}$ and ${\Omega}_{2},$ we insert $X\left(0\right),X\left(2\right),X\left(4\right)$ and $X\left(6\right)$ presented above into system Eq. (10) to get a system of four nonlinear algebraic equations. By solving this system analytically for ${\lambda}_{1},{\lambda}_{2},{\Omega}_{1}$ and ${\Omega}_{2},$ we obtain:
${\Omega}_{1}=\sqrt{\left(4+\sqrt{10}\right)a},{\Omega}_{2}=\sqrt{\left(4\sqrt{10}\right)a}.$
Therefore we can write the approximate periodic solution for the Eqs. (1314) as:
$\times \mathrm{cos}\left(\sqrt{\left(4\sqrt{10}\right)a}t\right).$
Here it should be noted that this approximate periodic solution is valid regardless of the oscillation amplitude $a$, i. e., for $0<a<+\infty $. In [20] the authors applied the parameter expansion method to Eqs. (1314) and the following approximate periodic solution has been obtained:
Eq. (25) has been also obtained by the authors in [21] and [23] by using the maxmin approach and He’s amplitudefrequency approach, respectively. The exact displacement of Eqs. (1314) is given by [21]:
where $B(*,*)$ denotes the Beta function.
In order to check the effectiveness of our approach we compare our result Eq. (24) with the approximate periodic solution obtained by the parameterexpansion and the maxmin methods, given by Eq. (25) and the exact displacement Eq. (26) in Figs. 14 at different values of $a=1,10,100,1000$. From these figures it can be concluded that the suggested approach with the CosineAT technique leads to approximate periodic solutions with good accuracy for all amplitudes.
Fig. 1. Comparison of the present approach, Eq. (24) with the parameterexpansion method [20], Eq. (25) and the exact solution (maxmin [21]), Eq. (26) for example 1, at $a=1$
Fig. 2. Comparison of the present approach, Eq. (24) with the parameterexpansion method [20], Eq. (25) and the exact solution (maxmin [21]), Eq. (26) for example 1, at $a=10$
Fig. 3. Comparison of the present approach, Eq. (24) with the parameterexpansion method [20], Eq. (25) and the exact solution (maxmin [21]), Eq. (26) for example 1, at $a=100$
Fig. 4. Comparison of the present approach, Eq. (24) with the parameterexpansion method [20], Eq. (25) and the exact solution (maxmin [21]), Eq. (26) for example 1, at $a=1000$
4.2. Example 2
Consider the nonlinear oscillator with a discontinuous term [2425]:
with the initial conditions:
By the same analysis as in the previous example, we can easily obtain the following recurrence scheme for $k\ge 0$:
$F\left(0\right)=\leftx\left(0\right)\right=a,$
$\left(k+1\right)\left(k+2\right)X\left(k+2\right)+X\left(k\right)+\u03f5\sum _{m=0}^{k}X\left(km\right)F\left(m\right)=0,$
$\sum _{m=0}^{k}\left(m+1\right)\left[F\left(m+1\right)F\left(km\right)X\left(m+1\right)X\left(km\right)\right]=0.$
Using the above relations by taking $N=6$, we obtain a system of algebraic equations for $k=0,\dots ,4$. By solving this system for the values of $X\left(2\right),\dots ,X\left(6\right)$, we get:
$X\left(6\right)=\frac{1}{720}\left(1+11\u03f5a+20{\u03f5}^{2}{a}^{2}+10{\u03f5}^{3}{a}^{3}\right)a.$
Hence we have the following truncated series solution up to ${t}^{6}$:
$\frac{1}{720}\left(1+11\u03f5a+20{\u03f5}^{2}{a}^{2}+10{\u03f5}^{3}{a}^{3}\right)a{t}^{6}.$
Here ${\Phi}_{6}\left(t\right)$ is also expressed as a polynomial in even powers of $t$ and as shown in the previous example we can deal with it by using the CosineAT technique. We also assume an approximate periodic solution in the form of Eq. (22). This leads to the approximate periodic solution:
$\times \mathrm{cos}\left(\sqrt{\frac{1}{2}\left(7+8\u03f5a\sqrt{25+64\u03f5a+40{\u03f5}^{2}{a}^{2}}\right)}t\right)$
$+\frac{\left(56\u03f5a+\sqrt{25+64\u03f5a+40{\u03f5}^{2}{a}^{2}}\right)a}{2\sqrt{25+64\u03f5a+40{\u03f5}^{2}{a}^{2}}}$
$\times \mathrm{cos}\left(\sqrt{\frac{1}{2}\left(7+8\u03f5a+\sqrt{25+64\u03f5a+40{\u03f5}^{2}{a}^{2}}\right)}t\right),$
which is also periodic for any given values of $\u03f5$ and $a$. The approximate periodic solution obtained by the variational iteration method [24] is given by:
Fig. 5. Comparison of the present approach, Eq. (32) with the variational iteration method, Eq. (33) and the numerical solution for example 2, at $\u03f5=0.1$ and $a=1$
Fig. 6. Comparison of the present approach, Eq. (32) with the variational iteration method, Eq. (33) and the numerical solution for example 2, at $\u03f5=0.2$ and $a=10$
Fig. 7. Comparison of the present approach, Eq. (32) with the variational iteration method, Eq. (33) and the numerical solution for example 2, at $\u03f5=0.3$ and $a=100$
Fig. 8. Comparison of the present approach, Eq. (32) with the variational iteration method, Eq. (33) and the numerical solution for example 2, at $\u03f5=0.01$ and $a=1000$
In Figs. 58 we compared our analytical approximate periodic solution in Eq. (32) with the variational iteration result Eq. (33) and the numerical solution obtained by Mathematica using the “NDSolve” command. It can be observed from these figures that our result in Eq. (32) provides excellent agreement with variational iteration method and the numerical solution regardless of the oscillation amplitude $a$.
4.3. Example 3
This example considers the following equation [24, 25]:
with the initial conditions:
where the parameter $\beta $ is considered real. Proceeding as above, we obtain for $k\ge 0$:
$F\left(0\right)=\leftx\left(0\right)\right=a,$
$\left(k+1\right)\left(k+2\right)X\left(k+2\right)+\beta \sum _{m=0}^{k}\sum _{l=0}^{m}X\left(km\right)X\left(ml\right)X\left(l\right)+\u03f5\sum _{m=0}^{k}X\left(km\right)F\left(m\right)=0,$
$\sum _{m=0}^{k}\left(m+1\right)\left[F\left(m+1\right)F\left(km\right)X\left(m+1\right)X\left(km\right)\right]=0.$
Using $N=6$, a system of algebraic equations is obtained for $k=0,\dots ,4$. On solving this system for $X\left(2\right),\dots ,X\left(6\right),$ yields:
$X\left(5\right)=0,X\left(6\right)=\frac{1}{720}\left(27{\beta}^{3}{a}^{3}+63\u03f5{\beta}^{2}{a}^{2}+46{\u03f5}^{2}\beta a+10{\u03f5}^{3}\right){a}^{4}.$
Therefore the approximate series solution can be written as:
$\frac{1}{720}\left(27{\beta}^{3}{a}^{3}+63\u03f5{\beta}^{2}{a}^{2}+46{\u03f5}^{2}\beta a+10{\u03f5}^{3}\right){a}^{4}{t}^{6},$
which is also expressed as a polynomial in even powers of $t$. The CosineAT technique leads to the approximate periodic solution:
where ${\mu}_{1},{\mu}_{2},\sigma $ and $\tau $ are given as:
$\sigma =432{\beta}^{4}{a}^{4}+1080\u03f5{\beta}^{3}{a}^{3}+937{\u03f5}^{2}{\beta}^{2}{a}^{2}+328{\u03f5}^{3}\beta a+40{\u03f5}^{4}.$
In [24] the approximate periodic solution was obtained by using the variational iteration method and given by:
It should be noted that the solution Eq. (39) is always periodic under the two conditions:
and:
Examples for the values of $\beta $ and $\u03f5$ that satisfy these conditions are shown graphically by Fig. 9 at different values of the amplitude $a$. The numerical results are depicted and compared with other solutions in Figs. 1013 at different values of $a$. The results show that at $a=1$ the Cosineafter treatment technique agrees with the other methods in a wider range of $t.$ However this agreement decreases with increasing $a$. It may be concluded that at large values of $a$ we have to increase the accuracy of our technique by increasing $n$ in Eq. (6).
Fig. 9. Examples for the values of $\beta $ and $\u03f5$
Fig. 10. Comparison of the present Cosine aftertreatment technique, CATTech., Eq. (39) with the variational iteration method, Eq. (41) for example 3, at $a=1$, $\beta =0.5$ and $\u03f5=4$
Fig. 11. Comparison of the present Cosine aftertreatment technique, CATTech., Eq. (39) with the variational iteration method, Eq. (41) for example 3, at $a=10$, $\beta =1$ and $\u03f5=3$
Fig. 12. Comparison of the present Cosine aftertreatment technique, CATTech., Eq. (39) with the variational iteration method, Eq. (41) for example 3, at $a=100$, $\beta =3$ and $\u03f5=4$
Fig. 13. Comparison of the present Cosine aftertreatment technique, CATTech., Eq. (39) with the variational iteration method, Eq. (41) for example 3, at $a=1000$, $\beta =3$ and $\u03f5=4$
5. Remarks
In several papers [2630] many authors showed that the DTM is nothing more than a disguised version of the Taylor series method. In this regard we agree with their point of view in the case when the DTM and Taylor series are applied to search for a series solution for a differential equation in which the discontinuous terms are not involved. However, it may be difficult to directly apply Taylor series method to obtain the series solutions for the present class of ordinary differential equations with discontinuous terms.
6. Conclusions
To our best knowledge, this is the first paper reported on the application of the differential transformation method to nonlinear oscillators with discontinuities. A new approach has been proposed to overcome the difficulty arising from the existence of the discontinuous terms. Three examples are solved and the obtained results demonstrate that the discontinuous function will not affect much the effectiveness and convenience of the differential transformation method.
References
 J. K. Zhou Differential Transformation and Its Applications for Electrical Circuits. Huazhong University Press, Wuhan, China, 1986. [CrossRef]
 Pukhov G. E. Differential Transformations and Mathematical Modelling of Physical Processes. Kiev, 1986. [CrossRef]
 Shing Huei Ho, Cha’o Kuang Chen Analysis of general elastically end restrained nonuniform beams using differential transform. Applied Mathematical Modelling, Vol. 22, 1998, p. 219234. [CrossRef]
 MingJyi Jang, ChiehLi Chen, YungChin Liy On solving the initialvalue problems using the differential transformation method. Appl. Math. Comput., Vol. 115, 2000, p. 145160. [CrossRef]
 Cha’o Kuang Chen, Shing Huei Ho Transverse vibration of a rotating twisted Timoshenko beam under axial loading using differential transform. International Journal of Mechanical Sciences, Vol. 41, 1999, p. 13391356. [CrossRef]
 Muhammet Koksal, Saadetdin Herdem Analysis of nonlinear circuits by using differential Taylor transform. Computers and Electrical Engineering, Vol. 28, 2002, p. 513525. [CrossRef]
 I. H. AbdelHalim Hassan Different applications for the differential transformation in the differential equations. Appl. Math. Comput., Vol. 129, 2002, p. 183201. [CrossRef]
 Aytac Arikoglu, Ibrahim Ozkol Solution of boundary value problems for integrodifferential equations by using differential transform method. Appl. Math. Comput., Vol. 168, 2005, p. 11451158. [CrossRef]
 A. S. V. Ravi Kanth, K. Aruna Solution of singular twopoint boundary value problems using differential transformation method. Phys. Lett. A, Vol. 372, 2008, p. 46714673. [CrossRef]
 Cha’o Kuang Chen, Shing Huei Ho Solving partial differential equations by twodimensional differential transform method. Appl. Math. Comput., Vol. 41, p. 171179. [CrossRef]
 MingJyi Jang, ChiehLi Chen, YungChin Liy Twodimensional differential transform for partial differential equations. Appl. Math. Comput., Vol. 121, 2001, p. 261270. [CrossRef]
 Abd ElHalim Ebaid Approximate periodic solutions for the nonlinear relativistic harmonic oscillator via differential transformation method. Commun. Nonlin. Sci. Numer. Simulat, Vol. 15, 2010, p. 19211927. [CrossRef]
 Fatma Ayaz Solutions of the system of differential equations by differential transform method. Appl. Math. Comput., Vol. 147, 2004, p. 547567. [CrossRef]
 ShihHsiang Chang, ILing Chang A new algorithm for calculating onedimensional differential transform of nonlinear functions. Appl. Math. Comput., Vol. 195, 2008, p. 799808. [CrossRef]
 A. S. V. Ravi Kanth, K. Aruna Twodimensional differential transform method for solving linear and nonlinear Schrodinger equations. Chaos, Solitons and Fractals, Vol. 41, Issue 5, 2009, p. 22772281. [CrossRef]
 A. S. V. Ravi Kanth, K. Aruna Differential transform method for solving the linear and nonlinear KleinGordon equation. Computer Physics Communications, Vol. 180, 2009, p. 708711. [CrossRef]
 Shaher Momani, Vedat Suat Erturk Solutions of nonlinear oscillators by the modified differential transform method. Computers and Mathematics with Applications, Vol. 55, 2008, p. 833842. [CrossRef]
 Abd Elhalim Ebaid A reliable aftertreament for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities. Commun. Nonlin. Sci. Numer. Simulat., Vol. 16, 2011, p. 528536. [CrossRef]
 M. Merdan, A. Gökdoan Solution of nonlinear oscillators with fractional nonlinearities by using the modified differential transformation method. Mathematical and Computational Applications, Vol. 16, Issue 3, 2011, p. 761772. [CrossRef]
 ShuQiang Wang, JiHuan He Nonlinear oscillator with discontinuity by parameterexpansion method. Chaos, Solitons and Fractals, Vol. 35, 2008, p. 688691. [CrossRef]
 DeQiang Zeng Nonlinear oscillator with discontinuity by the maxmin approach. Chaos, Solitons and Fractals, Vol. 42, 2009, p. 28852889. [CrossRef]
 D. Q. Zeng, Y. Y. Lee, C. K. Wong Analysis of a nonlinear oscillator with discontinuity. Comput. Math. Appl., Vol. 59, Issue 8, 2010, p. 25102515. [CrossRef]
 HuiLi Zhang Application of He’s amplitudefrequency formulation to a nonlinear oscillator with discontinuity. Comput. Math. Appl., Vol. 58, 2009, p. 21972198. [CrossRef]
 M. Rafei, D. D. Ganji, H. Daniali, H. Pashaei The variational iteration method for nonlinear oscillators with discontinuity. J. Sound. Vib., Vol. 305, 2007, p. 614620. [CrossRef]
 DaHua Shou Variational approach for nonlinear oscillators with discontinuities. Comput. Math. Appl., Vol. 58, 2009, p. 24162419. [CrossRef]
 C. Bervillier Status of the differential transformation method. Appl. Math. Comput., Vol. 218, 2012, p. 1015810170. [CrossRef]
 S. Abbasbandy, C. Bervillier Analytic continuation of Taylor series and the boundary value problems of some nonlinear ordinary differential equations. Appl. Math. Comput., Vol. 218, 2011, p. 21782199. [CrossRef]
 G. Corliss, Y. F. Chang Solving ordinary differential equations using Taylor series. ACM Trans. Math. Software, Vol. 8, 1982, p. 114. [CrossRef]
 R. Barrio, M. Rodriguez, A. Abad, F. Blesa Breaking the limits: the Taylor series method. Appl. Math. Comput., Vol. 217, 2011, p. 79407954. [CrossRef]
 D. H. Bailey, R. Barrio, J. M. Borwein Highprecision computation: mathematical physics and dynamics. Appl. Math. Comput., Vol. 218, 2012, p. 1010610121. [CrossRef]