# Wavelet-Galerkin Quasilinearization Method for Nonlinear Boundary Value Problems

Authors: Umer Saeed; Mujeeb ur Rehman
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

A numerical method is proposed by wavelet-Galerkin and quasilinearization approach for nonlinear boundary value problems. Quasilinearization technique is applied to linearize the nonlinear differential equation and then wavelet-Galerkin method is implemented to linearized differential equations. In each iteration of quasilinearization technique, solution is updated by wavelet-Galerkin method. In order to demonstrate the applicability of proposed method, we consider the various nonlinear boundary value problems.


---

## Body

## 1. Introduction

The Galerkin method [1] is a very well-known method for finding the numerical solutions of differential equations. According to wavelet-Galerkin method, connection coefficients are the inner products of Daubechies scaling functions and their derivatives, because we are taking Daubechies scaling functions as a Galerkin basis. The exact and explicit representations of the differential operators in orthonormal bases of compactly supported wavelets are described by Beylkin [2] and he also discussed the sparse representations of shift operators in orthonormal bases of compactly supported wavelets. Latto et al. [3] gave the connection coefficients for zeroth level of resolution; these connection coefficients are essentially based on an unbounded domain. Chen et al. [4] provided the way of calculating the connection coefficients on a bounded interval and these finite integrals play a vital role in the wavelet-Galerkin approximation of differential equations. Restrepo and Leaf [5] reviewed the inner products of Daubechies wavelets and their derivatives and used the connection coefficients for approximation of differential operators.Amaratunga et al. [6] implemented the wavelet-Galerkin technique for solving the one-dimensional counterpart of Helmholtz’s equation. Mishra and Sabina [7] used the wavelet-Galerkin method for solving linear, homogeneous boundary value problems with constant coefficients and compared the obtained solution with the exact solution by using a family of Daubechies wavelets and at different levels of resolution. Daubechies scaling functions as the Galerkin bases were used by Jianhua et al. [8] and the authors implemented the wavelet-Galerkin method for differential equations with a boundary layer. In order to implement the wavelet-Galerkin method, they considered the linear, inhomogeneous second order boundary value problem with constant coefficients. In [9–13] several applications of wavelet-Galerkin method are done and these applications are for linear boundary value problems. Motivated by the work of authors [6–13], we extended the wavelet-Galerkin method for the solution of nonlinear boundary value problems.The quasilinearization approach was introduced by Kalaba and Bellman [14, 15] as a generalization of the Newton-Raphson method [16] to solve individual or systems of nonlinear ordinary and partial differential equations. The quasilinearization approach is suitable for a general nonlinear ordinary or partial differential equations of any order.Jiwari [19] used a uniform Haar wavelet method with quasilinearization technique for the approximate solution of Burgers’ equation and compared the results with the solutions obtained by the other numerical methods and the exact solution. The same approach was used by Kaur et al. [20] for the solutions of nonlinear boundary value problems in which they treated the quadratic nonlinearity of unknown function. In [21], we extend the Haar wavelet quasilinearization technique for fractional nonlinear initial and boundary value problems.In this paper, we consider the compactly supported Daubechies scaling functions [22] as a Galerkin basis and propose a numerical method by combining wavelet-Galerkin method with quasilinearization technique for solving nonlinear boundary value problems. The method deals with not only quadratic nonlinearities but also various other forms of nonlinearities. To the best of our knowledge, the wavelet-Galerkin method, with Daubechies scaling functions as Galerkin basis, has not been implemented for numerical solutions of nonlinear differential equations. Illustrative problems show the advantage of the proposed method.The paper is arranged as follows: in Section2, we give a brief review of Daubechies’ wavelets, while in Section 3 we describe the two-term connection coefficients. In Section 4, we give the procedure of implementation of the wavelet-Galerkin method. In Section 5, the quasilinearization technique for dealing with ordinary differential equations is discussed, while in Section 6 we apply the wavelet-Galerkin method with quasilinearization technique to some nonlinear boundary value problems. Finally, in Section 7 we conclude our work.


## 2. Daubechies’ Wavelets

