Hyers–Ulam–Rassias Stability of Hermite’s Differential Equation

Abstract

In this paper, we studied the Hyers–Ulam–Rassias stability of Hermite’s differential equation, using Pachpatte’s inequality. We compared our results with those obtained by Blaga et al. Our estimation for $\left|z\left(x\right)-y\left(x\right)\right|$, where z is an approximate solution and y is an exact solution of Hermite’s equation, was better than that obtained by the authors previously mentioned, in some parts of the domain, especially in a neighborhood of the origin.

1. Introduction

The importance of studying Hermite’s equation results from the multiple applications that the equation has in theoretical physics, electrical engineering, economics, etc. This equation has the form:

$$y^{″}\left(x\right)-2x{y}^{\prime }\left(x\right)+\lambda y\left(x\right)=0,x\in \mathbb{R},y\in C^2(\mathbb{R},\mathbb{R}),\lambda \in \mathbb{R}.$$

If $\lambda =2n,n\in \mathbb{N}$, Hermite’s equation becomes:

$$y^{″}\left(x\right)-2x{y}^{\prime }\left(x\right)+2ny\left(x\right)=0,x\in \mathbb{R},y\in C^2(\mathbb{R},\mathbb{R}),n\in \mathbb{N},$$

(1)

and admits particular polynomial solutions $H_n$, called Hermite’s polynomials (see [1]), given by:

$$H_n\left(x\right)={\left(-1\right)}^n{e}^{x^2}{\left(e^{-x^2}\right)}^{\left(n\right)},x\in \mathbb{R}.$$

In what follows, we study the Hyers–Ulam–Rassias and generalized Hyers–Ulam–Rassias stability of Equation (1), using Pachpatte’s inequality. The problem was also studied in [2] by Blaga et al., as an application of a result regarding the Ulam stability of a linear system of differential equations with nonconstant coefficients. We compared our results with those obtained in [2]. Pachpatte’s next inequality is used in the paper.

Theorem 1

(Pachpatte [3], p. 39, Theorem 1.7.4.). Let $u,f,g$ be nonnegative continuous functions defined on ${\mathbb{R}}_{+},$ and h be a positive and nondecreasing continuous function defined on ${\mathbb{R}}_{+},$ for which the inequality:

$$u\left(x\right)\le h\left(x\right)+{\int }_0^{x}f\left(t\right)u\left(t\right)dt+{\int }_0^{x}f\left(t\right)\left[{\int }_0^{t}g\left(s\right)u\left(s\right)ds\right]dt,$$

holds for $x\in {\mathbb{R}}_{+}.$ Then:

$$u\left(x\right)\le h\left(x\right)\left(1+{\int }_0^{x}f\left(t\right)e^{{\int }_0^t\left(f\left(s\right)+g\left(s\right)\right)ds}dt\right),$$

for $x\in {\mathbb{R}}_{+}.$

We recall that Ulam stability was proposed by Ulam in [4]. Hyers [5] established the first result in this direction. Obloza [6] established the first result regarding the Ulam stability of differential equations and also investigated the relation between Lyapunov and Ulam stability. After that, the domain has developed rapidly. The Hyers–Ulam stability of linear differential equations of first order was studied for example in [7,8,9,10,11,12], of higher order in [13], linear difference equations in [14], and others in [15,16]. More on this topic can be found in the books [17,18]. In [19], linear differential equations of second order of the form:

$$y^{″}+\alpha y^{\prime }+\beta y=0\mathrm{and}{y}^{″}+\alpha y^{\prime }+\beta y=f\left(x\right),$$

$y\in C^2[a,b],f\in C[a,b],-\infty <a<b<\infty ,$ were studied.

In [20], the equation:

$$y^{″}+p\left(x\right)y^{\prime }+q\left(x\right)y+r\left(x\right)=0,$$

$y\in C^2(a,b),f\in C[a,b],p,q,r\in C(a,b),-\infty <a<b<\infty ,$ was investigated.

Various other types of differential equations of second order have been further studied, for example in [21], the equation:

$$y^{″}\left(x\right)+(1+\psi \left(x\right))y\left(x\right)=0,\psi \in C^1[0,\infty ),y\in C^2[0,\infty ),$$

in [22,23], the equation:

$$y^{″}\left(x\right)+\beta \left(x\right)y\left(x\right)=0,\beta \in C[a,b],y\in C^2[a,b],$$

in [22] with boundary conditions, in [23] with initial conditions, and in [24], the following equation:

$$y^{″}\left(x\right)+{\mu }^{2}y\left(x\right)=0,y\in C^2[a,b],-\infty <a<b<\infty ,$$

here using the Mahgoub transform.

In [25], the generalized Hyers–Ulam stability of the equation:

