JFSAJournal of Function Spaces and Applications0972-68022090-8997Hindawi Publishing Corporation38293210.1155/2012/382932382932Research ArticleSimple Harmonic Oscillator Equation and Its Hyers-Ulam StabilityJungSoon-MoKimByungbaeIsacGeorge1Mathematics SectionCollege of Science and TechnologyHongik University Jochiwon 339-701Republic of Koreahongik.ac.kr2012222012201206032008210420082012Copyright © 2012 Soon-Mo Jung and Byungbae Kim.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.

We solve the inhomogeneous simple harmonic oscillator equation and apply this result to obtain a partial solution to the Hyers-Ulam stability problem for the simple harmonic oscillator equation.

1. Introduction

Let X be a normed space over a scalar field 𝕂 and let I be an open interval, where 𝕂 denotes either or . Assume that a0,a1,,  an:I𝕂 are given continuous functions, g:IX is a given continuous function, and y:IX is an n times continuously differentiable function satisfying the inequality:an(t)y(n)(t)+an-1(t)y(n-1)(t)++a1(t)y(t)+a0(t)y(t)+g(t)ε,

for all tI and for a given ɛ>0. If there exists an n times continuously differentiable function y0:IX satisfyingan(t)y0(n)(t)+an-1(t)y0(n-1)(t)++a1(t)y0(t)+a0(t)y0(t)+g(t)=0,

and y(t)-y0(t)K(ɛ) for any tI, where K(ε) is an expression of ɛ with limɛ0K(ε)=0, then we say that the above differential equation has the Hyers-Ulam stability. For more detailed definitions of the Hyers-Ulam stability, we refer the reader to .

Obloza seems to be the first author who has investigated the Hyers-Ulam stability of linear differential equations (see [8, 9]). Here, we will introduce a result of Alsina and Ger (see ). If a differentiable function f:I satisfies the inequality |y(t)-y(t)|ɛ, where I is an open subinterval of , then there exists a solution f0:I of the differential equation y(t)=y(t) such that |f(t)-f0(t)|3ε for any tI. This result has been generalized by Takahasi et al. Indeed, it was proved in  that the Hyers-Ulam stability holds true for the Banach space valued differential equation y(t)=λy(t) (see also [12, 13]).

Moreover, Miura et al.  investigated the Hyers-Ulam stability of nth order linear differential equation with complex coefficients. They  also proved the Hyers-Ulam stability of linear differential equations of first order, y(t)+g(t)y(t)=0, where g(t) is a continuous function.

Jung also proved the Hyers-Ulam stability of various linear differential equations of first order . Moreover, he could successfully apply the power series method to the study of the Hyers-Ulam stability of Legendre differential equation (see ). Subsequently, the authors  investigated the Hyers-Ulam stability problem for Bessel differential equation by applying the same method.

In Section 2 of this paper, by using the ideas from [20, 21], we investigate the general solution of the inhomogeneous simple harmonic oscillator equation of the form:y(x)+ω2y(x)=m=0amxm,

where ω is a given positive number. Section 3 will be devoted to a partial solution of the Hyers-Ulam stability problem for the simple harmonic oscillator equation (2.1) in a subclass of analytic functions.

2. Inhomogeneous Simple Harmonic Oscillator Equation

A function is called a simple harmonic oscillator function if it satisfies the simple harmonic oscillator equation:y′′(x)+ω2y(x)=0.

The simple harmonic oscillator equation plays a great role in physics and engineering. In particular, it describes quantum particles confined in potential wells in quantum mechanics and the Hyers-Ulam stability of solutions of this equation is very important.

In this section, we define c0=c1=0 and for m1,c2m=i=0m-1(-1)m-i-1a2i(2i)!(2m)!ω2m-2i-2,c2m+1=i=0m-1(-1)m-i-1a2i+1(2i+1)!(2m+1)!ω2m-2i-2, where we refer to (1.3) for the am. We can easily check that these cm satisfy the followingam=(m+2)(m+1)cm+2+ω2cm, for any m{0,1,2,}.

Lemma 2.1.

(a) If the power series m=0amxm converges for all x(-ρ,ρ) with ρ>1, then the power series m=2cmxm with cm given in (2.2) satisfies the inequality |m=2cmxm|C1/(1-|x|) for some positive constant C1 and for any x(-1,1).