Daubechies [22, 23] constructed a family of compactly supported orthonormal wavelets. A wavelet system consists of a scaling function ϕ
(
x
) and a wavelet function ψ
(
x
). There are two important relations in wavelet theory, which we called two-scale relation. Consider
(1)
ϕ
(
x
)
=
∑
l
=
0
N
-
1
‍
p
l
ϕ
l
(
2
x
)
,
and the equation
(2)
ψ
(
x
)
=
∑
l
=
2
-
N
1
(
-
1
)
l
p
1
-
l
ϕ
l
(
2
x
)
,
where ϕ
l
(
·
)
:
=
ϕ
(
·
-
l
). Relations (1) and (2) are known as refinement relations. The coefficients p
l are called the wavelet filter coefficients. Relations (1) and (2) are also called refinement relations and N (an even integer) is the number of wavelet filter coefficients in the refinement relations. The supports of the scaling function ϕ
(
x
) and wavelet function ψ
(
x
) are [
0
,
N
-
1
] and [
1
-
N
/
2
,
N
/
2
], respectively. Daubechies [22] constructed wavelet filter coefficients p
l to satisfy the certain conditions. These conditions are linked with certain properties of scaling and wavelet functions [4]. Consider
(3)
∑
l
=
0
N
-
1
p
l
=
2
⟹
∫
-
∞
∞
ϕ
(
x
)
d
x
=
1
,
(4)
∑
l
=
0
N
-
1
p
l
p
l
-
2
m
=
2
δ
0
,
m
⟹
∫
-
∞
∞
ϕ
(
x
)
ϕ
(
x
-
m
)
d
x
=
δ
0
,
m
,
(5)
∑
l
=
2
-
N
1
(
-
1
)
l
p
1
-
l
p
l
-
2
m
=
0
⟹
∫
-
∞
∞
ϕ
(
x
)
ψ
(
x
-
m
)
d
x
=
0
,
for
any
integer
m
,
(6)
∑
l
=
0
N
-
1
(
-
1
)
l
l
m
p
l
=
0
⟹
∫
-
∞
∞
x
m
ψ
(
x
)
d
x
=
0
,
m
=
0,1
,
…
,
N
2
-
1
.Relation (3) shows that scaling functions have unit area and relations (4) and (5) indicate the orthonormality of ϕ and orthogonality of ϕ and ψ, respectively. Relation (6) shows thatmth moment of ψ is zero; that is, it has m vanishing moments, which implies that we can express the elements of the set 1
,
x
,
…
,
x
N
/
2
-
1 as a linear combination of ϕ
(
x
-
k
), integer translate of ϕ
(
x
).Daubechies wavelet has no explicit expressions for the scaling functionϕ
(
x
) and the wavelet function ψ
(
x
) at arbitrary x.The simplest Daubechies wavelet [24] is the Haar wavelet, which has explicit expression for calculating the scaling function ϕ
(
x
). It is also called D2, which means the Daubechies wavelet having two filter coefficients. It is also called db1, which means the Daubechies wavelet having one vanishing moment. db2 (D4) has four wavelet coefficients, that is, N
=
4, and two vanishing moments. Similarly, DN has N-coefficients and N
/
2-vanishing moments.


## 3. Two-Term Connection Coefficients

