Some Properties of Approximate Solutions of Linear Differential Equations

: In this paper, we will consider the Hyers-Ulam stability for the second order inhomogeneous linear differential equation, u (cid:48)(cid:48) ( x ) + α u (cid:48) ( x ) + β u ( x ) = r ( x ) , with constant coefﬁcients. More precisely, we study the properties of the approximate solutions of the above differential equation in the class of twice continuously differentiable functions with suitable conditions and compare them with the solutions of the homogeneous differential equation u (cid:48)(cid:48) ( x ) + α u (cid:48) ( x ) + β u ( x ) = 0. Several mathematicians have studied the approximate solutions of such differential equation and they obtained good results. In this paper, we use the classical integral method, via the Wronskian, to establish the stability of the second order inhomogeneous linear differential equation with constant coefﬁcients and we will compare our result with previous ones. Specially, for any desired point c ∈ R we can have a good approximate solution near c with very small error estimation.


Introduction
Consider the following form of the nth order differential equation, G u (n) , u (n−1) , . . . , u , u, where the function u : I → R is an n times continuously differentiable function. Given ε > 0 and n times continuously differentiable function u : I → R, we assume that G u (n) , u (n−1) , . . . , u , u, x ≤ ε (for all x ∈ I) (2) and we want to find a solution u 0 : I → R of the differential Equation (1) and a constant K > 0 such that |u(x) − u 0 (x)| ≤ Kε (for any x ∈ I).
If the above statement is true for any sufficiently small ε > 0 and any n times continuously differentiable function u : I → R and if the constant K is independent not only from ε but also from u, then the equation in (1) is said to possess the Hyers-Ulam stability property. The reader is referred to [1,2] for more details.
For a, b ∈ R, I denotes the interval [a, b] ⊂ R and we consider the second order inhomogeneous linear equation with the following form, where α and β are real-valued constants and r : I → R is a continuous function. We note that if we set r(x) ≡ 0 in (4), then the resulting differential equation, is called the homogeneous differential equation corresponding to (4). With various additional conditions, the Hyers-Ulam stability of Equation (4) was proved in [12][13][14][15][16][17]. In this paper, with weaker conditions, we will study the Hyers-Ulam stability of Equation (4) in three different cases with respect to the constant coefficients. In particular, we will try the classical integral method for error estimation. Moreover, using our results, we investigate some properties of approximate solutions of (4) in the class of the twice continuously differentiable functions, and then we compare the approximate solutions with the solutions of the homogeneous differential equation u (x) + αu (x) + βu(x) = 0 by estimating the norm of difference between them. Finally, in Section 5 we will conclude this paper.

