AAAAbstract and Applied Analysis1687-04091085-3375Hindawi Publishing Corporation53796310.1155/2009/537963537963Research ArticleOn Approximate Euler Differential EquationsJungSoon-Mo1MinSeungwook2EloePaul1Mathematics SectionCollege of Science and TechnologyHong-Ik UniversityChochiwon 339-701South Koreahongik.ac.kr2Division of Computer ScienceSangmyung University7 Hongji-dongJongno-guSeoul 110-743South Koreasmu.ac.kr20091210200920090907200907092009160920092009Copyright © 2009This 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.

We solve the inhomogeneous Euler differential equations of the form x2y(x)+αxy(x)+βy(x)=m=0amxm and apply this result to the approximation of analytic functions of a special type by the solutions of Euler differential equations.

1. Introduction

The stability problem for functional equations starts from the famous talk of Ulam and the partial solution of Hyers to the Ulam's problem (see [1, 2]). Thereafter, Rassias  attempted to solve the stability problem of the Cauchy additive functional equation in a more general setting.

The stability concept introduced by Rassias’ theorem significantly influenced a number of mathematicians to investigate the stability problems for various functional equations (see  and the references therein).

Assume that X and Y are a topological vector space and a normed space, respectively, and that I is an open subset of X. If for any function f:IY satisfying the differential inequality

an(x)y(n)(x)+an-1(x)y(n-1)(x)++a1(x)y(x)+a0(x)y(x)+h(x)ε for all xI and for some ε0, there exists a solution f0:IY of the differential equation

an(x)y(n)(x)+an-1(x)y(n-1)(x)++a1(x)y(x)+a0(x)y(x)+h(x)=0 such that f(x)-f0(x)K(ε) for any xI, where K(ε) is an expression of ε only, then we say that the above differential equation satisfies the Hyers-Ulam stability (or the local Hyers-Ulam stability if the domain I is not the whole space X). We may apply these terminologies for other differential equations. For more detailed definition of the Hyers-Ulam stability, refer to [1, 3, 5, 6, 811].

Obloza seems to be the first author who has investigated the Hyers-Ulam stability of linear differential equations (see [12, 13]). Here, we will introduce a result of Alsina and Ger (see ): if a differentiable function f:I is a solution of the differential inequality |y(x)-y(x)|ε, where I is an open subinterval of , then there exists a solution f0:I of the differential equation y(x)=y(x) such that |f(x)-f0(x)|3ε for any xI.

This result of Alsina and Ger has been generalized by Takahasi, Miura, and Miyajima: they proved in  that the Hyers-Ulam stability holds for the Banach space-valued differential equation y(x)=λy(x) (see also ).

Using the conventional power series method, the first author investigated the general solution of the inhomogeneous Hermite differential equation of the form

y′′(x)-2xy(x)+2λy(x)=m=0amxm under some specific condition, where λ is a real number and the convergence radius of the power series is positive. This result was applied to prove that every analytic function can be approximated in a neighborhood of 0 by a Hermite function with an error bound expressed by Cx2ex2 (see ).

In Section 2 of this paper, using power series method, we will investigate the general solution of the inhomogeneous Euler (or Cauchy) differential equation

x2y′′(x)+αxy(x)+βy(x)=m=0amxm, where α and β are fixed complex numbers and the coefficients am of the power series are given such that the radius of convergence is ρ>0. Moreover, using the idea from , we will approximate some analytic functions by the solutions of Euler differential equations.

In this paper, 0 denotes the set of all nonnegative integers.

2. General Solution of Inhomogeneous Euler Equations

The second-order Euler differential equation

x2y′′(x)+αxy(x)+βy(x)=0, which is sometimes called the second-order Cauchy differential equation, is one of the most famous differential equations and frequently appears in applications.

m2+(α-1)m+β=0 is called the auxiliary equation of the Euler differential equation (2.1), and every solution of (2.1) is of the form