In the present work, we are dealing with the second order nonlinear boundary value problems. We will be concerned with the 2-term connection coefficients. Two-term connection coefficients are defined as(7)
Ω
k
d
1
,
d
2
=
∫
-
∞
∞
ϕ
d
1
(
x
)
ϕ
d
2
(
x
-
k
)
d
x
.Taked times derivative of the Daubechies scaling function (1), by assuming that it is d times differentiable, to obtain
(8)
ϕ
d
(
x
)
=
2
d
∑
l
=
0
N
-
1
p
l
ϕ
l
d
(
2
x
)
,
where
ϕ
d
(
x
)
≔
d
d
ϕ
d
x
d
.Use (8) in (7) and changing variables, to obtain [3]
(9)
Ω
k
d
1
,
d
2
=
2
d
1
+
d
2
-
1
∑
m
,
l
p
m
p
l
-
2
k
+
m
∫
-
∞
∞
ϕ
d
1
(
x
)
ϕ
d
2
(
x
-
l
)
d
x
.Let⋀
d
1
,
d
2 be a column vector with 2
N
-
3 components which are connection coefficients: ⋀
d
1
,
d
2
=
[
Ω
k
d
1
,
d
2
]
k
=
1
:
2
N
-
3. Equation (9) gives a system of linear equations with ⋀
d
1
,
d
2 as unknown vector; we can write (9) in vector form as
(10)
1
2
d
1
+
d
2
-
1
⋀
d
1
,
d
2
=
T
⋀
d
1
,
d
2
,
where T is a square matrix of order 2
N
-
3; that is, T
k
:
l
=
∑
m
=
0
N
-
1
p
m
p
l
-
2
k
+
m, where indices k and l vary from 1 to 2
N
-
3. We use the substitution [13] y
=
2
j
x
,
C
l
=
2
j
/
2
c
l, throughout our work. This substitution provides a way of calculating the connection coefficients. Here, j and l are integers and 2
j is scaling factor. It corresponds to either the expansion (
j
>
0
) or the contraction (
j
<
0
) of the scaling or wavelet function.Define connection coefficients as(11)
Ω
k
0
,
d
=
2
d
j
∫
ϕ
(
0
)
(
y
)
ϕ
d
(
y
-
k
)
d
y
.Similarly, we can obtain(12)
1
2
d
-
1
⋀
0
,
d
=
T
⋀
0
,
d
.It is homogeneous system and thus does not have a unique nonzero solution. In order to make the system inhomogeneous, one equation is added which is derived from the moment equation of the scaling function [3] as follows:
(13)
x
d
=
∑
i
m
i
d
2
j
/
2
ϕ
i
(
2
j
x
)
,
where m
i
d is thedth moment of ϕ
j
,
i
(
x
)
:
=
ϕ
(
2
j
x
-
i
) and we can compute it by considering the orthonormality of ϕ
(
x
) that is
(14)
m
i
d
=
∫
x
d
2
j
/
2
ϕ
i
(
2
j
x
)
d
x
.Considering the substitutiony
=
2
j
x in (13), we arrive at
(15)
y
d
=
2
j
d
2
j
/
2
∑
i
m
i
d
ϕ
i
(
y
)
.Differentiate (15) d times to get
(16)
d
!
=
2
j
d
2
j
/
2
∑
i
m
i
d
ϕ
i
d
(
y
)
.Taking inner product on both sides of (16) with ϕ
(
y
-
p
),
(17)
d
!
=
2
j
/
2
∑
i
m
i
d
2
j
d
∫
ϕ
i
d
(
y
)
ϕ
p
(
y
)
d
y
,
or
(18)
d
!
2
-
j
/
2
=
∑
i
m
i
d
⋀
d
,
0
.Equation (14) implies
(19)
m
i
d
=
∫
y
d
2
j
d
2
j
/
2
ϕ
i
(
y
)
d
y
2
j
,
m
i
d
=
2
-
j
d
-
j
/
2
∫
y
d
ϕ
i
(
y
)
d
y
,
m
i
d
=
2
-
j
d
-
j
/
2
M
i
d
,
where M
i
d
=
∫
‍
x
d
ϕ
i
(
x
)
d
x is thedth moment of ϕ
i and M
0
0
=
1 implies the unit area under ϕ. Equation (19) shows thatdth moment of ϕ
j
,
i
(
x
) is equal to the 2
-
j
d
-
j
/
2 timesdth moment of ϕ
(
x
-
i
).Latto et al. [3] derive an explicit formula to compute the moments of ϕ
(
x
). Consider
(20)
M
i
j
=
1
2
(
2
j
-
1
)
∑
k
=
0
j
(
j
k
)
i
j
-
k
∑
l
=
0
k
-
1
(
k
l
)
M
0
l
(
∑
i
=
0
N
-
1
a
i
i
k
-
1
)
,
where a
i are Daubechies wavelet coefficients. Finally, we get the system for the calculation of connection coefficients. Consider
(21)
(
T
-
1
2
d
-
1
I
m
d
)
⋀
0
,
d
=
(
0
d
!
2
-
j
/
2
)
,
where m
d is a row vector with all the m
i
d.


## 4. Implementation of Wavelet-Galerkin Method

