Comparative Solution of Nonlinear Quintic Cubic Oscillator Using Modified Homotopy Perturbation Method

We use modified homotopy perturbation method to find out the solution of nonlinear cubic quintic equation. Besides this method solution of the problemwith the followingmethods is discussed, Energy BalanceMethod andHe’s Frequency Formulationmethod and then compare the results with each other and Global Error Method. The results show that these three methods are effective as global error method for nonlinear cubic quintic oscillator equation with multiple nonlinear terms and have different effects on the solution. In particular, the homotopy perturbation solution is quite surprising. A cubic quintic nonlinear oscillator is used as an example to compare the results.


Introduction
In the present era large numbers of engineering problems are nonlinear and it is almost difficult to solve them analytically.Perturbation methods are among the known methods to solve the nonlinear problems, which are based on the existence of small/large parameters so-called perturbation parameters.These so-called perturbation parameters are not present in many nonlinear problems so we cannot use perturbation methods.To overcome this difficulty researchers introduced such methods which do not require small parameters such as homotopy perturbation method (HPM) [1,2], book keeping parameter perturbation method [3], He's energy balance method (EBM) [4][5][6][7], global error minimization method [8], Variational Iterational method [4], and many others.In this paper, solution of nonlinear cubic quintic oscillatory equation is solved with iterative methods, the basic idea of the homotopy perturbation is used, and its applications in some nonlinear oscillator equations are studied.In most engineering problems, it is easy to find angular frequency of a solution to a nonlinear equation.He presented the method for the first time.These results are valid for weakly nonlinear systems as well as strongly nonlinear systems.In Section 2, solution with energy balance method is explained.Solution with He's frequency formulation method is explained in Section 3. In Section 4, HPM solution is explained and some conclusions are given in Section 5.
Being the trial function, substituting (5) into (4), we obtain the following residual equation: If, by chance, the exact solution had been chosen as the trial function, then it would be possible to make residual () zero for all values of "" by appropriate choice of .Since ( 5) is only an approximation to the exact solution, () cannot be made zero everywhere.

He's Frequency Formulation Method
Equation ( 1) can be written in the general form We use two trial functions  1 () = cos and  2 () =  cos , which are, respectively, the solutions of the following linear oscillator equations: where  is assumed to be the frequency of the nonlinear oscillator.The residuals for (8) are The original frequency-amplitude formulation reads [9] We just used the following formulation: Geng and Cai improved the formulation by choosing another location point [10]: The accuracy depends upon the chosen location point, and we have no general rule for the choice.
To illustrate this shortcoming, we consider the cubic quintic nonlinear equation [8]: The residual of the equations is where frequency can be calculated at 0 or /3: Now locating cos  = cos  =  (where 0 <  < 1,  = 1, for (15) and  = 1/2 for (16)), we obtain It shows that the accuracy depends on the point of location.
Advances in Mathematical Physics 3 3.1.The Method of Weighted Residual.We use the method of weighted residuals to overcome the shortcoming.To this end, we introduce two new residual variables R1 and R2 defined as We can approximately calculate  2 in the form Now, by using (20) for frequency formulation, we find the value of  2 for cubic quintic nonlinear oscillator: Its approximate frequency reads

Modified Homotopy Perturbation Method with Two Expanding Parameters
Consider the following nonlinear cubic quintic oscillator [8]: where  1 and  2 are constant.Equation ( 23) is exactly solvable and its solution is given in [11].According to homotopy perturbation method [1,2], homotopy equation can be constructed in the form where  is unknown frequency of oscillator to be determined.When  = 0, we have This can be approximately describing the solution properties.Embedding parameter  increases monotonically from zero to unit as the trivial problem.Equation ( 25) is continuously deformed to the original one.According to standard homotopy perturbation method [1,2], solution is expanded in a series of : Generally one iteration is enough to obtain an approximate solution with relatively high accuracy.If higher accuracy of solution is needed, homotopy equation is suggested in the following [1]: Solution is also expanded in the same form of ( 26) and coefficient 1 of the linear term in ( 27) is also expanded in series of : This method is also used in parameter-expansion method.In [6] a homotopy equation with an auxiliary term is suggested which vanishes for both cases with homotopy parameters  = 0 and  = 1 in constructed homotopy equation.
A detailed solution process is given and a relatively comprehensive survey on concepts, theory, and applications of homotopy perturbation method is reported [1][2][3]12].In this paper, we suggest a homotopy perturbation method with two expanding parameters.

Solution of Problem
Using HPM.In (23) there are two nonlinear terms, which might have different effects on the solution.In order to find the best match with different nonlinear terms, we can construct a homotopy equation in the form where  1 and  2 are homotopy parameters, Solution and coefficient 1 of linear term in (29) are expanded in the form This expansion is similar to modified Lindstedt-Poincare method with double series expansion [3].Substituting (30) and ( 31) into (29), collecting the same power of

Advances in Mathematical Physics
1   2 (,  = 0, 1, 2, 3, . ..), and setting coefficients zero, we can readily obtain the following equations: Initial conditions for   ( > 1) can be freely chosen, but it must hold that (0) =  and   (0) = 0 for final closed approximate solution.Solving (32), we have Substitution of  0 into (33) results in Requiring no secular term in  1 from (39), we have A particular solution for (39) is readily obtained as Substituting  1 ,  0 , and  1 into (34) one gets In order to remove a secular term from (42), consider Now obtain a particular solution for (42): Substituting the obtained results into (35), we obtain equation for  3 and  3 , respectively: Similarly, we can find out the values of  4 and  5 by using (36) and (37): The solution process continues in this way.Substituting the values of   ( = 1 − 5) into (31) and using  1 = 1 and  2 = 1, we have According to Nayfeh's result [13] for comparison This result is obtained under the assumption of small amplitude.Consequently, the above methods give the same result as that of Nayfeh [13].In case of  2 being zero it will become the well-known Duffing equation and its nonlinear angular frequency can be obtained.The result is the same as that obtained in [13] and its maximal relative error is less than 7.5%.

Conclusions
From Table 1 it is suggested that these methods not only are effective for the given nonlinear cubic quintic oscillatory equation but also can be extended to various nonlinear problems with multiple nonlinear terms.In particular, the solution procedure of modified homotopy perturbation method with two expanding parameters can be used as a paradigm for many other applications.

Table 1 :
Comparison of HPM, He's frequency, energy balance method, global error method, and exact solution.