Article International Journal of Modern Mathematical Sciences 2012 3(2): 63-76 International Journal of Modern Mathematical Sciences Journal homepage:wwwmodernscientificpresscom/journals/ijmmsaspx On Goursat Problems ISSN: 2166-286X Florida USA Muhammad Usman Tamour Zubair Umair Ali Syed Tauseef Mohyud-Din * Department of Mathematics HITEC University Taxila Cantt Pakistan * Author to whom correspondence should be addressed; Email: syedtauseefs@hotmailcom Article history: Received 29 May 2012 Received in revised form 1 August 2012 Accepted 8 August 2012 Published 13 August 2012 Abstract: In this paper we apply Homotopy Analysis Method (HAM) to find appropriate solutions of linear and nonlinear Goursat problems which are of utmost importance in applied and engineering sciences The proposed modification is the elegant coupling of Homotopy Analysis Method (HAM) Numerical results coupled with graphical representation explicitly reveal the complete reliability of the proposed algorithm Keywords: Homotopy analysis method linear and nonlinear Goursat problems exact solution Maple Mathematics Subject Classification (2000): 35Q79 1 Introduction The rapid development of nonlinear sciences [1-17] witnesses a wide range of analytical and numerical techniques by various scientists Most of the developed schemes have their limitations like limited convergence divergent results linearization discretization unrealistic assumptions and noncompatibility with the versatility of physical problems [1-11] In the similar context Liao [7-9] developed Homotopy Analysis Method (HAM) which is being applied on a wide range of nonlinear problems of physical nature see [1-17] and the references therein The basic motivation of present study is modification of traditional Homotopy Analysis Method (HAM) to tackle nonlinear partial differential equations It is observed that proposed technique is highly effective Moreover this method (HAM) is more user-friendly and it overcomes the complexities of selection of initial value Several examples are given which reveal the efficiency and reliability of the proposed algorithm
64 2 Analysis of Homotopy Analysis Method (HAM) We consider the following differential equation = 0 (1) where is a nonlinear operator denotes independent variables is an unknown function respectively For simplicity we ignore all boundary or initial conditions which can be treated in the similar way By means of generalizing the traditional Homotopy method Liao constructs the so called zero - order deformation equation 1 ; 0 = h; (2) where 01 is the embedding parameter h 0 is a nonzero parameter 0 is an auxiliary function is an auxiliary linear operator 0 is an initial guess of ; is a unknown function respectively It is important that one has great freedom to choose auxiliary things in HAM Obviously when and = 0 and = 1 it holds ;0 = ;1 = Respectively Thus as increases from 0 to 1 the solution ; varies from the initial guesses 0 to the solution Expanding in Taylor series with respect to we have ; = + where = 1 ;! at = 0 If the auxiliary linear operator the initial guess the auxiliary h and the auxiliary function are so properly chosen the above series converges at = 1 then we have Define the vector = + = { 0 1 2 0 } Differentiating equation (2) Times with respect to the embedding parameter and then setting = 0 and finally dividing them by! we obtain the %h-order deformation equation where and & ' = h( ' (3) ( 1 = 1 1 ; at = 0 1! 1 & = 0 1 = 1 > 1 Applying 1 both sides of (3) we get = & 1 + h 1 ( 1
65 In this way it is easily to obtain for m 1 at %h- order we have = when - we get an accurate approximation of the original equation(1) For the convergence of the above method we refer the reader Liao s work If equation (1) admits unique solution then this method will produce the unique solution If equation (1) does not possess unique solution the HAM will give a solution among many other (possible) solutions 3 Goursat Problems The Goursat problems arise in linear and nonlinear partial differential equations which mixed derivatives The slandered form of the Goursat problems is given by /% = 0/% / % 0 / 20 % 3 /0 = 4/ 0% = h% 40 = h0 = 00 The slandered form of the homogenous linear Goursat problems is given by /% = 0 0 / 20 % 3 /0 = 4/ 0% = h% 40 = h0 = 00 The slandered form of the inhomogeneous linear Goursat problems is given by /% = 0 + 5/% 0 / 20 % 3 /0 = 4/ 0% = h% 40 = h0 = 00 4 Numerical Application In this section we apply Homotopy analysis method (HAM) for linear and nonlinear Goursat problems Numerical results are very encouraging Example 41 We first consider the following homogenous linear Goursat problem /% = with initial conditions /0 = 6 7 0% = 6 8 00 = 1 To solve the given Equation by HAM we choose the linear operator with the property /%;9 = :; :7:8 </%;9= > + />? = 0
66 > + %>? = 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 = @ @ A/A% we now define a nonlinear operator as /%;9 = /%;9 78 /%;9 using the above definition we construct the zeroth-order deformation equation 1 9/%;9 0 /% = 9h/%;9/%;9 for 9 = 0 and 9 = 1 we can write /%;0 = /% /%;1 = /% Thus we obtain the %h order deformation equation /% & ' /% = h/%( ' with initial condition /0 = 0 00 = 1 where ( 1 = B 2 /% 1 1 C Now the solutions of the %h order deformation equation are /% = & 1 /% + 1 h/%( 1 1 we start with an initial approximation /% = 6 7 + 6 8 1 by means of the iteration formula as discuss above if h = 1 = 1 we can obtain directly the others components as the series solution is given by /% = 6 7 + 6 8 1 /% = 1 2 % 2 % + 268? /% = 1 2 %? + 26 8 % + 2 1 6 %E 1 24 %G 26 8 /% = 6 7 + 6 8 1 + 1 2 % 2 % + 268 1 2 %? + 26 8 % + 2 1 6 %E 1 24 %G 26 8 + and the close form solution is /% = 6 7I8
67 Figure 41 Depicts approximate and exact solutions Example 42 We first consider the following homogenous linear Goursat problem /% = 2 with initial conditions /0 6 7 0% 6 '?8 00 1 To solve the given Equation by HAM we choose the linear operator with the property /%;9 :; :7:8 </%;9= > />? 0 > %>? 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 @ @ A/A% we now define a nonlinear operator as /%;9 /%;9 78 2/%;9 using the above definition we construct the zeroth-order deformation equation 1 9/%;9 0 /% 9/%;9/%;9 for 9 0 and 9 1 we can write /%;0 /% /%;1 /% Thus we obtain the % order deformation equation /% & ' /% /%( ' with initial condition /0 0 00 1 where ( 1 B 2 /% 1 2 1 C
68 Now the solutions of the % order deformation equation are /% & 1 /% 1 /%( 1 + 1 we start with an initial approximation /% 6 7 6 '?8 1 by means of the iteration formula as discuss above if 1 1 we can obtain directly the others components as /% 6 7 6 '?8 1 /% 3 2 % %? 2%6 8 26 8 1 2 6'?8? /% 7 2% 2%? 26 8 % 66 8 6 '?8 2 3 %E the series solution is given by /% 6 7 6 '?8 1 3 2 % %? 2%6 8 26 8 1 2 6'?8 7 2% 2%? 26 8 % 66 8 6 '?8 2 3 %E and the close form solution is /%6 7'?8 Figure 42 Depicts approximate and exact solutions
69 Example 43 We first consider the following inhomogeneous linear Goursat problem /% % with initial conditions /0 = 6 7 0% = % + 6 8 00 = 1 To solve the given Equation by HAM we choose the linear operator with the property /%;9 = :; :7:8 </%;9= > + />? = 0 > + %>? = 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 = @ @ A/A% we now define a nonlinear operator as /%;9 = /%;9 78 /%;9 + % using the above definition we construct the zeroth-order deformation equation 1 9/%;9 0 /% = 9h/%;9/%;9 for 9 = 0 and 9 = 1 we can write /%;0 = /% /%;1 = /% Thus we obtain the %h order deformation equation /% & ' /% = h/%( ' with initial condition /0 = 0 00 = 1 where ( 1 = B 2 /% 1 1 + %C Now the solutions of the %h order deformation equation are /% = & 1 /% + 1 h/%( 1 1 we start with an initial approximation /% = 6 7 + % + 6 8 1 by means of the iteration formula as discuss above if h = 1 = 1 we can obtain directly the others components as 2 /% = 1 2 %2 + 2%6 % + 2 1 6 %3 1 24 %4 26 % 0 /% = 6 / + % + 6 % 1 1 /% = 1 2 % 2 % + 26%
70 the series solution is given by /% 6 7 % 6 8 1 1 2 %2 % 268 1 2 %? 2%6 8 2 1 6 %E 1 24 %G 26 8 and the close form solution is /%%6 7I8 Figure 43 Depicts approximate and exact solutions Example 44 We first consider the following inhomogeneous linear Goursat problem /% 4/%% 2 / 2 with initial conditions /06 7 0%6 8 001 To solve the given Equation by HAM we choose the linear operator with the property /%;9 :; :7:8 </%;9= > />? 0 > %>? 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 @ @A/A% we now define a nonlinear operator as /%;9/%;9 78 /%;94/%/? %? using the above definition we construct the zeroth-order deformation equation
71 1 9/%;9 0 /% = 9h/%;9/%;9 for 9 = 0 and 9 = 1 we can write /%;0 = /% /%;1 = /% Thus we obtain the %h order deformation equation /% & ' /% = h/%( ' with initial condition /0 = 0 00 = 1 where ( 1 = B 2 /% 1 + 1 4/% + / 2 % 2 C Now the solutions of the %h order deformation equation are /% = & 1 /% + 1 h/%( 1 1 we start with an initial approximation 0 /% = 6 / + 6 % 1 by means of the iteration formula as discuss above if h = 1 = 1 we can obtain directly the others components as 0 /% = 6 / + 6 % 1 1 /% = 1 18 % 18 9% + 186% + 9% 3 % 5 2 /% = 1 2 %2 + 2%6 % + 11 24 %4 7 180 %6 + 2 1 6 %3 26 % the series solution is given by /% = 6 7 + 6 8 1 + 1 18 % 18 9% + 1868 + 9% E % O 1 2 %? + 2%6 8 + 11 24 %G 7 180 %P + 2 1 6 %E 26 8 + and the close form solution is /% = %? /? + 6 7I8
72 Figure 44 Depicts approximate and exact solutions Example 45 We first consider the following inhomogeneous linear Goursat problem /% 3 / 3 3/ 2 % 3/% 2 % 3 with initial conditions /0 / 0% % 00 0 To solve the given Equation by HAM we choose the linear operator with the property /%;9 :; :7:8 </%;9= > />? 0 > %>? 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 @ @ A/A% we now define a nonlinear operator as /%;9 /%;9 78 E /%;9 / E 3/? % 3/%? % E using the above definition we construct the zeroth-order deformation equation 1 9/%;9 0 /% 9/%;9/%;9 for 9 0 and 9 1 we can write /%;0 /% /%;1 /% Thus we obtain the % order deformation equation /% & ' /% /%( ' with initial condition /0 0 00 1
73 Q where ( 1 B 2 /% 1 1 S0 R0 Q QR 1S / 3 3/ 2 % 3/% 2 % 3 C Now the solutions of the % order deformation equation are /% & 1 /% 1 /%( 1 + 1 we start with an initial approximation /% / % by means of the iteration formula as discuss above if 1 1 we can obtain directly the others components as 0 /% / % 1 /% 0 2 /% 0 the series solution is given by /% / % and the close form solution is /% / % Figure 45 Depicts approximate and exact solutions Example 46 We first consider the following inhomogeneous linear Goursat problem /% 2 6 2/ 6 2% 26 /%
74 with initial conditions /0 = 1 + 6 7 0% = 1 + 6 8 00 = 2 To solve the given Equation by HAM we choose the linear operator with the property /%;9 = :; :7:8 </%;9= > + />? = 0 > + %>? = 0 where > 1 and > 2 are the integral constants The inverse operator 1 is given by / % 0 0 1 = @ @ A/A% we now define a nonlinear operator as /%;9 = /%;9 78 +? /%;9 6?7 + 6?8 + 26 7I8 using the above definition we construct the zeroth-order deformation equation 1 9/%;9 0 /% = 9h/%;9/%;9 for 9 = 0 and 9 = 1 we can write /%;0 = /% /%;1 = /% Thus we obtain the %h order deformation equation /% & ' /% = h/%( ' with initial condition /0 = 0 00 = 1 where ( 1 = B 2 /% 1 + 1 U S=0 T=0 T U T 6 2/ + 6 2% + 26 /+% C Now the solutions of the %h order deformation equation are /% = & 1 /% + 1 h/%( 1 1 we start with an initial approximation 0 /% = 6 % + 6 / by means of the iteration formula as discuss above if h = 1 = 1 we can obtain directly the others components as 2 /% = 0 the series solution is given by 0 /% = 6 % + 6 / 1 /% = 0 /% = 6 8 + 6 7
75 and the close form solution is /%6 8 6 7 Figure 46 Depicts approximate and exact solutions 5 Conclusion In this paper we apply Homotopy analysis method on linear and nonlinear Goursat problems The result which obtained is very good and reliable as compared to other methods The advantages of HAM are illustrated It is clear that HAM is very powerful and efficient method to find the exact solution of wide range of linear and nonlinear problems References [1] S Abbasbandy Homotopy analysis method for generalized Benjamin-Bona-Mahony equation Z Angew MathPhys 59 (2008): 51-62 [2] S Abbasbandy and F S Zakaria Soliton solutions for the fifth-order K-dVequation with the homotopy analysis method Nonlinear Dyn 51 (2008): 83-87 [3] S Abbasbandy Homotopy analysis method for the Kawahara equation Nonlinear Anal (B) 11 (2010): 307-312
76 [4] S Abbasbandy E Shivanian and K Vajravelu Mathematical properties ofcurve in the framework of the homotopy analysis method Communications in Nonlinear Science and NumericalSimulation16 (2011): 4268-4275 [5] S Abbasbandy The application of homotopy analysis method to nonlinear equations arising in heat transfer Phys Lett A 360 (2006): 109 113 [6] S Abbasbandy The application of homotopy analysis method to solve a generalized Hirota- Satsuma coupled KdV equation Phys Lett A 361 (2007): 478 483 [7] S J Liao Beyond Perturbation: Introduction to the Homotopy Analysis Method CRC Press Boca Raton Chapman and Hall 2003 [8] S J Liao On the homotopy analysis method for nonlinear problems Appl Math Comput 147 (2004) :499 513 [9] S J Liao Comparison between the homotopy analysis method and homotopy perturbation method Appl Math Comput 169 (2005): 1186 1194 [10] M E Gurtin R C Maccamy On the diffusion of biological population Mathematical bioscience 33 (1977): 35-49 [11] W S C Gurney R M Nisbet The regulation of in homogenous populations journal of theoretical biology 52 (1975): 441-457 [12] Y G Lu Holder estimates of solution of biological population equations Applied Mathematics Letters 13(2000): 123-126 [13] Y Tan S Abbasbandy Homotopy analysis method for quadratic Riccati differential Equation Commun Nonlin Sci Numer Simul 13 (2008): 539 546 [14] T Hayat M Khan S Asghar Homotopy analysis of MHD flows of an Oldroyd 8 constant fluid Acta Mech 168 (2004): 213 232 [15] T Hayat M Khan Homotopy solutions for a generalized second-grade fluid past a porous Plate Nonlinear Dyn 42 (2005): 395 405 [16]S T Mohyud-Din and A Yildirim The numerical solution of three dimensional Helmholtz equation Chinese Physics Letters 27 (6) (2010) 060201 [17] A Yildirim and S T Mohyud-Din Analytical approach to space and time fractional Burger s equations Chinese Physics Letters 27 (9) (2010) 090501