yh(x)={c1xm1+c2xm2if  m1  and  m2  are  distinct  roots  of  (2.2),(c1+c2lnx)x(1-α)/2if  (1-α)2  is  a  double  root  of  (2.2), where c1 and c2 are complex constants (see [21, Section 2.7]).

Theorem 2.1.

Let α and β be complex constants such that no root of the auxiliary equation (2.2) is a nonnegative integer. If the radius of convergence of power series m=0amxm is at least ρ>0, then every solution y:(0,ρ) of the inhomogeneous Euler differential equation (1.4) can be expressed by y(x)=yh(x)+m=0amxmm2+(α-1)m+β for all x(0,ρ), where yh(x) is a solution of the Euler differential equation (2.1).

Proof.

Assume that y:(0,ρ) is a function given by (2.4), where yh(x) is a solution of the homogeneous Euler differential equation (2.1). We first prove that the function yp(x), defined by y(x)-yh(x), satisfies the inhomogeneous equation (1.4).

Indeed, we have x2yp(x)+αxyp(x)+βyp(x)=m=2m(m-1)amxmm2+(α-1)m+β+m=1αmamxmm2+(α-1)m+β+m=0βamxmm2+(α-1)m+β=m=0amxm, which proves that yp(x) is a particular solution of the inhomogeneous equation (1.4). Moreover, each power series appearing in the above equalities has the same radius of convergence as m=0amxm (which can be verified by using the ratio test). Since every solution to (1.4) can be expressed as a sum of a solution yh(x) of the homogeneous equation and a particular solution yp(x) of the inhomogeneous equation, every solution of (1.4) is certainly of the form (2.4).

We will now apply the ratio test to the power series in (2.4). Indeed, by setting cm=am/[m2+(α-1)m+β], we get |cm+1cm|=|m2+(α-1)m+β(m+1)2+(α-1)(m+1)+β||am+1am|, and hence we have limsupm|cm+1cm|=limsupm|am+1am|, which implies that the power series given in (2.4) has the same radius of convergence as power series m=0amxm, which is at least ρ. That is, y(x) given in (2.4) is well defined on its domain (0,ρ).

3. Approximate Euler Differential Equations

In this section, assume that α and β are complex constants and ρ is a positive constant. For a given K0, we denote by 𝒞K the set of all functions y:(0,ρ) with the properties (a) and (b):

y(x) is expressible by a power series m=0bmxm whose radius of convergence is at least ρ;

m=0|amxm|K|m=0amxm| for any x(0,ρ), where we set am=[m2+(α-1)m+β]bm for m0.

Let {bm} be a sequence of positive real numbers such that the radius of convergence of the series m=0bmxm is at least ρ, and let α and β satisfy either α1 and β>0 or β(α-1)2/4. If a function y:(0,ρ) is defined by y(x)=m=0bmxm, then y certainly belongs to 𝒞K with K1. So, the set 𝒞K is not empty if K1. In particular, if ρ is small and K is large, then 𝒞K is a large class of analytic functions y:(0,ρ).

We will now solve the approximate Euler differential equations in a special class of analytic functions, 𝒞K.

Theorem 3.1.

Let α and β be complex constants such that no root of the auxiliary equation (2.2) is a nonnegative integer. If a function y𝒞K satisfies the differential inequality |x2y(x)+αxy(x)+βy(x)|ε for all x(0,ρ) and for some ε0, then there exists a solution yh:(0,ρ) of the Euler differential equation (2.1) such that |y(x)-yh(x)|(supk01|k2+(α-1)k+β|)Kρρ-xε for any x(0,ρ).

Proof.

Since y belongs to 𝒞K, it follows from (a) and (b) that x2y′′(x)+αxy(x)+βy(x)=m=0[m2+(α-1)m+β]bmxm=m=0amxm for all x(0,ρ). By considering (3.1) and (3.3), we get |m=0amxm|ε for any x(0,ρ). This inequality, together with (b), yields that m=0|amxm|K|m=0amxm|Kε for each x(0,ρ).

Now, suppose that an arbitrary x(0,ρ) is given. Then we can choose an arbitrary constant ρ0(x,ρ). By Abel's formula (see [22, Theorem 6.30] ), we have m=0n|amρ0m||x/ρ0|m|m2+(α-1)m+β|=(m=0n|amρ0m|)|x/ρ0|n+1|(n+1)2+(α-1)(n+1)+β|+k=0n(m=0k|amρ0m|)|x/ρ0|k|k2+(α-1)k+β|M(k;ρ0;x), where we set M(k;ρ0;x)=1-|k2+(α-1)k+β(k+1)2+(α-1)(k+1)+β||xρ0|. Since M(k;ρ0;x)1 for any k0 and ρ0(x,ρ), it follows from (3.5) and (3.6) that m=0n|amρ0m||x/ρ0|m|m2+(α-1)m+β|Kε(|x/ρ0|n+1|(n+1)2+(α-1)(n+1)+β|+k=0n|x/ρ0|k|k2+(α-1)k+β|) If we let n in the above inequality, then we obtain m=0|amxm||m2+(α-1)m+β|Kεk=01|k2+(α-1)k+β||xρ0|k(supk01|k2+(α-1)k+β|)Kρ0ρ0-|x|ε for all x(0,ρ) and for any ρ0(x,ρ). Since ρ0/(ρ0-x)ρ/(ρ-x) as ρ0ρ, we get m=0|amxm||m2+(α-1)m+β|(supk01|k2+(α-1)k+β|)Kρρ-xε for every x(0,ρ).

Finally, it follows from (3.3), Theorem 2.1, and (3.10) that there exists a solution yh of the Euler differential equation (2.1) such that |y(x)-yh(x)||m=0amxmm2+(α-1)m+β|m=0|amxm||m2+(α-1)m+β|(supk01|k2+(α-1)k+β|)Kρρ-xε for any x(0,ρ).

4. An Example

We fix α=0, β=1/4 and suppose that a small ε>0 is given. We can choose a constant 0<ρ<1 such that

ε=m=1(m-1+14m)ρm.

We will consider the function y(x)=ln(1-x) which can be expressed by the power series m=1(-xm/m), whose radius of convergence is 1.

If we set b0=0 and bm=-1/m for any m, then it follows from (b) that

am={0for  m=0,-m+1-14mfor  m. Thus, we have

m=0|amxm|=m=1(m-1+14m)xm=|m=0amxm| for any x(0,ρ), which enables us to choose K=1. Thus, it holds that yC1.

Moreover, we have

|x2y′′(x)+14y(x)|=m=1(m-1+14m)xmm=1(m-1+14m)ρm=ε for all x(0,ρ). It now follows from Theorem 3.1 that there exist complex numbers A and B with

|ln(1-x)-(A+Blnx)x  |4ρρ-xm=1(m-1+14m)ρm for any x(0,ρ).

Acknowledgments

The authors would like to express their cordial thanks to the referees for their useful remarks and constructive comments which have improved the first version of this paper. This work was supported by National Research Foundation of Korea Grant funded by the Korean Government (no. 2009-0071206).

HyersD. H.On the stability of the linear functional equationProceedings of the National Academy of Sciences of the United States of America194127222224MR000407610.1073/pnas.27.4.222ZBL0061.26403UlamS. M.A Collection of Mathematical Problems1960New York, NY, USAInterscienceMR0120127RassiasTh. M.On the stability of the linear mapping in Banach spacesProceedings of the American Mathematical Society1978722297300MR50732710.2307/2042795ZBL0398.47040AokiT.On the stability of the linear transformation in Banach spacesJournal of the Mathematical Society of Japan195026466MR0040580ZBL0040.35501CzerwikS.Functional Equations and Inequalities in Several Variables2002SingaporeWorld ScientificMR1904790FortiG. L.Hyers-Ulam stability of functional equations in several variablesAequationes Mathematicae1995501-2143190MR133686610.1007/BF01831117ZBL0836.39007GăvruţaP.A generalization of the Hyers-Ulam-Rassias stability of approximately additive mappingsJournal of Mathematical Analysis and Applications19941843431436MR128151810.1006/jmaa.1994.1211ZBL0818.46043HyersD. H.IsacG.RassiasTh. M.Stability of Functional Equations in Several Variables1998Boston, Mass, USABirkhäuserMR1639801HyersD. H.RassiasTh. M.Approximate homomorphismsAequationes Mathematicae1992442-3125153MR118126410.1007/BF01830975ZBL0806.47056RassiasTh. M.On the stability of functional equations and a problem of UlamActa Applicandae Mathematicae200062123130MR177801610.1023/A:1006499223572ZBL0981.39014JungS.-M.Hyers-Ulam-Rassias Stability of Functional Equations in Mathematical Analysis2001Palm Harbor, Fla, USAHadronic PressMR1841182ObłozaM.Hyers stability of the linear differential equationRocznik Naukowo-Dydaktyczny. Prace Matematyczne199313259270MR1321558ZBL0964.34514ObłozaM.Connections between Hyers and Lyapunov stability of the ordinary differential equationsRocznik Naukowo-Dydaktyczny. Prace Matematyczne199714141146MR1489723ZBL1159.34332AlsinaC.GerR.On some inequalities and stability results related to the exponential functionJournal of Inequalities and Applications199824373380MR167190910.1155/S102558349800023XZBL0918.39009TakahasiS.-E.MiuraT.MiyajimaS.On the Hyers-Ulam stability of the Banach space-valued differential equation y=λyBulletin of the Korean Mathematical Society2002392309315MR1904665ZBL1011.34046MiuraT.JungS.-M.TakahasiS.-E.Hyers-Ulam-Rassias stability of the Banach space valued linear differential equations y=λyJournal of the Korean Mathematical Society20044169951005MR2095549ZBL1069.34079JungS.-M.Legendre's differential equation and its Hyers-Ulam stabilityAbstract and Applied Analysis20072007145641910.1155/2007/56419MR2365813ZBL1153.34306JungS.-M.Approximation of analytic functions by Hermite functionsBulletin des Sciences Mathematiques. In press2-s2.0-3664902524210.1016/j.bulsci.2007.11.001JungS.-M.Approximation of analytic functions by airy functionsIntegral Transforms and Special Functions20081912885891MR246709810.1080/10652460802321287ZBL1155.34002JungS.-M.An approximation property of exponential functionsActa Mathematica Hungarica20091241-2155163MR252062410.1007/s10474-009-8167-1KreyszigE.Advanced Engineering Mathematics19794thNew York, NY, USAJohn Wiley & SonsMR520566WadeW. R.An Introduction to Analysis20002ndUpper Saddle River, NJ, USAPrentice-Hall