Preliminaries
First, in the following proposition, we introduce the general solution of the second order inhomogeneous linear differential Equation (4). Proposition 1 ([18], §2.16). If a general solution u h : I → R of the homogeneous differential Equation (5) has the form u h (x) = d 1 u 1 (x) + d 2 u 2 (x) (d 1 and d 2 are arbitrary real constants), then a general solution u : I → R of the inhomogeneous linear differential Equation (4) has the form (a variation of constants or Duhamel's principle) where c is an arbitrarily chosen point of I and is the Wronskian of u 1 and u 2 .
The generalized Hyers-Ulam stability of the second order inhomogeneous linear differential equations with variable coefficients was investigated in ( [19], Corollary 3.3) in the class of twice continuously differentiable functions. In the following theorem, we introduce a version of ( [19], Corollary 3.3) for the constant coefficients case. The main results in the paper, Theorems 2-4, can be considered as corollaries to ([19], Corollary 3.3). In Theorems 5-7, these results are refined for r(x) = 0.
Let I be a (half)-open or closed subinterval of R and let c ∈ I. Theorem 1. Let r : R → R be a given continuous function, let α and β be real-valued constants. Suppose a general solution u h : I → R of the homogeneous differential Equation (5) is of the form u h (x) = d 1 u 1 (x) + d 2 u 2 (x), where d 1 and d 2 are arbitrary real-valued constants. Let u : I → R be a twice continuously differentiable function and satisfies the inequality for all x ∈ I, where the following integral exists with σ : I → [0, ∞), then there is a solution u 0 : I → R of differential Equation (4) satisfying for every x ∈ I, where c is arbitrarily given in I.

Differential Equations with Constant Coefficients
Now, we investigate the second order inhomogeneous linear differential Equation (4) with constant coefficients, where α and β are real-valued constants and r : R → R is a continuous function.
We denote by m 1 and m 2 the roots of the characteristic equation According to Theorem 1, the general solution of (4) is given as for α 2 − 4β > 0; or for for α 2 − 4β < 0, where c is a real number. In the following theorems, we investigate the Hyers-Ulam stability of the second order inhomogeneous differential Equation (4), with constant coefficients, defined on the whole space R.

Theorem 2.
Let r : R → R be a given continuous function, let α and β be given real-valued constants satisfying α 2 − 4β > 0, let c ∈ R, and let ε > 0 be given. If a twice continuously differentiable function u : R → R satisfies the inequality for every x ∈ R, then there exist real-valued constants d 1 and d 2 such that for all x ∈ R, where m 1 and m 2 are given in (6).
It further follows from the last inequality that for all x ∈ R, from which we get the inequality (10) by applying the l'Hospital's rule.
On the other hand, since u c is a solution of (4) with α 2 − 4β > 0, it follows from (7) that u c has the form given in the statement of our theorem.
We now consider the case when α 2 − 4β = 0, i.e., β = α 2 4 > 0 and prove the Hyers-Ulam stability of the inhomogeneous differential Equation (4). Theorem 3. Let r : R → R be a continuous function, let α = 0 be a real-valued constant, and let c ∈ R and ε > 0 be arbitrary constants. If a twice continuously differentiable function u : R → R satisfies the inequality for every x ∈ R, then there exist real-valued constants d 1 and d 2 such that for every x ∈ R, where λ is given in (6).
According to Theorem 1 with σ(x) = ε and I = R, there is a solution u c : R → R of (4) satisfying It further follows from the last inequality that for all x ∈ R. Finally, we know that u c is a solution of (4) with α 2 − 4β = 0. Therefore, in view of (8), u c has the form given in the statement of our theorem. Now, we are in the position to prove the Hyers-Ulam stability of the inhomogeneous differential Equation (4) for the case when α 2 − 4β < 0. Theorem 4. Let r : R → R be a continuous function, let α and β be real-valued constants satisfying α 2 − 4β < 0, and let c ∈ R and ε > 0 be arbitrary constants. If a twice continuously differentiable function u : R → R satisfies the inequality for all x ∈ R, then there exist real-valued constants d 1  for all x ∈ R, where λ and µ are given in (6).
Proof. Since α 2 − 4β < 0, we have u 1 (x) = e λx cos µx, u 2 (x) = e λx sin µx, and W(u 1 , u 2 )(x) = µe 2λx . According to Theorem 1 with σ(x) = ε and I = R, there is a solution u c : R → R of (4) satisfying for all x ∈ R. If we change the variable with τ = x − t, then we get On the other hand, since u c is a solution of (4) with α 2 − 4β < 0, it follows from (9) that u c has the form given in the statement of our theorem.
For any x ∈ R, let us denote by n(x) = [ µx π ] the integer satisfying If x ≥ 0, then we have If we replace λ by −λ in (15), then we get For an arbitrary x < 0, it follows from (16) that Remark 1. If x ≥ c in Theorem 4, then we use (15) and if x < c, then we use (17) to get

Approximation Properties
Let α and β be real-valued constants, and let the real vector space B(α; β) consist of those C 2 -functions u : R → R for which there exists a finite constant M (but depending on u) such that for every x ∈ R. It is clear that B(α; β) is a vector space over R. Let for all v 1 , v 2 ∈ B(α; β) and γ ∈ R, then B(α; β) is a vector space over R. Therefore, we may remark that B(α; β) is a large set enough to be a vector space.
In the next theorems, we study an approximation property of functions from B(α; β) by solutions of differential Equation (5).

Theorem 5.
Let α and β be real-valued constants satisfying α 2 − 4β > 0 and let c ∈ R be given. If u ∈ B(α; β), then there exist real-valued constants p and q such that where m 1 and m 2 are given in (6).
Proof. Since α 2 − 4β > 0, we have m 1 < m 2 , u 1 (x) = e m 1 x , and u 2 (x) = e m 2 x . Since u ∈ B(α; β), there is a constant ε > 0 such that (18) is true for every x ∈ R. On account of Theorem 2 with r(x) ≡ 0 and the inequality (11), there exist real-valued constants p and q such that for all x ∈ R. Thus, we have which ends our proof.
When α 2 − 4β = 0, for the proof of the following theorem, we apply Theorem 3.

Theorem 6.
Let α = 0, β = α 2 4 , and c ∈ R be real-valued constants. If y ∈ B(α; β), then there exist real-valued constants p and q such that where λ is given in (6). (18) is true for every x ∈ R. On account of Theorem 3 with r(x) ≡ 0, there exist real-valued constants p and q such that Moreover, we apply the l'Hospital's rule to get which completes the proof.
Finally, we investigate an approximation property of the approximate solutions for each function in B(α; β) by some solution of the differential Equation (5) when α 2 − 4β < 0.

Theorem 7.
Let α and β be real-valued constants satisfying α 2 − 4β < 0 and let c ∈ R be a given constant. If u ∈ B(α; β), then there exist real-valued constants p and q such that u(x) − pe λx cos µx − qe λx sin µx = o |x − c| as x → c, where λ and µ are given in (6).
Proof. Since α 2 − 4β < 0, we have u 1 (x) = e λx cos µx and u 2 (x) = e λx sin µx. By the definition (14) of n(x), we know that Since u ∈ B(α; β), a finite constant M > 0 exists such that the inequality in (18) is satisfied for every x ∈ R. For both cases c − π µ < x < c and c ≤ x < c + π µ , we apply Theorem 4 with r(x) ≡ 0 and Remark 1 and we consider the last argument to see that there exist real-valued constants p and q such that for all c − π µ < x < c + π µ . Furthermore, we apply the l'Hospital's rule to calculate which completes our proof.

Conclusions
In the main theorems, we investigated the Hyers-Ulam stability of second order inhomogeneous linear differential Equation (4), with constant coefficients, defined on the whole space R, while previous papers have investigated the local Hyers-Ulam stability of the second order linear differential equations defined on the bounded intervals like I = (a, b) or I = [a, b] with −∞ < a < b < +∞ under some additional conditions.
Ghaemi et al. [13] recently proved a theorem concerning the Hyers-Ulam stability of the exact second order linear differential equations which yields the following inequality instead of (10) under the additional assumptions that m = α+ √ α 2 −4β 2 > 0 and the domain of u and u c is a bounded interval I = (a, b) (see [13], Corollary 2.3). Our results in (10) contain more generalized cases than those in [13]. Moreover, the inequality (19) does not seem to be better than the inequality (10), while we can find an approximate solution with very small error near c for any desired point c ∈ R. Cîmpean and Popa [12] proved Hyers-Ulam stability of the linear differential equations with constant coefficients. From corollaries in [12], one can obtain the following inequalities (for either α 2 − 4β > 0 and β = 0), 4ε α 2 (for α 2 − 4β ≤ 0 and α = 0).
However, our results in Theorems 2-4 are a better approximation than theirs. Moreover, by our results in Theorems 5-7, for any desired c ∈ R, we can obtain an approximation with very small error.
Finally, compared to the results in [14][15][16] we used very classical integral method and have better error estimation than theirs. While we can find an approximate solution generally with very small error near c for any desired point c ∈ R, in their estimates this is only possible in special cases.

Remark 2.
It would be much nicer if we could provide an estimate of the form |u(x) − u c (x)| ≤ C ε µ x−c 0 e λt sin µtdt , x ∈ I for some constant C, instead of (13) in Theorem 4. We will investigate this problem in next paper.