Consider the following form of boundary value problem:(22)
y
′′
(
x
)
+
a
y
′
(
x
)
+
b
y
(
x
)
=
f
(
x
)
,
y
(
0
)
=
α
,
y
(
1
)
=
β
,
where a, b, α, and β are real constants. A trial solution for (22) is
(23)
y
(
x
)
≅
∑
l
c
l
2
j
/
2
ϕ
(
2
j
x
-
l
)
,
l
∈
Z
.Use (23) in (22) to obtain
(24)
d
2
d
x
2
∑
l
c
l
2
j
/
2
ϕ
(
2
j
x
-
l
)
+
a
d
d
x
∑
l
c
l
2
j
/
2
ϕ
(
2
j
x
-
l
)
+
b
∑
l
c
l
2
j
/
2
ϕ
(
2
j
x
-
l
)
-
f
(
x
)
≠
0
,
l
∈
Z
.For simplicity, use substitutiony
=
2
j
x
,
C
l
=
2
j
/
2
c
l. Also, we have d
n
/
d
x
n
=
2
n
j
d
n
/
d
y
n
,
n
∈
N.Now (24) implies
(25)
∑
l
C
l
2
2
j
d
2
d
y
2
ϕ
(
y
-
l
)
+
a
∑
l
C
l
2
j
d
d
y
ϕ
(
y
-
l
)
+
b
∑
l
C
l
ϕ
(
y
-
l
)
-
f
(
y
2
j
)
≠
0
,
l
∈
Z
.Multiplyingϕ
p
(
y
), on both sides of (25), and integrating, we get
(26)
∑
l
C
l
2
2
j
∫
ϕ
(
y
-
p
)
d
2
d
y
2
ϕ
(
y
-
l
)
d
y
+
a
∑
l
C
l
2
j
∫
ϕ
(
y
-
p
)
d
d
y
ϕ
(
y
-
l
)
d
y
+
b
∑
l
C
l
∫
ϕ
(
y
-
p
)
ϕ
(
y
-
l
)
d
y
=
∫
ϕ
(
y
-
p
)
f
(
y
2
j
)
d
y
,
where f
(
x
)
=
∑
i
=
0
m
b
i
x
i, is a polynomial of degree m in x.The orthonormality of Daubechies wavelets implies(27)
∫
ϕ
(
y
-
p
)
ϕ
(
y
-
l
)
d
y
=
δ
p
,
l
.Now use (11) and (27) in (26); we have
(28)
∑
l
C
l
(
Ω
l
-
p
0,2
+
a
Ω
l
-
p
0,1
+
b
δ
p
,
l
)
=
∑
i
=
0
m
b
i
∫
y
i
2
j
i
ϕ
(
y
-
p
)
d
y
,
or
(29)
∑
l
C
l
(
Ω
l
-
p
0,2
+
a
Ω
l
-
p
0,1
+
b
δ
p
,
l
)
=
∑
i
=
0
m
b
i
2
j
i
M
p
i
.Treatment of the boundary conditions [25] is as follows.Conditionsy
(
0
)
=
α and y
(
1
)
=
β imply
(30)
∑
l
C
l
ϕ
l
(
0
)
=
α
,
∑
l
C
l
ϕ
l
(
1
)
=
β
.Since inner product ofϕ
l
(
y
) and ϕ
p
(
y
) gives (27), therefore (30) implies
(31)
∑
l
C
l
δ
l
,
p
(
0
)
=
α
,
∑
l
C
l
δ
l
,
p
(
1
)
=
β
.


## 5. Quasilinearization

