A Numerical-Computational Technique for Solving. Transformed Cauchy-Euler Equidimensional. Equations of Homogeneous Type

Transcription

1 Advanced Studies in Theoretical Physics Vol. 9, 015, no., 85-9 HIKARI Ltd, A Numerical-Computational Technique for Solving Transformed Cauchy-Euler Equidimensional Equations of Homogeneous Type S. O. Edeki *, A. A. Opanuga, H. I. Okagbue and G. O. Akinlabi 1 Department of Mathematical Sciences Covenant University, Canaanland, Otta, Nigeria * Corresponding author S. A. Adeosun Department of Mathematical Sciences Crescent University, Abeokuta, Nigeria A. S. Osheku Department of Mathematics University of Ibadan, Ibadan, Nigeria Copyright 015 S. O. Edeki, A. A. Opanuga, H. I. Okagbue, G. O. Akinlabi, S. A. Adeosun, and A. S. Osheku. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract This work considers the solution of Transformed-Cauchy-Euler differential equations via Differential Transform method (DTM). For illustration and application of the method s efficiency and reliability, two examples of order and 3 are solved. The results agreed with the exact solution obtained via Laplace transform method. Keywords: Cauchy-Euler differential equations, DTM, Laplace Transform method, Series solution

2 86 S. O. Edeki et al. 1. Introduction A Cauchy-Euler equation, also known as Euler s equation named after Augustine Louis Cauchy ( ), is a linear homogeneous ordinary differential equation with variable coefficients. This type of equation commonly appears as second order equation in physics and engineering applications such as Laplace s equation in polar coordinates and many other aspects of applied sciences [1]. One of the characteristic features of Cauchy-Euler equation is that the order of each derivative in the equation equals the power of the independent variable, say, x in the corresponding coefficient. This type of equations can be reduced to linear equation with constant coefficients by using appropriate substitution [1-]. Various numerical methods have been proposed to solve ordinary differential equations. The Differential transformation method (DTM) adopted in this work was first proposed by Zhou [3] to solve both linear and nonlinear initial value problems in electric circuits analysis. Other applications of Differential Transform method can be found in [4-10]. It is also very useful when dealing with integro-differential equations encountered in finance [11]. In this work, a numerical-computational method is applied to solve Cauchy-Euler Equidimensional Equations. To illustrate this, second and third order Cauchy-Euler differential equations of the transformed type are solved.. The Cauchy-Euler Equidimensional Equation A Cauchy-Euler differential equation is a differential equation of the form: a x y a x y... a xy a y 0 n n ( n) n1 ( n1) 1 n1 1 (1) ( m) where a 0, a 1,..., a n are constants and y represents an mth derivative of the unknown function of yx ( ). Considering a second order differential equation of the form: ( ) ( ) 0 () x x ax t b The solution can be obtained through change of variables by applying Euler s formula: x e t, ( t ln x & xdt dx) (3) 1 ( x) x ( t) and (4) ( x) x ( ( t) ( t)) (5)

3 Numerical-computational technique 87 substituting the results in (3-5) in () yields a second order differential equation with constant coefficients of the form: ( ( ) ( )) ( ) 0 (6) 1 x x t t axx t b This on simplification gives: ( t) ( t) a ( t) b 0 (7) ( t) ( a 1) ( t) b 0 (8) Equations (7-8) above can be solved using the characteristic polynomial: ( a1) b 0 (9) With roots 1 and which give the general solution depending on the type of roots it has (or w.r.t. (3) resp.) 1 1 Case1: ( ) t t t a e a e (or ( x) a x a x ) (10) Case: ( ) t t t a e a te ( or ( x) a x a (ln x) x ) (11) t ( t) e ( a1 cos t a sin t) 1 or ( x) x ( a1 cos( ln x) a sin( ln x)) Case (3) (1) 3. Analysis of Differential Transformation Method In illustrating this method, we consider an arbitrary one dimensional function wxwhose ( ) differential transform is defined as: k 1 d w( x) Wk ( ) k k! dx x 0 (13) In the equation above, wx ( ) is the original function and Wk ( ) is the transformed function. The Taylor series about a point x 0 is given as k k x d w( x) wk ( ) k! dx (14) k k0 x0 The inverse differential transform is expressed as:

4 88 S. O. Edeki et al. w( x) W ( k) x k (15) k0 The following theorems are deduced from (13) (15). The proofs of these are found in [8] and in other related research papers. i. If w( x) v1( x) v( x), then W( k) V1( k) V( k) ii. If w( x) v( x), then W( k) V ( k), where is a constant ( ) iii. If w( x) dv x, then W( k) ( k 1) V ( k 1) dx ( ) iv. If w( x) d v x, then W( k) ( k 1)( k ) V ( k ) dx n If w( x) x, then W( k) ( k n), where is the kronecker delta. 4. Numerical Results 1, for kn 0, for kn In this subsection, we apply differential transform method to solve some examples of Cauchy-Euler differential equations of the transformed type. Problem 1 []: Consider the following second order Cauchy-Euler differential equation: 3 x y ( x) 4 xy( x) y 0, y(1), y(1) 1 (16) Using the above substitution, (16) is transformed to a second order differential equation with constant coefficients: 3 y ( t) 7 y( t) y 0, y(0), y(0) 1 (17) The theoretical solution of the equation is: 1 9 t 1 3 t t y( t) e e, ( x e ) (18) 5 5 Taking the differential transform of (17), we obtain: 3( k 1)( k ) Y( k ) 7( k 1) Y( k 1) Y( k) (19) Y(0), Y(1) 1 (0) Therefore, the recurrence relation in (19-0) gives the following:

