@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 On Solution of Regular Singular Ordinary Boundary Value Problem Luma. N. M. Tawfiq Heba. W. Rasheed Dept. of Mathematics /College of Education for Pure Science (Ibn Al-Haitham)/ University of Baghdad Received in: 20 October 2011 , Accepted in: 7 Desember 2011 Abstract This paper devoted to the analysis of regular singular boundary value problems for ordinary differential equations with a singularity of the different kind , we propose semi - analytic technique using two point osculatory interpolation to construct polynomial solution, and discussion behavior of the solution in the neighborhood of the regular singular points and its numerical approximation. Many examples are presented to demonstrate the applicability and efficiency of the methods. Finally , we discuss behavior of the solution in the neighborhood of the singularity point which appears to perform satisfactorily for singular problems. Kay ward : Singular boundary value problems, ODE, BVP. 249 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 Introduction In the study of nonlinear phenomena in physics, engineering and other sciences, many mathematical models lead to singular two-point boundary value problems (SBVP) associated with nonlinear second order ordinary differential equations (ODE) . In mathematics, a singularity is in general a point at which a given mathematical object is not defined, or a point of an exceptional set where it fails to be well-behaved in some particular way, such as Many problems in varied fields as thermodynamics, electrostatics, physics, and statistics give rise to ordinary differential equations of the form : y″ = f(t, y, y') , 0 < x < 1 , (1) On some interval of the real line with some boundary conditions. A two-point BVP associated to the second order differential equation (1) is singular if one of the following situations occurs : a and/or b are infinite; f is unbounded at some x0∈ [0,1] or f is unbounded at some particular value of y or y′ [1] . How to solve a linear ODE of the form : A(x)y′′ + B(x)y′ + C(x)y = 0 , (2) The first thing we do is, rewrite the ODE as : y′′ + P(x)y′ + Q(x)y = 0 , (3) where, of course, P(x) = B(x) / A(x) , and Q(x) = C(x) / A(x) . there are two types of a point x0 ∈ [0,1] : Ordinary Point and Singular Point. Also, there are two types of Singular Point : Regular and Irregular Points, A function y(x) is analytic at x0 if it has a power series expansion at x0 that converges to y(x) on an open interval containing x0 .A point x0 is an ordinary point of the ODE (3), if the functions P(x) and Q(x) are analytic at x0. Otherwise x0 is a singular point of the ODE, i.e. P(x) = P0 + P1(x-x0) + P2(x-x0)2+…….. = i i i xxp )( 0 0 −∑ ∞ = , (4) Q(x) = Q0 + Q1(x-x0) + Q2(x-x0)2+………. = i i i xxq )( 0 0 −∑ ∞ = , (5) If A , B and C are polynomials then a point x0 such that A(x0) ≠ 0 is an ordinary point. On the other hand if P(x) or Q(x) are not analytic at x0 then x0 is said to be a singular. A singular point x0 of the ODE (3) is a regular singular point of the ODE if the functions xP(x) and x2Q(x) are analytic at x0. Otherwise x0 is an irregular singular point of the ODE . [2] L.F. Shampine in [3] give other definition, which illustrated by the following : If limx→x0 (x – x0)P(x) finite and limx→x0 (x – x0)2 Q(x) finite , (6) that is, if both (x –x0)P(x) and (x –x0)2Q(x) posses a Taylor series at x0, then x0 is called a regular singular point, otherwise x0 is an irregular singular point . If A, B and C are polynomials and suppose A(x0) = 0, then x0 is a regular singular point if : limx → x0 (x - x0)( B/ A) and limx → x0 (x - x0)2(C/A) are finite , (7) Now, we state the following theorem without proof which gives us a useful way of testing if a singular point is regular. Theorem 1 [4] If the limx→0 P(x) and limx→0 Q(x) exist, are finite, and are not 0 then x = 0 is a regular singular point. If both limits are 0, then x = 0 may be a regular singular point or an ordinary point. If either limit fails to exists or is ±∞ then x = 0 is an irregular singular point . There are four kinds of singularities : • The first kind is the singularity at one of the ends of the interval [0,1] ; 250 | Mathematics http://en.wikipedia.org/wiki/Mathematics http://en.wikipedia.org/wiki/Set_(mathematics) http://en.wikipedia.org/wiki/Well-behaved @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 • The second kind is the singularity at both ends of the interval [0,1] • The third kind is the case of a singularity in the interior of the interval; • The forth and final kind is simply treating the case of a regular differential equation on an infinite interval. In this paper, we focus of the first three kinds . 2. Solution of Second Order SBVP In this section we suggest semi analytic technique to solve second order SBVP as following, we consider the SBVP : xm y''+ f( x, y, y' ) = 0 , (8a) gi( y(0) , y(1) ,y'(0), y'(1) ) = 0 , i = 1 , 2 , (8b) where f , g1, g2 are in general nonlinear functions of their arguments . The simple idea behind the use of two-point polynomials is to replace y(x) in problem (8), or an alternative formulation of it, by a P2n+1 which enables any unknown boundary values or derivatives of y(x) to be computed . The first step therefore is to construct the P2n+1 , to do this we need the Taylor coefficients of y (x) at x = 0 : y = a 0 + a 1 x + ∑ ∞ =2i a i x i , (9) where y(0)= a0 ,y'(0)= a1 ,y"(0) / 2! =a2 ,…, y(i)(0) / i! = ai , i= 3, 4,… then insert the series forms (9) into (8a) and equate coefficients of powers of x to obtain a2 . Also we need Taylor coefficient of y(x) about x = 1 : y = b 0 + b 1 (x-1) + ∑ ∞ =2i b i (x-1) i , (10) where y(1) = b0 , y'(1) =b1 , y"(1) / 2! =b2,…, y(i)(1) / i! =bi ,i = 3,4,… then insert the series form (10) into (8a) and equate coefficients of powers of (x-1) to obtain b2 ,then derive equation (8a) with respect to x and iterate the above process to obtain a3 and b3 , now iterate the above process many times to obtain a4 ,b4 ,then a5 ,b5 and so on, that is ,we can get ai and bi for all i ≥ 2( the resulting equations can be solved using MATLAB to obtain ai and bi for all i ≥ 2 ) , the notation implies that the coefficients depend only on the indicated unknowns a0, a1, b0, b1, we get two of these four unknown by the boundary condition .Now, we can construct a P2n+1(x) from these coefficients ( aisۥ and bisۥ ) by the following : P2n+1 = ∑ = n i 0 { ai Qi(x) + (-1)i bi Qi(1-x) } , (11) where ( x j / j!)(1-x) 1+n ∑ − = jn s 0       + s sn xs = Q j (x) / j! we see that (11) have only two unknowns from a0 ,b0 ,a1 and b1 to find this, we integrate equation (8a) on [0, x] to obtain : xmy'(x) – mxm-1y(x) + m(m–1) ∫ x 0 xm-2y(x) dx + ∫ x 0 f(x, y, y') dx = 0 , (12a) and again integrate equation (12a) on [0, x] to obtain : xmy(x) –2m ∫ x 0 xm-1y(x) dx +m(m-1)∫ x 0 (1-x)xm-2y(x)dx+∫ x 0 (1-x)f(x,y,y') = 0 , (12b) Putting x = 1 in (12) then gives : b1 – mb0 + m(m-1) ∫ 1 0 xm-2 y(x) dx + ∫ 1 0 f(x, y, y') dx = 0 , (13a) and 251 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 b0–2m∫ 1 0 xm-1y(x)dx +m(m-1)∫ 1 0 (1-x)xm-2 y(x) dx + ∫ 1 0 (1-x)f(x, y, y')dx = 0 , (13b) Use P2n+1 as a replacement of y(x) in ( 13 ) and substitute the boundary conditions (8b) in (13) then, we have only two unknown coefficients b1, b0 and two equations (13) so, we can find b1, b0 for any n by solving this system of algebraic equations using MATLAB, so insert b0 and b1 into (11) , thus (11) represents the solution of (8) . Extensive computations have shown that this generally provides a more accurate polynomial representation for a given n . 3. Examples In this section, many examples will be given to illustrate the efficiency, accuracy , implementation and utility of the suggested method. The bvp4c solver of MATLAB has been modified accordingly so that it can solve some class of SBVP as effectively as it previously solved nonsingular BVP. Example 1 Consider the following SBVP : x2 y"(x) – 9xy'(x) + 25y(x) = 0 , 0 ≤ x ≤ 1 With BC : y'(0) = 0 , y(1) = 1. The exact solution is : y(x) = x5 It is clear that x = 0, is regular singular point and it is singularity of first kind . Now, we solve this example using semi - analytic technique , From equations (11) we have : P9(x) = x5 For more details ,table(1) give the results for different nodes in the domain, for n = 4, i.e. P9 and errors obtained by comparing it with the exact solution. Figure (1) gives the accuracy of the suggested method. Example 2 Consider the following SBVP: (1- x2) y" + x y' + y = 0 , 0 ≤ x ≤ 1 With BC : y'(0) = 0 , y'(1) = - y(1) . It is clear that x = 1, is regular singular point and it is singularity of first kind . Now, we solve this example using semi-analytic technique ,From equation (11) we have : if n = 2 ,we have P5 as follows : P5 = - 0.0109890110 x5 + 0.0627943485 x4 - 0.3296703297x3 - 0.9293563579x2 + x + 1.858712715855573. Higher accuracy can be obtained by evaluating higher n, now take n = 3, i.e., P7 = - 0.0017069701x7 +0.0098130167x6 - 0.0293029872 x5 +0.0766016098x4 - 0.3333333333 x3 - 0.9288921481 x2 + x + 1.8577842962 . Now, increase n, to get higher accuracy , let n = 4, i.e. , P9 = -0.0003876355x9 +0.0025102583x8 -0.007761001x7 +0.0167512244x6 - 0.033114847889001 x5 + 0.077409158926178 x4 - 0.3333333333x3 – 0.9289099071x2 + x + 1.8578198142 . For more details ,table (2) gives the results of different nodes in the domain, for n = 2, 3, 4. Also, figure (2) illustrate suggested method for n= 4. Example3 Consider the following SBVP: – x2 y"– 2xy' + 2y = – 4x2 , 0 ≤ x ≤ 1 With BC : y(0) = y(1) = 0 and exact solution is y = x2 – x It is clear that x = 0, is regular singular point and it is singularity of first kind . Now, we solve this example using semi-analytic technique , From equation (11) we have : P9 = x2 – x . For more details ,table (3) give the results of different nodes in the domain, for n = 4. Also, figure (3) illustrate suggested method for n= 4. 4. Behavior of the solution in the neighborhood of the singularity x= 0 252 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 Our main concern in this section will be the study of the behavior of the solution in the neighborhood of singular point . Consider the following SIVP : y′′(x) + ((N − 1) / x ) y′(x) = f(y) , N ≥ 1 , 0 < x < 1 , (14) y(0) = y0 , limx→0+ x y′(x) = 0 , (15) where f(y) is continuous function . As the same manner in [6], let us look for a solution of this problem in the form : y(x) = y0 − C xk (1 + o(1)) , (16) y′(x) = − C k xk−1(1 + o(1)) , y′′(x) = − C k (k − 1) xk−2(1 + o(1)) , x → 0+ where C is a positive constant and k > 1. If we substitute (16) in (14) we obtain : C = (1/ k) (f(y0) /N )k−1 , (17) In order to improve representation (16) we perform the variable substitution : y(x) = y0 − C xk (1 + g(x)) , (18) we easily obtain the following result which is similar to the results in [6]. Theorem 2 For each y0 > 0, problem (14), (15) has, in the neighborhood of x = 0, a unique solution that can be represented by : y(x, y0) = y0 − C xk (1 + g xk + o(xk) ) , where k, C and g are given by (17) and (18), respectively. We see that these results are in good agreement with the ones obtained by the method in [6], they are also consistent with the results presented in [7]. In order to estimate the convergence order of the suggested method at x = 0, we have carried out several experiments with different values of n and used the formula : cy0 = − log2 ( |y0n3 − y0n2| / |y0n2− y0n1 | ) , (19) where y0ni is the approximate value of y0 obtained with ni ,ni = 1,2, 3, 4,… References 1. Robert ,L.B. and S. C. Courtney , (1996) " Differential Equations A Modeling perspective "United States of America ,. .2. Rachůnková,I. ; Staněk ,S. ,and Tvrdý , M. (2008) " Solvability of Nonlinear Singular Problems for Ordinary Differential Equations ",New York, USA, . 3. Shampine ,L.; F. Kierzenka, J.; and Reichelt ,M. W. (2000)" Solving Boundary Value Problems for Ordinary Differential Equations in Matlab with bvp4c", 4. Howell , K.B.(2009) " Ordinary Differential Equations" , USA , Spring 5. Burden ,L. R. and Faires , J. D. (2001). " Numerical Analysis " , Seventh Edition. 6. Morgado L. and Lima ,P. (2009) " Numerical methods for a singular boundary value problem with application to a heat conduction model in the human head", Proceedings of the International Conference on Computational and Mathematical Methods in Science and Engineering, CMMSE . 7. Abukhaled, M.; Khuri,S.A. and Sayfy , A. (2011) A NUMERICAL APPROACH FOR SOLVING A CLASS OF SINGULAR BOUNDARY VALUE PROBLEMS ARISING IN PHYSIOLOGY ", INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING , 8, No.2,:353–363, 253 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 Table (1): The result of the method for P9 of example 1 P9(x) 0 a0 5 b1 Error | y(x) – P9| Osculatory interpolation P9 exact solution y(x) 0 0 0 0 0 0.000010000000000 0.000010000000000 0.1 0 0.000320000000000 0.000320000000000 2 0 0 0.002430000000000 0.002430000000000 .3 0 0 0.010240000000000 0.010240000000000 .4 0 0 0.031250000000000 0.031250000000000 .5 0 0 0.077760000000000 0.077760000000000 .60 0 0.168070000000000 0.168070000000000 .70 0 0.327680000000000 0.327680000000000 .80 0 0.590490000000000 0.590490000000000 .90 0 1 1 1 Table (2): The result of the method for n = 2, 3, 4 of example 2 P9 P7 P5 1.857819814228270 1.857784296228232 1.858712715855573 a0 1.650983731093794 1.650963483915429 1.65149136577708 b0 P9 P7 P5 xi 1.857819814228270 1.857784296228232 1.858712715855573 0 1.948204807591012 1.948169418187120 1.949095651491374 0.1 2.018110988145184 2.018075735440336 2.018998053375218 0.2 2.065775139074470 2.065740050052486 2.066651475667217 0.3 2.089560915079904 2.089526555233489 2.090411805337528 0.4 2.087939057875687 2.087906467474147 2.088746075353169 0.5 2.059472007087712 2.059442143114885 2.060219277864830 0.6 2.002801634885403 2.002774819034724 2.003481177393688 0.7 1.916639282152951 1.916615117147540 1.917253124018220 0.8 1.799757582816419 1.799735452392664 1.800314866561014 0.9 1.650983731088339 1.650963483906875 1.651491365775585 1 Table( 3): The result of the method for n = 4 of example 3 P9 1- a1 1 b1 |9y-P| P9 y: exact 0 0.210000000000000- 0.210000000000000- .30 0 .2400000000000000- .2400000000000000- .60 0 .0900000000000000- .0900000000000000- .90 254 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 Table( 4): Comparison between suggested and other method given in [5]of example 3 Errors |y(x) –P9| P9(x) using Osculatory interpolation y2(x) using cubic spline y1(x) using piecewise linear algorithm exact solution y(x) xi 0 .2100000000000- .2100000000000- .2123333333330- 0.210000000000- .30 0 .2400000000000- .2400000000000- .2413333333330- .2400000000000- 0.6 0 .0900000000000- 0.090000000000- .0903333333330- .0900000000000- 0.9 Fig.( 1):Comparison between the exact and semi-analytic solution P9 of example1 Fig.(2): illustrate suggested method for n= 4,i.e.,P9 of example 2 . Fig. (3): illustrate suggested method for n= 4,i.e.,P9 of example 3 . 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 the solution at n=4 x y exact p9 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 2.05 2.1 the solution at n=4 x y p9 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -0.25 -0.2 -0.15 -0.1 -0.05 0 the solution at n=4 x y exact P9 255 | Mathematics @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ÚÓ‘Ój�n€a@Î@Úœäñ€a@‚Ï‹»‹€@·rÓ:a@Âig@Ú‹©@Ü‹1a26@@ÖÜ»€a@I2@‚b«@H2013 Ibn Al-Haitham Jour. for Pure & Appl. Sci. Vol. 26 (2) 2013 النظامیة الشاذة االعتیادیةمسائل القیم الحدودیة حول حل لمى ناجي محمد توفیق ھبة ولید رشید جامعة بغداد / )ابن الھیثم( للعلوم الصرفة كلیة التربیة / الریاضیات علوم قسم 2011كانون الثاني 7قبل البحث في 2011الثاني تشرین 20استلم البحث في : الخالصة للمعادالت التفاضلیة االعتیادیة الھدف من ھذا البحث عرض دراسة تحلیلیة لمسائل القیم الحدودیة النظامیة الشاذة بوصفة النقطتین للحصول على الحل يالتقنیة شبھ التحلیلیة باستخدام االندراج التماسي ذأننا نقترح اذوبأنواع مختلفة وسھولة أداء الطریقة المقترحة و أخیرا ناقشنا سلوك الكافیة ،حدود، كذلك ناقشنا عدد من األمثلة لتوضیح الدقة ، متعددة و اقترحنا صیغة جدیدة مطورة لتخمین الخطأ تساعد في تقلیل .الحل في جوار النقاط الشاذة و إیجاد الحل التقریبي لھا لحسابات العملیة وإظھار النتائج بشكل مرضي فیما یخص المسائل الشاذة .ا عتیادیة ، مسائل القیم الحدودیة ه ، معادالت تفاضیلھ امسائل القیم الحدودیة الشاذ الكلمات المفتاحیة : 256 | Mathematics Fig. (3): illustrate suggested method for n= 4,i.e.,P9 of example 3 .