The quasilinearization approach is a generalized Newton-Raphson technique for functional equations [26, 27]. It converges quadratically to the exact solution, if there is convergence at all, and it has monotone convergence.Let us consider the nonlinearnth order differential equation as follows:
(32)
y
′′
(
x
)
=
f
(
y
(
x
)
,
x
)
,
with the boundary conditions
(33)
y
(
a
)
=
α
,
y
(
b
)
=
β
,
a
≤
x
≤
b
.Here,f is a function of the function y
(
x
). Choose an initial approximation of the function y
(
x
), let us say y
0
(
x
); it may be y
0
(
x
)
=
α, for a
≤
x
≤
b. The function f can now be expanded around the function y
0
(
x
) by the use of the Taylor series. Consider
(34)
f
(
y
(
x
)
,
x
)
≈
f
(
y
0
(
x
)
,
x
)
+
(
y
(
x
)
-
y
0
(
x
)
)
f
y
0
(
y
0
(
x
)
,
x
)
,
where second and higher order terms are ignored. Using (34) in (32), we get
(35)
y
′′
(
x
)
≈
f
(
y
0
(
x
)
,
x
)
+
(
y
(
x
)
-
y
0
(
x
)
)
f
y
0
(
y
0
(
x
)
,
x
)
.Solving (35) for y
(
x
), call it y
1
(
x
) and expand (32) about y
1
(
x
)
(36)
y
′′
(
x
)
≈
f
(
y
1
(
x
)
,
x
)
+
(
y
(
x
)
-
y
1
(
x
)
)
f
y
1
(
y
1
(
x
)
,
x
)
,we obtain a third approximation fory
(
x
), call it y
2
(
x
). Assume that the problem converges and continue the procedure for obtaining desired accuracy. Recurrence relation is of the form
(37)
y
r
+
1
′′
(
x
)
=
f
(
y
r
(
x
)
,
x
)
+
(
y
r
+
1
(
x
)
-
y
r
(
x
)
)
f
y
r
(
y
r
(
x
)
,
x
)
,
where y
r
(
x
) is known and can be used for obtaining y
r
+
1
(
x
). Equation (37) is always a linear differential equation and boundary conditions are
(38)
y
r
+
1
(
a
)
=
α
,
y
r
+
1
(
b
)
=
β
.Now consider the nonlinear second order differential equation of the form [26]
(39)
y
′′
(
x
)
=
f
(
y
′
(
x
)
,
y
(
x
)
,
x
)
.Here, the first derivativey
′
(
x
) can be considered as another function and (39) implies
(40)
y
r
+
1
′′
(
x
)
=
f
(
y
r
′
(
x
)
,
y
r
(
x
)
,
x
)
+
(
y
r
+
1
(
x
)
-
y
r
(
x
)
)
f
y
r
(
y
r
′
(
x
)
,
y
r
(
x
)
,
x
)
+
(
y
r
+
1
′
(
x
)
-
y
r
′
(
x
)
)
f
y
r
′
(
y
r
′
(
x
)
,
y
r
(
x
)
,
x
)
,
with the same boundary conditions
(41)
y
r
+
1
(
a
)
=
α
,
y
r
+
1
(
b
)
=
β
.Similarly, one can follow the same procedure for higher order nonlinear differential equations to obtain the recurrence relation(42)
L
n
y
r
+
1
(
x
)
=
f
(
y
r
(
x
)
,
y
r
′
(
x
)
,
…
,
y
r
n
-
1
(
x
)
,
x
)
+
∑
j
=
0
n
-
1
‍
(
y
r
+
1
j
(
x
)
-
y
r
j
(
x
)
)
×
f
y
j
(
y
r
(
x
)
,
y
r
′
(
x
)
,
…
,
y
r
n
-
1
(
x
)
,
x
)
,
where n is order of the differential equation. Equation (42) is always a linear differential equation and can be solved recursively, where y
r
(
x
) is known and one can use it to get y
r
+
1
(
x
).In order to test the wavelet-Galerkin method with quasilinearization technique, four different nonlinear problems are considered.


## 6. Applications