5 Numerical-computational technique 89 1 Y(), for k 0, 5 Y(3), for k 1, 18 9 Y(4), for k Y(5), for k Hence, k t y( t) Y( k) t t t t t t, ( x e ) (1) k0 Table I: Tabular Solutions of problem1 t DTM EXACT ABS ERR Figure1: Solution graphs of problem1 Remark 1: In table I, a numerical comparison between the exact solution and the 5 approximate solution (up to the term containing t via the DTM) along with the corresponding absolute errors (ABS ERR) is presented. Also, graphical look of such solutions is in figure1, with series1 & series for the DTM & the exact solutions respectively. Problem []: Consider the following third order Cauchy-Euler differential equation: ( ) 10 ( ) 0 0 0, (1) 0, (1) 1, (1) 1 3 x y x x y x xy y y y () Equation () is transformed to a second order differential equation with constant coefficient: y ( t) 7 y ( t) 8y 0 0, y(0) 0, y(0) 1, y (0) 1 (3)

6 90 S. O. Edeki et al. The theoretical solution of () is: 3 t 1 10t 8 t t y( t) e e e, ( x e ) (4) The differential transform of (3) is: with ( k 1)( k )( k 3) Y( k 3) [8( k 1) Y( k 1) 0 Y( k) 7( k 1)( k ) Y( k )] 1 Y(0) 0, Y(1) 1, Y() (5) Therefore, the recurrence relation in (5) gives the following: 35 Y(3), for k 0, 6 93 Y(4), for k 1, 4 k y( t) Y( k) t t t t t 6 4 (6) k0 Table II: Tabular Solutions of problem t EXACT DTM ABS ERR E E E E E E E E , , Series1 Series Figure: Solution graphs of problem Remark 1: In table II, a numerical comparison between the exact solution and the 5 approximate solution (up to the term containing t via the DTM) along with the corresponding absolute errors (ABS ERR) is presented. Also, graphical look of such solutions is in figure, with series1 & series for the DTM & the exact solutions respectively.

7 Numerical-computational technique Concluding Remarks In this paper, Differential transform method has been applied successfully to solve transformed Cauchy-Euler differential equations, whose results can be linked easily with the original Cauchy-Euler s differential equations. The results obtained are in strong agreement with the exact solution. Differential transform method is easy to apply to any kind of differential equation. It does not require any form of linearization, perturbation or initial guess, and hence, it is recommended as an efficient and a reliable numerical-computational method for higher order Cauchy-Euler differential equations in applied sciences. Acknowledgements.The authors wish to express sincere thanks to the anonymous reviewers and/or referees for their positive and constructive comments towards the improvement of the paper. References [1] E. Kreyszig, Advanced Engineering Mathematics, Willey, ISBN: , (006). [] Martha L. Abell and James P. Braselton (010): Introductory Differential Equation with Boundary Value Problems. ISBN: [3] J. K. Zhou Differential Transformation and its application for electrical circuits, Huarjung University press, Wuuhahn, China, (1986). [4] Arikoglu A, Ozkol Solution of Difference Equations by using differential transform method, Appl. Math. Comput. 173 (1), (006): [5] S. O. Edeki, A. A. Opanuga and H. I. Okaghue On Iterative Techniques for Numerical Solutions of linear and nonlinear Differential Equations, J. Math. Comput. Sci. 4(4), (014) [6] K. R. Raslan, A. Biswas and Zain F. Abu Sheer Differential Transform method for solving partial differential equations with variable coefficients, International Journal of Physical Sciences. 7(9), (01): [7] A. A. Opanuga, S. O. Edeki, H. I. Okaghue, G. O. Akinlabi, A. S. Osheku and B. Ajayi On Numerical Solution of Systems of Ordinary Differential Equations by Numerical-Analytical method, Applied Mathematical Sciences, 8 (164), (014):

8 9 S. O. Edeki et al. [8] S. O. Edeki, H. I Okaghue, A. A. Opanuga, S. A. Adeosun, A Semi-Analytical Method for Solutions of a Certain Class of Second Order Ordinary Differential Equations, Applied Mathematics, 5, (014): [9] Z. Odibat, Differential Transform method for solving Volterra integral equation with separable kernels, Math. Comput. Modelling 48, (008): [10] J. Biazar, M. Eslami, Differential Transform method for Quadratic Riccati Differential Equation, International Journal of Nonlinear Science, 9(4), (010): [11] S. O. Edeki, I. Adinya, O. O. Ugbebor, The Effect of Stochastic Capital Reserve on Actuarial Risk Analysis via an Integro-differential equation IAENG Int. J. Appl. Math. 44 (), (014): Received: January 10, 015; Published: February 6, 015