(b) If the power series m=0amxm converges for all x(-ρ,ρ) with ρ1, then for any positive ρ0<ρ, the power series m=2cmxm with cm given in (2.2) satisfies the inequality |m=2cmxm|C2 for any x[-ρ0,ρ0] and for some positive constant C2 which depends on ρ0. Since ρ0 is arbitrarily close to ρ, this means that m=2cmxm is convergent for all x(-ρ,ρ).

Proof.

(a) Since the power series m=0amxm is absolutely convergent on its interval of convergence, with x=1, m=0am converges absolutely, that is, m=0|am|<M1 by some number M1.

We know that |a2i|(2i)!(2m)!ω2m-2i-2=|a2i|2m(2m-1)ω(2m-2)ω(2i+1){|a2i|2m(2m-1)(for  0<ω1)|a2i|2m(2m-1)ω[ω](for  ω>1)|a2i|  2m(2m-1)max{1,ωω}, since for ω>1, each factor of the form ω/ in the summand is either less than 1 if >[ω], or is bigger than or equal to 1 if [ω], where [ω] denotes the largest integer less than or equal to ω. Thus, we obtain |c2m|i=0m-1|a2i|(2i)!(2m)!ω2m-2i-2i=0m-1|a2i|2m(2m-1)max{1,ωω}max{1,ωω}i=0|ai|max{M1,M1ωω}C1.

Similarly, we have |a2i+1|(2i+1)!(2m+1)!ω2m-2i-2|a2i+1|2m(2m-1)max{1,ωω},   and |c2m+1|C1 for all m1.

Therefore, we get |m=2cmxm|m=2|cm||xm|C1m=2|xm|C11-|x|, for every x(-1,1).

(b) The power series m=0amxm is absolutely convergent on its interval of convergence, and, therefore, for any given ρ0<ρ, the series m=0|amxm| is convergent on [-ρ0,ρ0] and m=0|am||x|mm=0|am|ρ0mM2 for any x[-ρ0,ρ0]. Now, it follows from (2.2), (2.4), (2.6), and (2.8) that |m=2cmxm|m=1|c2m|ρ02m+m=1|c2m+1|ρ02m+1m=1ρ02mi=0m-1|a2i|(2i)!(2m)!ω2m-2i-2+m=1ρ02m+1i=0m-1|a2i+1|(2i+1)!(2m+1)!ω2m-2i-2m=1i=0m-1|a2i|ρ02imax{1,ωω}2m(2m-1)+  m=1i=0m-1|a2i+1|ρ02i+1max{1,ωω}2m(2m-1)m=1M2max{1,ωω}2m(2m-1)+m=1M2max{1,ωω}2m(2m-1)max{M2,M2ωω}m=11m(2m-1)  32max{M2,M2ωω}C2, for any x[-ρ0,ρ0].

Lemma 2.2.

Suppose that the power series m=0amxm converges for all x(-ρ,ρ) with some positive ρ. Let ρ1=min{1,ρ}. Then, the power series m=2cmxm with cm given in (2.2) is convergent for all x(-ρ1,ρ1). Further, for any positive ρ0<ρ1, |m=2cmxm|C for any x[-ρ0,ρ0] and for some positive constant C which depends on ρ0.

Proof.

The first statement follows from the latter statement. Therefore, let us prove the latter statement. If ρ1, then ρ1=ρ. By Lemma 2.1(b), for any positive ρ0<ρ=ρ1, |m=2cmxm|C2 for each x[-ρ0,ρ0] and for some positive constant C2 which depends on ρ0.

If ρ>1, then by Lemma 2.1(a), for any positive ρ0<1=ρ1, we get |m=2cmxm|C11-|x|C11-ρ0max{C11-ρ0,C2}C, for all x[-ρ0,ρ0] and for some positive constant C which depends on ρ0.

Using these definitions and the lemmas above, we will now show that m=2cmxm is a particular solution of the inhomogeneous simple harmonic oscillator equation (1.3).

Theorem 2.3.

Assume that ω is a given positive number and the radius of convergence of the power series m=0amxm is ρ>0. Let ρ1=min{1,ρ}. Then, every solution y:(-ρ1,ρ1) of the simple harmonic oscillator equation (1.3) can be expressed by y(x)=yh(x)+m=2cmxm,where yh(x) is a simple harmonic oscillator function and cm are given by (2.2).

Proof.

We show that m=2cmxm satisfies (1.3). By Lemma 2.2, the power series m=2cmxm is convergent for each x(-ρ1,ρ1).

Substituting m=2cmxm for y(x) in (1.3) and collecting like powers together, it follows from (2.2) and (2.3) that (with c0=c1=0) y′′(x)+ω2y(x)=m=0[(m+2)(m+1)cm+2+ω2cm]xm=m=0amxm, for all x(-ρ1,ρ1).

Therefore, every solution y:(-ρ1,ρ1) of the inhomogeneous simple harmonic oscillator equation (1.3) can be expressed by y(x)=yh(x)+m=2cmxm,where yh(x) is a simple harmonic oscillator function.

3. Partial Solution to Hyers-Ulam Stability Problem

In this section, we will investigate a property of the simple harmonic oscillator equation (2.1) concerning the Hyers-Ulam stability problem. That is, we will try to answer the question whether there exists a simple harmonic oscillator function near any approximate simple harmonic oscillator function.

Theorem 3.1.

Let y:(-ρ,ρ) be a given analytic function which can be represented by a power series m=0bmxm whose radius of convergence is at least ρ>0. Suppose there exists a constant ɛ>0 such that |y′′(x)+ω2y(x)|ɛ, for all x(-ρ,ρ) and for some positive number ω. Let ρ1=min{1,ρ}. Define am=(m+2)(m+1)bm+2+ω2bm for all m{0,1,2,} and suppose further that m=0|amxm|K|m=0amxm|, for all x(-ρ,ρ) and for some constant K. Then, there exists a simple harmonic oscillator function yh:(-ρ1,ρ1) such that |y(x)-yh(x)|Cε, for all x[-ρ0,ρ0], where ρ0<ρ1 is any positive number and C is some constant which depends on ρ0.

Proof.

We assumed that y(x) can be represented by a power series and y′′(x)+ω2y(x)=m=0amxm also satisfies m=0|amxm|K|m=0amxm|Kɛ, for all x(-ρ,ρ) from (3.1).

According to Theorem 2.3, y(x) can be written as yh(x)+m=2cmxm for all x(-ρ1,ρ1), where yh is some simple harmonic oscillator function and cm are given by (2.2). Then by Lemmas 2.1 and 2.2 and their proofs (replace M1 and M2 with Kɛ in Lemma 2.1), |y(x)-yh(x)|=|m=2cmxm|Cɛ for all x[-ρ0,ρ0], where ρ0<ρ1 is any positive number and C is some constant which depends on ρ0.

Actually from the proof of Lemma 2.1, with both M1 and M2 replaced by Kɛ, we find C1=max{Kɛωω,Kɛ} and C2=3/2C1. Further from the proof of Lemma 2.2, we have Cɛ=max{C11-ρ0,C2}=max{Kɛ1-ρ0ωω,Kɛ1-ρ0,32Kɛωω,32Kɛ},   we find C=max{K1-ρ0ωω,K1-ρ0,32Kωω,32K},which completes the proof of our theorem.

4. Example

In this section, we show that there certainly exist functions y(x) which satisfy all the conditions given in Theorem 3.1. We introduce an example related to the simple harmonic oscillator equation (1.3) for ω=1/4.

Let yh(x) be a simple harmonic oscillator function for ω=1/4 and let y:(-1,1) be an analytic function given byy(x)=yh(x)+ɛm=0x2m4m+1, where ɛ is a positive constant. (We can easily show that the radius of convergence of the power series m=0x2m/4m+1 is 2). Then, we havey′′(x)+116y(x)=m=0amxm,

wheream={4m2+12m+92m+6(for  m{0,2,4,})0(for  m{1,3,5,}),m=0|amxm|=|m=0amxm|

for any x(-1,1). So we can here choose K=1.

Furthermore, we get|y(x)+116y(x)|m=08ɛ(m+1)(2m+1)4m+3|x|2m+m=0ɛ4m+3|x|2mm=015ɛ32|x|2m2m+m=0ɛ43|x|2m4mm=015ɛ3212m+m=0ɛ4314m<ɛ,and it follows from (4.1) that|y(x)-yh(x)|=|ɛm=0x2m4m+1|ɛ4m=014m=Cɛ,

for all x(-1,1), where we set C=1/3.

CzerwikS.Functional Equations and Inequalities in Several Variables2002River Edge, NJ, USAWorld Scientific10.1142/97898127781161904790HyersD. H.On the stability of the linear functional equationProceedings of the National Academy of Sciences of the United States of America194127222224000407610.1073/pnas.27.4.222ZBL0061.26403HyersD. H.IsacG.RassiasT. M.Stability of Functional Equations in Several Variables1998Boston, Mass, USABirkhäuser Boston Inc.1639801HyersD. H.RassiasT. M.Approximate homomorphismsAequationes Mathematicae1992442-312515310.1007/BF018309751181264ZBL0806.47056JungS.-M.Hyers-Ulam-Rassias Stability of Functional Equations in Mathematical Analysis2001Palm Harbor, Fla, USAHadronic Press Inc.1841182RassiasT. M.On the stability of the linear mapping in Banach spacesProceedings of the American Mathematical Society197872229730010.2307/2042795507327ZBL0398.47040UlamS. M.Problems in Modern Mathematics1964John Wiley & Sons0280310ObłozaM.Hyers stability of the linear differential equationRocznik Naukowo-Dydaktyczny1993132592701321558ZBL0964.34514ObłozaM.Connections between Hyers and Lyapunov stability of the ordinary differential equationsRocznik Naukowo-Dydaktyczny1997141411461489723ZBL1159.34332AlsinaC.GerR.On some inequalities and stability results related to the exponential functionJournal of Inequalities and Applications19982437338010.1155/S102558349800023X1671909ZBL0918.39009TakahasiS.-E.MiuraT.MiyajimaS.On the Hyers-Ulam stability of the Banach space-valued differential equation y=λyBulletin of the Korean Mathematical Society200239230931510.4134/BKMS.2002.39.2.3091904665ZBL1011.34046MiuraT.On the Hyers-Ulam stability of a differentiable mapScientiae Mathematicae Japonicae200255117241885772MiuraT.JungS.-M.TakahasiS.-E.Hyers-Ulam-Rassias stability of the Banach space valued linear differential equations y=λyJournal of the Korean Mathematical Society2004416995100510.4134/JKMS.2004.41.6.9952095549ZBL1069.34079MiuraT.MiyajimaS.TakahasiS.-E.Hyers-Ulam stability of linear differential operator with constant coefficientsMathematische Nachrichten2003258909610.1002/mana.2003100882000046ZBL1039.34054MiuraT.MiyajimaS.TakahasiS.-E.A characterization of Hyers-Ulam stability of first order linear differential operatorsJournal of Mathematical Analysis and Applications2003286113614610.1016/S0022-247X(03)00458-X2009624ZBL1045.47037JungS.-M.Hyers-Ulam stability of linear differential equations of first orderApplied Mathematics Letters200417101135114010.1016/j.aml.2003.11.0042091847ZBL1061.34039JungS.-M.Hyers-Ulam stability of linear differential equations of first order. IIApplied Mathematics Letters200619985485810.1016/j.aml.2005.11.0042240474ZBL1125.34328JungS.-M.Hyers-Ulam stability of linear differential equations of first order. IIIJournal of Mathematical Analysis and Applications2005311113914610.1016/j.jmaa.2005.02.0252165468ZBL1087.34534JungS.-M.Hyers-Ulam stability of a system of first order linear differential equations with constant coefficientsJournal of Mathematical Analysis and Applications2006320254956110.1016/j.jmaa.2005.07.0322225975ZBL1106.34032JungS.-M.Legendre's differential equation and its Hyers-Ulam stabilityAbstract and Applied Analysis200720071410.1155/2007/56419564192365813ZBL1153.34306KimB.JungS.-M.Bessel's differential equation and its Hyers-Ulam stabilityJournal of Inequalities and Applications2007200782164010.1155/2007/216402366337ZBL1132.39023