In this section, we solve some nonlinear ordinary differential equations by the wavelet-Galerkin method along with quasilinearization technique and compare the results with those obtained by other methods and exact solution.Example 1.
We consider the nonlinear two-point boundary value problem [28] as follows:
(43)
y
′′
(
x
)
+
y
2
(
x
)
=
x
4
+
2
,
0
<
x
<
1
,
subject to the boundary conditions y
(
0
)
=
0
,
y
(
1
)
=
1. The exact solution is y
(
x
)
=
x
2. Apply quasilinearization to (43); we get
(44)
y
r
+
1
′′
(
x
)
+
2
y
r
(
x
)
y
r
+
1
(
x
)
=
y
r
2
(
x
)
+
2
+
x
4
,
with the boundary conditions y
r
+
1
(
0
)
=
0, y
r
+
1
(
1
)
=
1.
Applying wavelet-Galerkin method to (44), we have
(45)
∑
l
C
l
(
Ω
l
-
p
0,2
+
2
y
r
(
x
)
δ
p
,
l
)
=
y
r
2
(
x
)
+
∑
i
=
0
4
b
i
2
j
i
M
p
i
,
and boundary conditions imply
(46)
∑
l
C
l
δ
l
,
p
(
0
)
=
0
,
∑
l
C
l
δ
l
,
p
(
1
)
=
1
,
where b
=
[
2,0
,
0,0
,
1
], with the initial approximation y
0
(
x
)
=
0.
We solved (43) by using D
16 and fixed the level of resolution J
=
13. Figure 1 shows the exact solution and approximate solution by proposed method at first and second iteration; that is, y
0
(
x
) is the initial approximation, and by using y
0
(
x
) we get y
1
(
x
), that is, the solution at first iteration, and then y
1
(
x
) is used to get y
2
(
x
), which is the solution of (43) at second iteration. Table 1 is used to compare the approximate solution by proposed method at second iteration for J
=
13 and exact solution. We may get more accurate results while increasing level of resolution at higher iteration.Table 1
Comparison of exact solutiony
exact and solution by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration: level of resolution J
=
13 and D
16 is used as Galerkin bases.
Third  iteration
J
=
13
x
y
NEW
y
exact
Absolute error
0.1
1.0022
e
-
2
1.0
e
-
2
2.2423
e
-
5
0.2
4.0045
e
-
2
4.0
e
-
2
4.4852
e
-
5
0.3
9.0067
e
-
2
9.0
e
-
2
6.7318
e
-
5
0.4
1.6009
e
-
1
1.6
e
-
1
8.9901
e
-
5
0.5
2.5011
e
-
1
2.5
e
-
1
1.1275
e
-
4
0.6
3.6014
e
-
1
3.6
e
-
1
1.3613
e
-
4
0.7
4.9016
e
-
1
4.9
e
-
1
1.6038
e
-
4
0.8
6.4019
e
-
1
6.4
e
-
1
1.8602
e
-
4
0.9
8.1021
e
-
1
8.1
e
-
1
2.1369
e
-
4Comparison of exact solution and solution by wavelet-Galerkin method with quasilinearization technique atJ
=
13, for different iterations, and we used D
16.
(a)
(b)
(c)
(d)Example 2.
We consider the nonlinear Bratu’s boundary value problem as follows:(47)
-
y
′′
(
x
)
=
λ
e
y
(
x
)
,
y
(
0
)
=
0
,
y
(
1
)
=
0
.
The quasilinear form of (47) is
(48)
y
r
+
1
′′
(
x
)
+
λ
e
y
r
(
x
)
y
r
+
1
(
x
)
=
-
λ
e
y
r
(
x
)
(
1
-
y
r
(
x
)
)
,
with the boundary conditions y
r
+
1
(
0
)
=
0, y
r
+
1
(
1
)
=
0.
Wavelet-Galerkin method for (48) implies
(49)
∑
l
C
l
(
Ω
l
-
p
0,2
+
λ
e
y
r
(
x
)
δ
p
,
l
)
=
-
λ
e
y
r
(
x
)
(
1
-
y
r
(
x
)
)
,
and, from boundary conditions, we have
(50)
∑
l
C
l
δ
l
,
p
(
0
)
=
0
,
∑
l
C
l
δ
l
,
p
(
1
)
=
0
,
with the initial approximation y
0
(
x
)
=
0.
Bratu’s boundary value problem is solved by usingD
12 as Galerkin bases and at level of resolution J
=
13. Decomposition method [17] is already implemented on (47). We compared our results with the results obtained by decomposition method and exact solution. Our results are more accurate as compared to decomposition method [17] at λ
=
1 and λ
=
2 as shown in Tables 2 and 3, respectively. We used the MATLAB command of one-dimensional data interpolation using spline to get the values at x
=
0.1
,
x
=
0.2
,
…
,
x
=
0.9 and plot the exact and approximate solutions at these points for λ
=
1 and λ
=
2 as shown in Figure 2.Table 2
Comparison of exact solutiony
exact, solution by decomposition method y
DM, and solution by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration and D
12 is used as Galerkin bases.
λ
=
1
Second iteration
J
=
13
x
y
DM [17]
y
NEW
y
exact
0.1
0.0471616875
0.0498386633
0.0498467900
0.2
0.0871680000
0.0891737644
0.0891899350
0.3
0.1177614375
0.1175850602
0.1176090956
0.4
0.1369920000
0.1347586231
0.1347902526
0.5
0.1435546875
0.1405003497
0.1405392142
0.6
0.1369920000
0.1347446015
0.1347902526
0.7
0.1177614375
0.1175571769
0.1176090956
0.8
0.0871680000
0.0891323318
0.0891899350
0.9
0.0471616875
0.0497841328
0.0498467900Table 3
Comparison of exact solutiony
exact, solution by decomposition method y
DM, and solution by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration and D
12 is used as Galerkin bases.
λ
=
2
Second iteration
J
=
13
x
y
DM [17]
y
NEW
y
exact
0.1
0.1144107440
0.1143861795
0.0991935000
0.2
0.2064191156
0.2063705472
0.1917440000
0.3
0.2738793116
0.2738079343
0.2679915000
0.4
0.3150893646
0.3149970528
0.3183360000
0.5
0.3289524214
0.3288416969
0.3359375000
0.6
0.3150893646
0.3149632892
0.3183360000
0.7
0.2738793116
0.2737413240
0.2679915000
0.8
0.2064191156
0.2062728271
0.1917440000
0.9
0.1144107440
0.1142597336
0.0991935000Comparison of exact solution and solution by wavelet-Galerkin method with quasilinearization technique at level of resolutionJ
=
13 and we use D
12 as Galerkin bases and at λ
=
1 and λ
=
2, respectively.
(a)
(b)Example 3.
Consider the nonlinear Troesch’s boundary value problem as follows:(51)
y
′′
(
x
)
-
λ
sinh
⁡
⁡
(
λ
y
(
x
)
)
=
0
,
0
≤
x
≤
1
,
y
(
0
)
=
0
,
y
(
1
)
=
1
.
The quasilinearized form of (51) is
(52)
y
r
+
1
′′
(
x
)
-
λ
2
cosh
⁡
⁡
(
λ
y
r
(
x
)
)
y
r
+
1
(
x
)
=
λ
sinh
⁡
⁡
(
λ
y
r
(
x
)
)
-
λ
2
y
r
(
x
)
cosh
⁡
⁡
(
λ
y
r
(
x
)
)
,
where 0
≤
x
≤
1, with the boundary conditions y
r
+
1
(
0
)
=
0, y
r
+
1
(
1
)
=
1.
Implementation of wavelet-Galerkin method to (52) implies
(53)
∑
l
C
l
(
Ω
l
-
p
0,2
-
λ
2
cosh
⁡
⁡
(
λ
y
r
(
x
)
)
δ
p
,
l
)
=
λ
sinh
⁡
⁡
(
λ
y
r
(
x
)
)
-
λ
2
y
r
(
x
)
cosh
⁡
⁡
(
λ
y
r
(
x
)
)
,
and boundary conditions lead to
(54)
∑
l
C
l
δ
l
,
p
(
0
)
=
0
,
∑
l
C
l
δ
l
,
p
(
1
)
=
1
,
with the initial approximation y
0
(
x
)
=
0. Tables 4 and 5 and Figure 3 represent the solution of (51) at second iteration. We use D
8 as Galerkin bases to find the solution of (51) at different level of resolutions and at λ
=
0.5 and λ
=
1 as shown in Tables 4 and 5, respectively. Solutions by proposed method are compared with variational iteration method [18] and with exact solution. Our results are in high agreement with exact solution and better than variational iteration method [18].Table 4
Comparison of exact solutiony
exact, solution by variational iteration method y
VIM, and solutions by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration, at different level of resolutions, and we used D
8 as Galerkin bases.
λ
=
0.5
Second   iteration
J
=
11
J
=
13
x
y
VIM [18]
y
NEW
y
NEW
y
exact
0.1
0.100042
0.095994
0.095956
0.095177
0.2
0.200334
0.192228
0.192152
0.190634
0.3
0.301128
0.288944
0.288829
0.286653
0.4
0.402677
0.386385
0.386231
0.383523
0.5
0.505241
0.484798
0.484606
0.481537
0.6
0.609082
0.584436
0.584204
0.581002
0.7
0.714470
0.685556
0.685284
0.682235
0.8
0.821682
0.788424
0.788112
0.785572
0.9
0.931008
0.893316
0.892962
0.891367
1.0
1.042740
1.000592
1.000194
1.000000Table 5
Comparison of exact solutiony
exact, solution by variational iteration method y
VIM, and solutions by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration, at different level of resolutions, and we used D
8 as Galerkin bases.
λ
=
1
Seconditeration
J
=
11
J
=
13
x
y
VIM [18]
y
NEW
y
NEW
y
exact
0.1
0.100167
0.084715
0.084674
0.081797
0.2
0.201339
0.170280
0.170197
0.164531
0.3
0.304541
0.257558
0.257433
0.249167
0.4
0.410841
0.347444
0.347275
0.336732
0.5
0.521373
0.440880
0.440666
0.428347
0.6
0.637362
0.538877
0.538616
0.525274
0.7
0.760162
0.642538
0.642226
0.628971
0.8
0.891287
0.753089
0.752722
0.741168
0.9
1.032460
0.871929
0.871495
0.863970
1.0
1.185650
1.000718
1.000226
1.000000Comparison of exact solution and solution by wavelet-Galerkin method with quasilinearization technique at level of resolutionJ
=
13 and D
8 is used as Galerkin bases and at λ
=
0.5 and λ
=
1, respectively.
(a)
(b)Example 4.
Consider the following nonlinear boundary value problem:(55)
y
′′
(
x
)
+
y
′
(
x
)
+
y
3
(
x
)
+
y
(
x
)
=
f
(
x
)
,
subject to the boundary conditions y
(
0
)
=
0
,
y
(
1
)
=
0, where f
(
x
)
=
2
+
2
x
+
x
2
-
20
x
3
-
5
x
4
-
x
5
+
x
6
-
3
x
9
+
3
x
12
-
x
15 and exact solution is given by
(56)
y
(
x
)
=
x
2
-
x
5
.
Applying the quasilinearization technique to (55), we get
(57)
y
r
+
1
′′
(
x
)
+
y
r
+
1
′
(
x
)
+
(
1
+
3
y
r
2
(
x
)
)
y
r
+
1
(
x
)
=
f
(
x
)
+
2
y
r
3
(
x
)
,
with the boundary conditions y
r
+
1
(
0
)
=
0
,
y
r
+
1
(
1
)
=
0.
Now we apply the wavelet-Galerkin method to (57) as follows:
(58)
∑
l
C
l
(
Ω
l
-
p
0,2
+
Ω
l
-
p
0,1
+
(
1
+
3
y
r
2
(
x
)
)
δ
p
,
l
)
=
2
y
r
3
(
x
)
+
∑
i
=
0
15
b
i
2
j
i
M
p
i
,
and boundary conditions imply
(59)
∑
l
C
l
δ
l
,
p
(
0
)
=
0
,
∑
l
C
l
δ
l
,
p
(
1
)
=
0
,
where b
=
[
2,2
,
1
,
-
20
,
-
5
,
-
1,1
,
0,0
,
-
3,0
,
0,3
,
0,0
,
-
1
], with the initial approximations y
0
(
x
)
=
0.
Numerical solution by proposed method for (55) at second iteration and level of resolution J
=
13 are stable and accurate as shown in Table 6. To get the more accurate results increase the iteration or level of resolution or both. Figure 4 shows that approximate solution converges to the exact solution while increasing iterations. Here, D
10 is used as Galerkin bases.Table 6
Comparison of exact solutiony
exact and solution by the wavelet-Galerkin method with quasilinearization technique y
NEW at second iteration: level of resolutions J
=
13 and D
10 is used as Galerkin bases.
Second   iteration
J
=
13
x
y
NEW
y
exact
Absolute error
0.1
0.0100320512
0.0099900000
4.2051152889
e
-
005
0.2
0.0397617227
0.0396800000
8.1722744304
e
-
005
0.3
0.0876864075
0.0875700000
1.1640750572
e
-
004
0.4
0.1499017545
0.1497600000
1.4175454098
e
-
004
0.5
0.2189016605
0.2187500000
1.5166050983
e
-
004
0.6
0.2823782628
0.2822400000
1.3826276219
e
-
004
0.7
0.3220219443
0.3219300000
9.1944320423
e
-
005
0.8
0.3123213602
0.3123199999
1.3602317094
e
-
006
0.9
0.2193634808
0.2195099999
1.4651915492
e
-
004Comparison of exact solution and solution by wavelet-Galerkin method with quasilinearization technique atJ
=
13, for different iterations, and we use D
10 as Galerkin bases.
(a)
(b)
(c)
(d)


## 7. Conclusion

It is shown that the proposed method, wavelet-Galerkin method with quasilinearization technique, gives stable and accurate results when applied to different nonlinear boundary value problems. The proposed method provides better and more accurate results as compared to variational iteration method and decomposition method, as shown in Tables2–5. Also results are in good agreement with exact solutions. Figures 1 and 4 show that approximate solution converges to the exact solution while iterations are increased and absolute error goes down. The main advantage of the proposed method is that the different type of nonlinearities can be easily handled.



---

*Source: 868934-2014-06-16.xml*