$$y^{″}+p\left(x\right)y^{\prime }+q\left(x\right)y=f\left(x\right),y\in C^2\left(I\right),p,q,f\in C\left[I\right],I\subset \mathbb{R},$$

was studied, if there exists a solution $y_1\in C^2\left(I\right)$ of the equation:

$$y^{″}+p\left(x\right)y^{\prime }+q\left(x\right)y=0.$$

As an application, some remarks regarding generalized the Hyers–Ulam–Rassias stability of Hermite’s equation are given. In [26], the equation:

$$y^{″}\left(x\right)+f(x,y\left(x\right))=0,\mathrm{if}\left|f(x,y(x\left)\right)\right|\le A{\left|y\right|}^{\alpha },\mathrm{and}f(x,0)=0,$$

was investigated, using the Gronwall lemma. In [27], the equations:

$$y^{″}\left(x\right)+f(x,y\left(x\right))=0,$$

and:

$$y^{″}\left(x\right)+f(x,y\left(x\right),y^{\prime }\left(x\right))=0,$$

$y\in C^2\left(I\right),f\in C\left(I\right),I\subset \mathbb{R},$ were studied, using the Bihari lemma. The results obtained in [27] generalize those established in [26].

The outline of the paper is the following: In Section 2, we present the stability notions and prove several auxiliary results useful in the following sections (Remarks 1–4). The first main result (Theorem 2) is given in the next section and concerns the Hyers–Ulam–Rassias stability of Hermite’s Equation (2). The second main result (Theorem 3), regarding the generalized Hyers–Ulam–Rassias stability of Hermite’s Equation (2), is proven in Section 4. A graphical representation illustrating the difference between our result in Theorem 2 and that established by Blaga et al. in [2] is shown in Figure 1.

2. Preliminary Notions and Results

In what follows, we consider $n\in \mathbb{N},I=\left[0,n+1\right],\alpha \in \mathbb{R},\beta \in \mathbb{R}$, the equation:

$$y^{″}\left(x\right)-2x{y}^{\prime }\left(x\right)+2ny\left(x\right)=0,x\in I,y\in C^2(I,\mathbb{R}),$$

(2)

and the conditions:

$$\left\{\begin{array}{c}y\left(0\right)=\alpha \\ y^{\prime }\left(0\right)=\beta \end{array}\right..$$

(3)

Let $\epsilon >0.$ Let $\phi \in C^2(I,{\mathbb{R}}_{+})$. We also consider the inequalities:

$$\left|z^{″}\left(x\right)-2x{z}^{\prime }\left(x\right)+2nz\left(x\right)\right|\le \epsilon ,x\in I,z\in C^2(I,\mathbb{R}),$$

(4)

and:

$$\left|z^{″}\left(x\right)-2x{z}^{\prime }\left(x\right)+2nz\left(x\right)\right|\le \epsilon \phi \left(x\right),x\in I,z\in C^2(I,\mathbb{R}).$$

(5)

Definition 1.

Equation (2) is called Hyers–Ulam–Rassias stable if there is a real number $c>0$ and a function $\psi :I\to \left[0,\infty \right)$ so that for any solution z of the inequality (4), satisfying (3), there is a solution y of Equation (2), satisfying (3), such that:

$$\left|z\left(x\right)-y\left(x\right)\right|\le c·\epsilon ·\psi \left(x\right),\forall x\in I.$$

Definition 2.

Equation (2) is called generalized Hyers–Ulam–Rassias stable if there is a real number $c>0$ and a function $\psi :[a,b]\to \left[0,\infty \right)$ so that for any solution z of the inequality (5), satisfying (3), there is a solution y of Equation (2), satisfying (3), such that:

$$\left|z\left(x\right)-y\left(x\right)\right|\le c·\epsilon ·\psi \left(x\right),\forall x\in I.$$

In Definitions 1 and 2, z is called an approximate solution and y is called an exact solution of Hermite’s Equation (2).

Remark 1.

A function $z\in C^2(I,\mathbb{R})$ is a solution of (4) if and only if there exists a function $g\in C^2(I,\mathbb{R})$ such that:

 1.

$\left|g\left(x\right)\right|\le \epsilon ,\forall x\in I;$

 2.

$z^{″}\left(x\right)=2x{z}^{\prime }\left(x\right)-2nz\left(x\right)+g\left(x\right),\forall x\in I.$

Remark 2.

If a function $z\in C^2(I,\mathbb{R})$ is a solution of (4), then z is a solution of the inequality:

$$\left|z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x-2{\int }_0^{x}tz\left(t\right)dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt\right|\le \epsilon \frac{x^2}{2},\forall x\in I.$$

Indeed, integrating $z^{″}\left(x\right)=2x{z}^{\prime }\left(x\right)-2nz\left(x\right)+g\left(x\right),$ from zero to x, we have:

$$z^{\prime }\left(x\right)-z^{\prime }\left(0\right)={\int }_0^{x}2s{z}^{\prime }\left(s\right)ds-2n{\int }_0^{x}z\left(s\right)ds+{\int }_0^{x}g\left(s\right)ds;$$

hence:

$$z^{\prime }\left(x\right)-z^{\prime }\left(0\right)=2xz\left(x\right)-{\int }_0^{x}2z\left(s\right)ds-2n{\int }_0^{x}z\left(s\right)ds+{\int }_0^{x}g\left(s\right)ds,$$

and:

$$z^{\prime }\left(x\right)-z^{\prime }\left(0\right)=2xz\left(x\right)-2\left(n+1\right){\int }_0^{x}z\left(s\right)ds+{\int }_0^{x}g\left(s\right)ds.$$

Integrating again from zero to x, we have:

$$z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x=2{\int }_0^{x}tz\left(t\right)dt-2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt+{\int }_0^x\left[{\int }_0^{t}g\left(s\right)ds\right]dt.$$

Hence, we obtain:

$$\begin{array}{cc}\hfill & \left|z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x-2{\int }_0^{x}tz\left(t\right)dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt\right|\hfill \\ \hfill & \le {\int }_0^x\left[{\int }_0^t\left|g\left(s\right)\right|ds\right]dt\le \epsilon \frac{x^2}{2}.\hfill \end{array}$$

Remark 3.

A function $z\in C^2(I,\mathbb{R})$ is a solution of (5) if and only if there exists a function $g\in C^2(I,\mathbb{R})$ such that:

 1.

$\left|g\left(x\right)\right|\le \epsilon \phi \left(x\right),\forall x\in I;$

 2.

$z^{″}\left(x\right)=2x{z}^{\prime }\left(x\right)-2nz\left(x\right)+g\left(x\right),\forall x\in I,$

where $\phi \in C^2(I,{\mathbb{R}}_{+})$.

Analogous to Remark 2, we can prove the following result.

Remark 4.

If a function $z\in C^2(I,\mathbb{R})$ is a solution of (5), then z is a solution of the inequality:

$$\left|z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x-2{\int }_0^{x}tz\left(t\right)dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt\right|\le \epsilon {\int }_0^x\left[{\int }_0^t\phi \left(s\right)ds\right]dt,$$

$\forall x\in I.$

3. Hyers–Ulam–Rassias Stability of Hermite’s Equation

We present below a result regarding the Hyers–Ulam–Rassias stability of Hermite’s Equation (2), with the conditions (3). We also make a comparison between our results and those set out in [2].

Theorem 2.

Let $n\in \mathbb{N},I=\left[0,n+1\right],\alpha \in \mathbb{R},\beta \in \mathbb{R}.$ Then:

 1.

Equation (2) has a unique solution with:

$$\left\{\begin{array}{c}y\left(0\right)=\alpha \\ y^{\prime }\left(0\right)=\beta \end{array}\right.;$$

 2.

Hermite’s Equation (2), with the conditions (3), is Hyers–Ulam–Rassias stable.

Proof. 

1. It is a well-known result from differential equations of order two, since the coefficients of Equation (2), $p\left(x\right)=-2x,q\left(x\right)=2n,$ are continuous functions on I (see [28]);

2. Let z be a solution of (4) and y be the solution of (2) such that:

$$\left\{\begin{array}{c}z\left(0\right)=y\left(0\right)=\alpha \\ z^{\prime }\left(0\right)=y^{\prime }\left(0\right)=\beta \end{array}\right..$$

From Remark 2, we have:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \left|z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x-2{\int }_0^{x}tz\left(t\right)dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt\right|\hfill \\ \hfill & +2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt\hfill \\ \hfill & \le \epsilon \frac{x^2}{2}+2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt.\hfill \end{array}$$

Since $t\le x\le n+1$, we have:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \epsilon \frac{x^2}{2}+2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt\hfill \\ \hfill & \le \epsilon \frac{x^2}{2}+{\int }_0^{x}2\left(n+1\right)\left|z\left(t\right)-y\left(t\right)\right|dt+{\int }_0^{x}2\left(n+1\right)\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt.\hfill \end{array}$$

Using now Pachpatte’s inequality from Theorem 1, we obtain:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \epsilon \frac{x^2}{2}\left(1+{\int }_0^{x}2\left(n+1\right)e^{{\int }_0^t\left(2n+3\right)ds}dt\right)=\epsilon \frac{x^2}{2}\left(1+{\int }_0^{x}2\left(n+1\right)e^{\left(2n+3\right)t}dt\right)\hfill \\ \hfill & =\epsilon \frac{x^2}{2}\left(1+\frac{2n+2}{2n+3}\left(e^{\left(2n+3\right)x}-1\right)\right)\hfill \\ \hfill & =\epsilon \frac{x^2}{2}\left(\frac{1}{2n+3}+\frac{2n+2}{2n+3}{e}^{\left(2n+3\right)x}\right).\hfill \end{array}$$

□

In the following, using the graphical representation from Figure 1, we make a comparison between the estimation $\left|z\left(x\right)-y\left(x\right)\right|$ obtained by Blaga et al. in [2] and that obtained by us in Theorem 2. For $n=0$ we have $I=[0,1]$, and we compare the functions $f_1\left(x\right)=\sqrt{\pi }{e}^{x^2+x}$, which appeared in [2], and $f_2\left(x\right)=x^2\left(\frac{1}{3}+\frac{2}{3}{e}^{3x}\right),$ which appears in Theorem 2. In Figure 1 we graphically represent both functions. It can be seen that in some parts of the domain, our estimation was better than the one obtained in [2], especially in a neighborhood of the origin.

For other values of n, we have similar results.

4. Generalized Hyers–Ulam–Rassias Stability of Hermite’s Equation

We establish now a result regarding the generalized Hyers–Ulam–Rassias stability of Hermite’s Equation (2), with the conditions (3).

Theorem 3.

Let $n\in \mathbb{N},I=\left[0,n+1\right],\alpha \in \mathbb{R},\beta \in \mathbb{R}.$ If:

 (i)

$\exists {\lambda }_{\phi }>0$ such that ${\int }_0^x\left[{\int }_0^t\phi \left(s\right)ds\right]dt\le {\lambda }_{\phi }\phi \left(x\right),\forall x\in I;$

 (ii)

$\phi \in C^2(\left[0,\infty \right),{\mathbb{R}}_{+})$ is increasing;

Then Hermite’s Equation (2), with the conditions (3), is generalized Hyers–Ulam–Rassias stable.

Proof. 

Let z be a solution of (5) and y be the solution of (2) such that:

$$\left\{\begin{array}{c}z\left(0\right)=y\left(0\right)=\alpha \\ z^{\prime }\left(0\right)=y^{\prime }\left(0\right)=\beta \end{array}\right..$$

From Remark 4, we have:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \left|z\left(x\right)-z\left(0\right)-z^{\prime }\left(0\right)x-2{\int }_0^{x}tz\left(t\right)dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^{t}z\left(s\right)ds\right]dt\right|\hfill \\ \hfill & +2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt\hfill \\ \hfill & \le \epsilon {\int }_0^x\left[{\int }_0^t\phi \left(s\right)ds\right]dt+2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt.\hfill \end{array}$$

Since $t\le x\le n+1$, we have:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \epsilon {\lambda }_{\phi }\phi \left(x\right)+2{\int }_0^{x}t\left|z\left(t\right)-y\left(t\right)\right|dt+2\left(n+1\right){\int }_0^x\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt\hfill \\ \hfill & \le \epsilon {\lambda }_{\phi }\phi \left(x\right)+{\int }_0^{x}2\left(n+1\right)\left|z\left(t\right)-y\left(t\right)\right|dt+{\int }_0^{x}2\left(n+1\right)\left[{\int }_0^t\left|z\left(t\right)-y\left(t\right)\right|ds\right]dt.\hfill \end{array}$$

Using now Pachpatte’s inequality from Theorem 1, we obtain:

$$\begin{array}{cc}\hfill \left|z\left(x\right)-y\left(x\right)\right|& \le \epsilon {\lambda }_{\phi }\phi \left(x\right)\left(1+{\int }_0^{x}2\left(n+1\right)e^{{\int }_0^t\left(2n+3\right)ds}dt\right)=\epsilon \frac{x^2}{2}\left(1+{\int }_0^{x}2\left(n+1\right)e^{\left(2n+3\right)t}dt\right)\hfill \\ \hfill & =\epsilon {\lambda }_{\phi }\phi \left(x\right)\left(1+\frac{2n+2}{2n+3}\left(e^{\left(2n+3\right)x}-1\right)\right)\hfill \\ \hfill & =\epsilon {\lambda }_{\phi }\phi \left(x\right)\left(\frac{1}{2n+3}+\frac{2n+2}{2n+3}{e}^{\left(2n+3\right)x}\right).\hfill \end{array}$$

□

5. Conclusions

In this paper, we studied Hyers–Ulam–Rassias and generalized Hyers–Ulam–Rassias stability of Hermite’s equation, using a new approach, namely Pachpatte’s inequality. The problem was also considered in [2] by Blaga et al., as a particular case of a general result regarding the Ulam stability of a linear system of differential equations with nonconstant coefficients. We remark, using the graphical representation from Figure 1, that our estimation for $\left|z\left(x\right)-y\left(x\right)\right|$, between an approximate solution and the exact solution of Equation (2), in some parts of the domain, was better than the one obtained in [2], especially in the neighborhood of the origin.