Research

# On singular nonlinear distributional and impulsive initial and boundary value problems

Seppo Heikkilä

Author Affiliations

Department of Mathematical Sciences, University of Oulu, BOX 3000, FIN-90014, Oulu, Finland

Boundary Value Problems 2011, 2011:24  doi:10.1186/1687-2770-2011-24

 Received: 29 April 2011 Accepted: 16 September 2011 Published: 16 September 2011

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

#### Purpose

To derive existence and comparison results for extremal solutions of nonlinear singular distributional initial value problems and boundary value problems.

#### Main methods

Fixed point results in ordered function spaces and recently introduced concepts of regulated and continuous primitive integrals of distributions. Maple programming is used to determine solutions of examples.

#### Results

New existence results are derived for the smallest and greatest solutions of considered problems. Novel results are derived for the dependence of solutions on the data. The obtained results are applied to impulsive differential equations. Concrete examples are presented and solved to illustrate the obtained results.

MSC: 26A24, 26A39, 26A48, 34A12, 34A36, 37A37, 39B12, 39B22, 47B38, 47J25, 47H07, 47H10, 58D25

##### Keywords:
distribution; primitive; integral; regulated; continuous; initial value problem; boundary value problem; singular; distributional

### 1 Introduction

In this paper, existence and comparison results are derived for the smallest and greatest solutions of first and second order singular nonlinear initial value problems as well as second order boundary value problems.

Recently, similar problems are studied in ordered Banach spaces, e.g., in [1-4], by converting problems into systems of integral equations, integrals in these systems being Bochner-Lebesgue or Henstock-Kurzweil integrals. A novel feature in the present study is that the right-hand sides of the considered differential equations comprise distributions on a compact real interval [a, b]. Every distribution is assumed to have a primitive in the space R [ a , b ] of those functions from [a, b] to ℝ which are left-continuous on (a, b], right-continuous at a, and which have right limits at every point of (a, b). With this presupposition, the considered problems can be transformed into integral equations which include the regulated primitive integral of distributions introduced recently in [5].

The paper is organized as follows. Distributions on [a, b], their primitives, regulated primitive integrals and some of their properties, as well as a fixed point lemma are presented in Section 2. In Section 3, existence and comparison results are derived for the smallest and greatest solutions of first order initial value problems.

A fact that makes the solution space R [ a , b ] important in applications is that it contains primitives of Dirac delta distributions δλ, λ ∈ (a, b). This fact is exploited in Section 4, where results of Section 3 are applied to impulsive differential equations. The continuous primitive integral of distributions introduced in [6] is also used in these applications.

Existence of the smallest and greatest solutions of the second order initial and boundary value problems, and dependence of these solutions on the data are studied in Sections 5 and 6. Applications to impulsive problems are also presented.

Considered differential equations may be singular, distributional and impulsive. Differential equations, initial and boundary conditions and impulses may depend functionally on the unknown function and/or on its derivatives, and may contain discontinuous nonlinearities. Main tools are fixed point theorems in ordered spaces proved in [7] by generalized monotone iteration methods. Concrete problems are solved to illustrate obtained results. Iteration methods and Maple programming are used to determine solutions.

### 2 Preliminaries

Distributions on a compact real interval [a, b] are (cf. [8]) continuous linear functionals on the topological vector space D of functions φ : ℝ → ℝ possessing for every j ∈ ℕ0 a continuous derivative φ(j) of order j that vanishes on ℝ\(a, b). The space D is endowed with the topology in which the sequence (φk) of D converges to φ D if and only if φ k ( j ) φ ( j ) uniformly on (a, b) as k → ∞ and j ∈ ℕ0. As for the theory of distributions, see, e.g., [9,10].

In this paper, every distribution g on [a, b] is assumed to have a primitive, i.e., a function G R [ a , b ] whose distributional derivative G' equals to g, in the function space

R [ a , b ] = { G : lim t s + G ( s ) e x i s t s , lim s t - G ( s ) = G ( t ) i f a s < t b , a n d G ( a ) : = lim s a + G ( s ) } . (2.1)

The value 〈g, φ〉 of g at φ D is thus given by

g , φ = G , φ = - G , φ = - a b G ( t ) φ ( t ) d t .

Such a distribution g is called RP integrable. Its regulated primitive integral is defined by

r s t g : = G ( t ) G ( s ) , a s t b , where G is a primitive of g in [ a , b ] . (2.2)

As noticed in [5], the regulated primitive integral generalizes the wide Denjoy integral, and hence also Riemann, Lebesgue, Denjoy and Henstock-Kurzweil integrals.

Denote by A R [ a , b ] the set of those distributions on [a, b] that are RP integrable on [a, b]. If g A R [ a , b ] , then the function t r a t g is that primitive of g which belongs to the set

P R [ a , b ] = { G R [ a , b ] : G ( a ) = 0 } .

It can be shown (cf. [5]) that a relation ≼, defined by

f g i n A R [ a , b ] i f a n d o n l y i f r a t f r a t g f o r a l l t [ a , b ] , (2.3)

is a partial ordering on A R [ a , b ] . In particular,

f = g i n A R [ a , b ] i f a n d o n l y i f r a t f = r a t g f o r a l l t [ a , b ] . (2.4)

Given partially ordered sets X = (X, ≤) and Y = (Y, ≼), we say that a mapping f : X Y is increasing if f(x) ≼ f(y) whenever x y in X, and order-bounded if there exist f± Y such that f- f (x) ≼ f+ for all x X.

The following fixed point result is a consequence of [11], Theorem A.2.1, or [7], Theorem 1.2.1 and Proposition 1.2.1.

Lemma 2.1. Given a partially ordered set P = (P, ≤), and its order interval [x-, x+] = {x P : x- x x+}, assume that a mapping G : [x-, x+] → [x-, x+] is increasing, and that each well-ordered chain of the range G[x-, x+] of G has a supremum in P and each inversely well-ordered chain of G[x-, x+] has an infimum in P. Then G has the smallest and greatest fixed points, and they are increasing with respect to G.

Remarks 2.1. Under the hypotheses of Lemma 2.1, the smallest fixed point x* of G is by [[7], Theorem 1.2.1] the maximum of the chain C of [x-, x+] that is well ordered, i.e., every nonempty subset of C has the smallest element, and that satisfies

( I ) x - = min C , a n d i f x - < x , t h e n x C i f a n d o n l y i f x = sup G [ { y C : y < x } ] .

The smallest elements of C are Gn(x-), n ∈ ℕ0, as long as Gn(x-) = G(Gn-1(x-)) is defined and Gn-1(x-) < Gn(x-), n ∈ ℕ. If Gn-1(x-) = Gn (x-) for some n ∈ ℕ, there is the smallest such an n, and x* = Gn-1(x-) is the smallest fixed point of G in [x-, x+]. If x ω = sup n G n ( x - ) is defined in P and is a strict upper bound of {Gn(x-)}n∈ℕ, then xω is the next element of C. If xω = G(xω), then x* = xω, otherwise the next elements of C are of the form Gn(xω), n ∈ ℕ, and so on.

The greatest fixed point x* of G is the minimum of the chain D of [x-, x+] that is inversely well ordered, i.e., every nonempty subset of D has the greatest element, and that has the following property:

( I I ) x + = max D , a n d i f x < x + , t h e n x D i f a n d o n l y i f x = inf G [ { y D : x < y } ] .

The greatest elements of D are n-fold iterates Gn(x+), as long as they are defined and Gn(x+) < Gn-1(x+). If equality holds for some n ∈ ℕ, then x* = Gn-1(x+) is the greatest fixed point of G in [x-, x+].

### 3 First order initial value problems

In this section, existence and comparison results are derived for the smallest and greatest solutions of first order initial value problems. Denote by L l o c 1 ( a , b ] , -∞< a < b < ∞, the space of locally Lebesgue integrable functions from the half-open interval (a, b] to ℝ. L l o c 1 ( a , b ] is ordered a.e. pointwise, and its a.e. equal functions are identified.

Given p : [a, b] → ℝ+, consider the initial value problem (IVP)

( p u ) = g ( u ) , lim t a + ( p u ) ( t ) = c ( u ) , (3.1)

where c(u) ∈ ℝ, and g ( u ) A R [ a , b ] . We are looking for solutions of (3.1) from the set

S = { u L l o c 1 ( a , b ] : lim t a + ( p u ) ( t ) e x i s t s , a n d p u R [ a , b ] } , (3.2)

Definition 3.1. We say that a function u S is a subsolution of the IVP (3.1) if

( p u ) g ( u ) , a n d lim t a + ( p u ) ( t ) c ( u ) . (3.3)

If reversed inequalities hold in (3.3), we say that u is a supersolution of (3.1). If equalities hold in (3.3), then u is called a solution of (3.1).

We shall first transform the IVP (3.1) into an integral equation.

Lemma 3.1. Given c(u) ∈ ℝ, u L l o c 1 ( a , b ] and p : [a, b] → ℝ+, assume that 1 p L l o c 1 ( a , b ] , and that g ( u ) A R [ a , b ] . Then u is a solution of the IVP (3.1) in S if and only if u is a solution of the following integral equation:

u ( t ) = 1 p ( t ) c ( u ) + r a t g ( u ) , t ( a , b ] . (3.4)

Proof: Assume that u is a solution of (3.1) in S. The definition of S and (3.1) ensure by (2.2) that

r r t g ( u ) = r r t ( p u ) = ( p u ) ( t ) ( p u ) ( r ) , a < r t < b .

Allowing r tend to a+ and applying the initial condition of (3.1) we see that (3.4) is valid. Conversely, let u be a solution of (3.4). According to (3.4) we have

( p u ) ( t ) = c ( u ) + r a t g ( u ) , t ( a , b ] . (3.5)

This equation implies that u S, that the initial condition of (3.1) is valid, and that

( p u ) = g ( u ) .

Thus, u is a solution of the IVP (3.1) in S. □

Our first existence and comparison result for the IVP (3.1) reads as follows.

Theorem 3.1. Assume that g : L l o c 1 ( a , b ] A R [ a , b ] is increasing, that p : [a, b] → ℝ+, that 1 p L l o c 1 ( a , b ] , and that the IVP (3.1) has a subsolution u- and a supersolution u+ in S satisfying u- u+. Then (3.1) has the smallest and greatest solutions within the order interval [u-, u+] of S. Moreover, these solutions are increasing with respect to g and c.

Proof: Define a mapping G : L l o c 1 ( a , b ] L l o c 1 ( a , b ] by

G ( u ) ( t ) : = 1 p ( t ) c ( u ) + r a t g ( u ) , t ( a , b ] . (3.6)

Because g is increasing, it follows from (2.3) and (3.6) that G is increasing. Applying (2.3), [[5], Theorem 7] and Definition 3.1 we see that if u L l o c 1 ( a , b ] and u- u u+, then

( p u - ) ( t ) - c ( u ) ( p u - ) ( t ) - lim r a + ( p u - ) ( r ) = lim r a + r r t ( p u - ) (1) = r a t ( p u - ) r a t g ( u - ) r a t g ( u ) , a < t b . (2) (3)

Thus

u - ( t ) 1 p ( t ) c ( u ) + r a t g ( u ) = G ( u ) ( t ) , t ( a , b ] .

Similarly, it can be shown that G(u)(t) ≤ u+(t) for each t ∈ (a, b]. Thus, G maps the order interval [u-, u+] of L l o c 1 ( a , b ] into [u-, u+]. Let W be a well-ordered or an inversely well-ordered chain in G[u-, u+]. It follows from [[1], Proposition 9.36] and its dual that sup W and inf W exist in L l o c 1 ( a , b ] .

The above proof shows that the operator G defined by (3.6) satisfies the hypotheses of Lemma 2.1 when P = L l o c 1 ( a , b ] . Thus G has the smallest fixed point u* and the greatest fixed point u* in [u-, u+]. These fixed points are the smallest and greatest solutions of the integral equation (3.4) in [u-, u+]. This result and Lemma 3.1 imply that u* and u* belong to S, and they are the smallest and greatest solutions of the IVP (3.1) in [u-, u+]. Moreover, u* and u* are by Lemma 2.1 increasing with respect to G. This result implies by (2.3) and (3.6) the last conclusion of Theorem. □

The following result is a consequence of Theorem 3.1.

Proposition 3.1. Assume that mappings g : L l o c 1 ( a , b ] A R [ a , b ] and c : L l o c 1 ( a , b ] are increasing and order-bounded, that p : [a, b] → ℝ+, and that 1 p L l o c 1 ( a , b ] . Then, the IVP (3.1) has in S the smallest and greatest solutions that are increasing with respect to g and c.

Proof: Because g and c are order-bounded, there exist g A R [ a , b ] and c± ∈ ℝ such that g-g(x) ≼ g+ and c- c(x) ≤ c+ for all x L l o c 1 ( a , b ] . Denote

u ( t ) = 1 p ( t ) c + r a t g , t ( a , b ] .

Then u± S, and

( p u - ) = g - g ( x ) g + = ( p u + ) f o r a l l x L l o c 1 ( a , b ] ,

and

lim t a + ( p u - ) ( t ) = c - c ( x ) c + = lim t a + ( p u + ) f o r a l l x L l o c 1 ( a , b ] .

Thus u- is a subsolution and u+ is a supersolution of (3.1), whence the IVP (3.1) has by Theorem 3.1 the smallest solution u* and the greatest solution u* in the order interval [u-, u+] of S.

If u S is any solution of (3.1), then

1 p ( t ) c - + r a t g - 1 p ( t ) c ( u ) + r a t g ( u ) 1 p ( t ) c + + r a t g + , t ( a , b ] ,

or equivalently,

u - ( t ) u ( t ) u + ( t ) , t ( a , b ] .

Consequently, u ∈ [u-, u+], whence u* and u* are the smallest and greatest of all the solutions of (3.1) in S. □

In the next proposition, the Henstock-Kurzweil integral K can be replaced by any of the integrals called Riemann, Lebesgue, Denjoy and wide Denjoy integrals.

Proposition 3.2. Assume that g(x) is RP integrable on [a, b] for every x L l o c 1 ( a , b ] , and that

r a t g ( x ) = i = 1 n H i ( t ) K a t f i ( x ) + H 0 ( t ) , t [ a , b ] , (3.7)

Where H 0 P R [ a , b ] , and for each i = 1,..., n, Hi : [a, b] → [0, ∞) has right limits on [a, b), is left-continuous on (a, b], and fi : [a, b] → ℝ satisfies the following hypotheses.

(fi1) fi(x) is Henstock-Kurzweil integrable on [a, b] for every x L l o c 1 ( a , b ] .

(fi2) There exist Henstock-Kurzweil integrable functions

f
i, f ̄ i : [ a , b ] such that K a t f - i K a t f i ( x ) K a t f i ( y ) K a t f ̄ i , t [ a , b ] whenever x y in L l o c 1 ( a , b ] .

If c : L l o c 1 ( a , b ] is increasing and order-bounded, then the IVP (3.1) has in S the smallest and greatest solutions that are increasing with respect to fi and c.

Proof: The hypotheses imposed above ensure by (2.3) and (3.7) that g is an increasing mapping from L l o c 1 ( a , b ] to the order interval [g-, g+] of A R [ a , b ] , where

r a t g = i = 1 n H i ( t ) K a t f ¯ i + H 0 ( t ) , r a t g + = i = 1 n H i ( t ) K a t f ¯ i + H 0 ( t ) , t [ a , b ] .

Thus the conclusions follow from Proposition 3.1.

Example 3.1. Assume that

r a t g ( x ) = H 1 ( t ) K a t f 1 ( x ) + H 0 ( t ) , t [ 0 , b ] , (3.8)

where b ≥ 1, H 0 R [ 0 , b ] , H1 is the Heaviside step function, i.e., H 1 ( t ) = 1 , t > 0 , 0 , t 0 ,

f 1 ( x ) ( t ) = 1 1 0 5 1 0 5 arctan 1 2 1 x ( t ) - H 0 ( t ) p ( t ) d t sin 1 t - 1 t sgn sin 1 t cos 1 t , t ( 0 , b ] , x L l o c 1 ( 0 , b ] , [ z ] = max { n : n z } a n d s g n ( z ) = z | z | , z 0 , 0 , z = 0 .

Note, that the greatest integer function [·] occurs in the function f1(x). Prove that the IVP

( p u ) = g ( u ) , lim t 0 + ( p u ) ( t ) = 0 , (3.9)

where p(t) = t, t ∈ [0, b], has the smallest and greatest solutions, and calculate them.

Solution: Problem (3.9) is of the form (3.1), where c(u) = a = 0 and p(t) ≡ t. The hypotheses (f11) and (f21) are valid when

K 0 t f - 1 = - 2 t | sin 1 t | + H 0 ( t ) , t ( 0 , b ] , K 0 t f ̄ 1 = 2 t | sin 1 t | + H 0 ( t ) , t ( 0 , b ] .

Thus the IVP (3.9) has by Proposition 3.2 the smallest and greatest solutions. They are the smallest and greatest fixed points of the mapping G defined by

G ( x ) ( t ) : = 1 t r 0 t g ( x ) , t ( 0 , b ] , x L l o c 1 ( 0 , b ] . (3.10)

G is an increasing mapping from L l o c 1 ( 0 , b ] , to its order interval [u-, u+], where

u ( t ) : = 2 | sin 1 t | + H 0 ( t ) t , t ( 0 , b ] .

Calculating the successive approximations Gn(u±) we see that G7(u±) = G8(u±). This means by Remark 2.1 that u* = G7(u-) and u* = G7(u+) are the smallest and greatest fixed points of G in [u-, u+]. According to the proof of Proposition 3.1, u* and u* are also the smallest and greatest solutions, of the initial value problem (3.9) in S. The exact expressions of u* and u* are:

u * ( t ) = - arctan 6 7 2 2 9 1 0 0 0 0 | sin 1 t | + H 0 ( t ) t , t ( 0 , b ] , u * ( t ) = arctan 1 6 8 0 7 2 5 0 0 | sin 1 t | + H 0 ( t ) t , t ( 0 , b ] .

### 4 Applications to impulsive problems

In this section, we assume that Λ is a well-ordered subset of (a, b). Let δλ, λ ∈ Λ, denote the translation of Dirac delta distribution for which r a t δ λ = H ( t λ ) , t a, where H is the Heaviside step function. Consider the singular distributional Cauchy problem

( p u ) = λ Λ I ( λ , u ) δ λ + f ( u ) , lim t a + ( p u ) = c ( u ) , (4.1)

where p : [a, b] → ℝ+ and 1 p L l o c 1 ( a , b ] . The values of f are distributions on [a, b], and the values of I are real numbers.

Definition 4.1. By a solution of (4.1), we mean such a function u S that satisfies (4.1), for which p · u is continuous on [a, b]\Λ, and has impulses

Δ ( p u ) ( λ ) : = ( p u ) ( λ + ) - ( p u ) ( λ ) = I ( λ , u ) , λ Λ .

In the study of (4.1), the regulated primitive integral is replaced by the continuous primitive integral presented in [6]. A distribution g on [a, b] is called distributionally Denjoy (DD) integrable on [a, b], denote g A C [ a , b ] , if g has a continuous primitive, i.e., g is a distributional derivative of a function G C[a, b]. The continuous primitive integral of g is defined by

c s t g = G ( t ) G ( s ) , a s t b .

A C [ a , b ] is a proper subset of A R [ a , b ] , and for every g A C [ a , b ] its continuous and regulated primitive integrals are equal. As shown in [6], A C [ a , b ] contains functions that are wide Denjoy integrable, and hence also Riemann, Lebesgue, Denjoy and Henstock-Kurzweil integrable on [a, b]. On the other hand, distributional derivatives of nowhere differentiable Weierstrass function and almost everywhere differentiable Cantor function are distributionally but not wide Denjoy integrable.

It can be shown (cf. [6]) that relation , defined by

f g i f a n d o n l y i f c a t f c a t g f o r a l l t [ a , b ] , (4.2)

is a partial ordering on A C [ a , b ] .

Transformation of the Cauchy problem (4.1) into an integral equation is presented in the following lemma.

Lemma 4.1. Assume that u S, that f ( u ) A C [ a , b ] , and that λ Λ | I ( λ , u ) | < . Then u is a solution of (4.1) if and only if

u ( t ) = 1 p ( t ) c ( u ) + λ Λ I ( λ , u ) H ( t - λ ) + c a t f ( u ) , t ( a , b ] . (4.3)

Proof: Assume first that u S satisfies (4.3). Because Λ is well-ordered, it follows that if λ ∈ Λ and λ < sup Λ, then H(t - λ) = 1 on (λ, S(λ)], where S(λ) = min{μ ∈ Λ : λ < μ}. This property implies that if the function v : (a, b] → ℝ is defined by

v ( t ) = 1 p ( t ) c ( u ) + λ Λ I ( λ , u ) H ( t - λ ) , t ( a , b ] , (4.4)

then the function p · v is constant on every interval (λ, S(λ)], Λ ∋ Λ < sup Λ, on [a, min Λ], and on (sup Λ, b] if sup Λ < b. In particular, p v R [ a , b ] , and the distributional derivative of p · v is

( p v ) = λ Λ I ( λ , u ) δ λ . (4.5)

Thus

( p u ) = ( p v ) + f ( u ) = λ Λ I ( λ , u ) δ λ + f ( u ) .

Since t c a t f ( u ) is continuous on [a, b], then p · u is continuous on [a, b]\Λ. Because

( p v ) ( t ) - ( p v ) ( λ ) = I ( λ , u ) H ( t - λ ) = I ( λ , u ) , λ Λ , t ( λ , S ( λ ) ] ,

then

Δ ( p u ) ( λ ) = ( p u ) ( λ + ) - ( p u ) ( λ ) = ( p v ) ( λ + ) - ( p v ) ( λ ) = I ( λ , u ) , λ Λ .

Moreover lim t a + ( p u ) ( t ) = c ( u ) , so that u is a solution of the IVP (4.1).

Assume next that u S is a solution of (4.1). Denoting

z ( t ) = u ( t ) - v ( t ) , t [ a , b ] ,

where v is defined by (4.4), it follows from (4.1) and (4.5) that

( p z ) = f ( u ) , lim t a + ( p z ) = 0 .

Because f(u) is DD integrable on [a, b], then

( p z ) ( t ) = c a t f ( u ) , t [ a , b ] .

Thus

( p u ) ( t ) = ( p z ) ( t ) + ( p v ) ( t ) = c ( u ) + λ Λ I ( λ , u ) H ( t - λ ) + c a t f ( u ) , t [ a , b ] ,

or equivalently, (4.3) holds. □

Noticing that the IVP (4.1) is a special case of the Cauchy problem (3.1), where

g ( u ) = λ Λ I ( λ , u ) δ λ + f ( u ) , (4.6)

the results of Section 3 can be applied to study the IVP (4.1). The following result is a consequence of Proposition 3.1.

Proposition 4.1. The distributional IVP (4.1) has the smallest and greatest solutions that are increasing with respect to f and c, if f : L l o c 1 ( a , b ] A C [ a , b ] and c : L l o c 1 ( a , b ] are increasing and order-bounded, if p : [a, b] → ℝ+, if 1 p L l o c 1 ( a , b ] , and if I : Λ × L l o c 1 ( a , b ] has the following properties.

(I) λ Λ | I ( λ , x ) | M < for all x L l o c 1 ( a , b ] , and x I(λ,x) is increasing for all λ ∈ Λ.

Proof: The given hypotheses imply that (4.6) defines a mapping g : L l o c 1 ( a , b ] A R [ a , b ] that is increasing and order-bounded. Thus, the IVP (3.1) has by Proposition 3.1 the smallest solution u* and the greatest solution u* in S, and they are increasing with respect to g and c. By Lemma 4.1, u* and u* are the smallest and greatest solutions of the IVP (4.1), and they are increasing with respect to f, and c, since g is increasing with respect to f. □

The initial value problem

d d t ( p ( t ) u ( t ) ) = q ( t , u ( t ) , u ) a . e . o n [ a , b ] , lim t a + ( p ( t ) u ( t ) ) = c ( u ) , (4.7)

combined with the impulsive property:

Δ ( p u ) ( λ ) = ( p u ) ( λ + ) - ( p u ) ( λ ) = I ( λ , u ) , λ Λ , (4.8)

form a special case of the IVP (4.1) when f is the Nemytskij operator associated with the function q : [ a , b ] × × L l o c 1 ( a , b ] by

f ( x ) : = q ( , x ( ) , x ) , x L l o c 1 ( a , b ] .

Considering distributions δλ as generalized functions t αδ (t - λ), t ∈ [a, b], we can rewrite the system (4.7), (4.8) as

d d t ( p ( t ) u ( t ) ) = λ Λ I ( λ , u ) δ ( t - λ ) + q ( t , u ( t ) , u ) a . e . o n [ a , b ] , lim t a + ( p ( t ) u ( t ) ) = c ( u ) . (4.9)

For instance, Proposition 4.1 implies the following result:

Corollary 4.1. The impulsive Cauchy problem (4.9) has the smallest and greatest solutions which are increasing with respect to q and c, if c : L l o c 1 ( a , b ] is increasing and order- bounded, and if the hypotheses (I) and the following hypotheses are valid.

(q0) q(·, x(·); x) is Henstock-Kurzweil integrable on [a, b] for every x L l o c 1 ( a , b ] .

(q1) K a t q ( s , x ( s ) , x ) d s K a t q ( s , y ( s ) , y ) d s for all t ∈ [a, b] whenever × y in L l o c 1 ( a , b ] .

(q2) There exist Henstock-Kurzweil integrable functions q± : [a, b] → ℝ such that K a t q - ( s ) d s K a t q ( s , x ( s ) , x ) d s K a t q + ( s ) d s for all x L l o c 1 ( a , b ] and t ∈ [a, b].

Example 4.1. Determine the smallest and greatest solutions of the IVP

t u ( t ) + u ( t ) = 1 1 0 4 [ 1 0 4 arctan ( u ( 1 ) ) ] δ t - 1 2 + q ( t , u ) a . e . o n [ 0 , 1 ] , lim t 0 + ( t u ( t ) ) = 0 , (4.10)

when q is defined by

{ q ( t , x ) = [ 10 4 tanh ( 2 5 1 ( x ( s ) d s ) ] 10 4 h ( t ) , x L l o c 1 ( 0 , 1 ] , t ( 0 , 1 ] , where h ( t ) = | cos ( 1 t ) | + 1 t sgn ( cos ( 1 t ) ) sin ( 1 t ) , t ( 0 , 1 ] . (4.11)

'[·]' denotes, as before, the greatest integer function, and 'sgn' the sign function.

Solution: The IVP (4.10) is a special case of (4.6), when a = 0, b = 1, c(u) = 0, p(t) = t, t ∈ [0, 1], I ( 1 2 , x ) = 1 1 0 4 [ 1 0 4 arctan ( x ( 1 ) ) ] , and Λ = { 1 2 } . The validity of the hypotheses of Corollary 4.1 is easy to verify. Thus, the IVP (4.10) has the smallest and greatest solutions. These solutions are the smallest and greatest fixed points of G : L l o c 1 ( 0 , 1 ] L l o c 1 ( 0 , 1 ] , defined by

G ( x ) ( t ) = 1 t 1 1 0 4 [ 1 0 4 arctan ( x ( 1 ) ) ] H t - 1 2 + K 0 t q ( t , x ) , x L l o c 1 ( 0 , 1 ] , t ( 0 , 1 ] . (4.12)

Calculating the successive approximations

y n + 1 = G ( y n ) , y 0 = x - a n d z n + 1 = G ( z n ) , z 0 = x + , w h e r e x ( t ) = 2 t H t - 1 2 + 1 2 K 0 t h ( s ) d s = 2 t H t - 1 2 1 t cos ( 1 t ) , t ( 0 , 1 ] ,

it turns out that ( y n ) n = 0 1 7 is strictly increasing, that ( z n ) n = 0 1 6 is strictly decreasing, that y17 = G(y17), and that z16 = G(z16). Thus u* = y17 and u* = z16 are by Remark 2.1 the smallest and greatest solutions of (4.1) with c(u) = 0. The exact formulas of u* and u* are

u * ( t ) = - 4 4 3 9 H ( t - 1 2 ) 5 0 0 0 t - 6 3 1 3 1 0 0 0 0 cos ( 1 t ) , t ( 0 , 1 ] , u * ( t ) = 2 2 1 9 H ( t - 1 2 ) 2 5 0 0 t + 6 3 1 1 1 0 0 0 0 cos ( 1 t ) , t ( 0 , 1 ] .

Remarks 4.1. The function (t, x) α q(t, x), defined in (4.11), has the following properties.

• It is Henstock-Kurzweil integrable, but it is not Lebesgue integrable with respect to the independent variable t if x ≠ 0, because h is not Lebesgue integrable on [0,1].

• Its dependence on the variables t and x is discontinuous, since the signum function sgn, the greatest integer function [·], and the function h are discontinuous.

• Its dependence on the unknown function x is nonlocal, since the integral of function x appears in the argument of the tanh-function.

• Its dependence on x is not monotone, since h attains positive and negative values in an infinite number of disjoint sets of positive measure. For instance, y*(t) > y*(t) for all t ∈ (0, 1], but the difference function t α q(t, y*) -q(t, y*) is neither nonnegative-valued nor Lebesgue integrable on [0, 1].

Notice also that in Example 4.1 dependence of the function I ( 1 2 , x ) = 1 1 0 4 [ 1 0 4 arctan ( x ( 1 ) ) ] on x is discontinuous.

### 5 Second order initial value problems

We shall study the second order initial value problem in this section

( p u ) = f ( u , u ) , lim t a + ( p u ) ( t ) = c ( u , u ) , lim t a + u ( t ) = d ( u , u ) , (5.1)

where f : L l o c 1 ( a , b ] 2 A R [ a , b ] , c , d : L l o c 1 ( a , b ] 2 , p : [a, b] → ℝ+, -∞ < a < b < ∞.

We are looking for the smallest and greatest solutions of (5.1) from the set

Y = { u : ( a , b ] : u L l o c 1 ( a , b ] , lim t a + u ( t ) a n d lim t a + ( p u ) ( t ) e x i s t , a n d p u R [ a , b ] } . (5.2)

The IVP (5.1) can be converted to a system of integral equations which does not contain derivatives.

Lemma 5.1. Assume that p : [a, b] →ℝ+, that 1 p L 1 [ a , b ] and that f ( u , v ) A R [ a , b ] for all u , v L l o c 1 ( a , b ] . Then u is a solution of the IVP (5.1) in Y if and only if (u, u') = (u, v), where ( u , v ) L l o c 1 ( a , b ] 2 is a solution of the system

u ( t ) = d ( u , v ) + a t v ( s ) d s , t ( a , b ] , v ( t ) = 1 p ( t ) c ( u , v ) + r a t f ( u , v ) , t ( a , b ] . (5.3)

Proof: Assume that u is a solution of the IVP (5.1) in Y , and denote

v ( t ) = u ( t ) , t ( a , b ] . (5.4)

The differential equation, the initial conditions of (5.1), the definition (5.2) of Y and the notation (5.4) imply that

r a t f ( u , v ) = lim r a + r r t f ( u , v ) = lim r a + r r t ( p v ) = ( p v ) ( t ) lim r a + ( p v ) ( r ) ) = ( p v ) ( t ) c ( u , v ) , t ( a , b ] ,

and

u ( t ) - d ( u , v ) = lim r a + ( u ( t ) - u ( r ) ) = lim r a + r t u ( s ) d s (1) = a t u ( s ) d s = a t v ( s ) d s , t ( a , b ] . (2) (3)

Thus, the integral equations of (5.3) hold.

Conversely, let (u, v) be a solution of the system (5.3) in L l o c 1 ( a , b ] . The first equation of (5.3) implies that u is a.e. differentiable and v = u', and that the second initial condition of (5.1) is fulfilled. Since v = u', it follows from the second equation of (5.3) that

( p u ) ( t ) = c ( u , u ) + r a t f ( u , u ) , t ( a , b ] . (5.5)

The equation (5.5) implies that p · u' belongs to R [ a , b ] , and that the differential equation and first initial condition of (5.1) hold. Thus u is a solution of the IVP (5.1) in Y. □

Assume that Lloc(a, b] is ordered a.e. pointwise, that Y is ordered pointwise, and that the functions p, f, c and d satisfy the following hypotheses:

Our main existence and comparison result for the IVP (5.1) reads as follows.

Theorem 5.1. Assume that p : [a, b] +, that 1 p L 1 [ a , b ] , and that the mappings f : L l o c 1 ( a , b ] A R [ a , b ] and c , d : L l o c 1 ( a , b ] are increasing and order-bounded. Then, the IVP (5.1) has the smallest and greatest solutions in Y, and they are increasing with respect to f, c and d.

Proof: The hypotheses imposed on f, c and d imply that the following conditions are valid.

(f0) f(u, v) is RP integrable on [a, b] for every ( u , v ) L l o c 1 ( a , b ] , and there exist such h - , h + A R [ a , b ] that h- f (u1, v1) ≼ f (u2, v2) ≼ h+ for all u i , v i L l o c 1 ( a , b ] , i = 1, 2, u1 u2 and v1 v2.

(c0) c± ∈ ℝ, and c- c(u1, v1) ≤ c(u2, v2) ≤ c+ whenever u i , v i L l o c 1 ( a , b ] , i = 1, 2, u1 u2 and v1 v2.

(d0) d± ∈ ℝ, and d- ≤ d(u1, v1) ≤ d(u2, v2) ≤ d+ whenever u i , v i L l o c 1 ( a , b ] , i = 1, 2, u1 u2 and v1 v2.

Assume that P = L l o c 1 ( a , b ] 2 is ordered componentwise. We shall first show that the vector-functions x+, x- given by

x ( t ) : = d + a t 1 p ( s ) c + r a s h d s , 1 p ( t ) c + r a t h (5.6)

define functions x± P. Since 1/p is Lebesgue integrable and the functions t c + a t h belong to R [ a , b ] , then the second components of x± belong to L l o c 1 ( a , b ] . This result implies that the first components of x± are defined and continuous, whence they belong to L l o c 1 ( a , b ] .

Similarly, by applying also the given hypotheses one can verify that the relations

G 1 ( u , v ) ( t ) : = d ( u , v ) + K a t v ( s ) d s , t ( a , b ] , G 2 ( u , v ) ( t ) : = 1 p ( t ) c ( u , v ) + r a t f ( u , v ) , t J , (5.7)

define an increasing mapping G = (G1, G2) : [x-, x+] → [x-, x+].

Let W be a well-ordered chain in the range of G. The sets W1 = {u : (u, v) ∈ W} and W2 = {v : (u, v) ∈ W} are well-ordered and order-bounded chains in L l o c 1 ( a , b ] . It then follows from [[1], Proposition 9.36] that the supremums of W1 and W2 exist in L l o c 1 ( a , b ] . Obviously, (sup W1, sup W2) is the supremum of W in P. Similarly one can show that each inversely well-ordered chain of the range of G has the infimum in P.

The above proof shows that the operator G = (G1, G2) defined by (5.7) satisfies the hypotheses of Lemma 2.1, and therefore G has the smallest fixed point x* = (u*,v*) and the greatest fixed point x* = (u*, v*). It follows from (5.7) that (u*, v*) and (u*, v*) are solutions of the system (5.3). According to Lemma 5.1, u* and u* belong to Y and are solutions of the IVP (5.1).

To prove that u* and u* are the smallest and greatest of all solutions of (5.1) in Y , let u Y be any solution of (5.1). In view of Lemma 5.1, (u, v) = (u, u') is a solution of the system (5.3). Applying the hypotheses (f0), (c0) and (d0) it is easy to show that x = (u, v) ∈ [x-, x+], where x± are defined by (5.6). Thus x = (u, v) is a fixed point of G = (G1, G2) : [x-, x+] → [x-, x+], defined by (5.7). Because x* = (u*, v*) and x* = (u*, v*) are the smallest and greatest fixed points of G, then (u*, v*) (u, v) (u*, v*). In particular, u* u u*, whence u* and u* are the smallest and greatest of all solutions of the IVP (5.1).

The last assertion is an easy consequence of the last conclusion of Lemma 2.1 and the definition (5.7) of G = (G1, G2). □

Consider next the the following special case of (5.1) where the values of f are combined with impulses and a Henstock-Kurzweil integrable function:

f ( u , v ) ( t ) = λ Λ I ( λ , u , v ) ( λ ) δ ( t - λ ) + q ( t , u , v ) , u , v L l o c 1 ( a , b ] , t [ a , b ] .

In this case problem (5.1) can be rewritten as

d d t ( p ( t ) u ( t ) ) = λ Λ I ( λ , u , u ) δ ( t - λ ) + q ( t , u , u ) a . e . o n [ a , b ] , lim t a + ( p ( t ) u ( t ) ) = c ( u , u ) , lim t a + u ( t ) = d ( u , u ) . (5.8)

The next result is a consequence of Theorem 5.1.

Corollary 5.1. Assume that p : [a, b] → ℝ+, 1 p L 1 [ a , b ] , that functions c , d : L l o c 1 ( a , b ] 2 are increasing and order-bounded, and that the mappings q : [ a , b ] × L l o c 1 ( a , b ] 2 and I : Λ × L l o c 1 ( a , b ] 2 satisfies the following hypotheses.

(q1) q(·, x) is Henstock-Kurzweil integrable on [a, b] for all x L l o c 1 ( a , b ] 2 .

(q2) There exist Henstock-Kurzweil integrable functions q± : [a, b] → ℝ such that K a t q ( s ) d s K a t q ( s , x ) d s K a t q ( s , y ) d s K a t q + ( s ) d s , t [ a , b ] , whenever × ≤ y in L l o c 1 ( a , b ] .

(I) λ Λ | I ( λ , x ) | M < for all x L l o c 1 ( a , b ] 2 , and × α I(λ, x) is increasing for all λ ∈ Λ.

Then, the impulsive IVP (5.8) has the smallest and greatest solutions that are increasing with respect to q, c and d.

Example 5.1. Determine the smallest and greatest solutions of the following singular impulsive IVP.

{ d d t ( t u ( t ) ) = tanh ( [ 20 u [ 1 ] + 10 u [ 1 ] ] 100 ) δ ( t 1 2 ) + [ 1 2 ( u ( s ) + u ( s ) ) d s ] 1 + | [ 1 2 ( u ( s ) + u ( s ) ] | d d t ( t sin 1 t + t ) a .e . on  [ 0 , 3 ] , lim t 0 + t u ( t ) = [ u ( 1 ) ] 1 + | [ u ( 1 ) ] | , lim t 0 + u ( t ) = [ u ( 1 ) ] 1 + | [ u ( 1 ) ] | . (5.9)

Solution: System (5.9) is a special case of (5.8) by setting a = 0, b = 3, p ( t ) = t , Λ = { 1 2 } , and q, c, d and I are given by

q ( t , u , v ) = + [ 1 2 ( u ( s ) + v ( s ) ) d s ] 1 + | [ 1 2 ( u ( s ) + v ( s ) ) d s ] | d d t t sin 1 t + t , c ( u , v ) = [ v ( 1 ) ] 1 + | [ v ( 1 ) ] | , d ( u , v ) = [ u ( 1 ) ] 1 + | [ u ( 1 ) ] | , I ( 1 2 , u , v ) = tanh ( [ 2 0 u [ 1 ] + 1 0 v [ 1 ] ] 1 0 0 ) . (5.10)

It is easy to verify that the hypotheses of Corollary 5.1 hold. Thus (5.9) has the smallest and greatest solutions. The functions x- and x+ defined by (5.6) can be calculated, and their first components are:

u + ( t ) = 1 - 2 2 π 3 + 2 t + 2 t t 3 sin 1 t + 4 t 3 cos 1 t (1) + 4 2 π 3 F r e s n e l S 2 π t + 2 t t 3 + ( 2 t - 2 ) H t - 1 2 , (2) u - = - u + , (3) (4)

where

F r e s n e l S ( x ) = 0 x sin π 2 t 2 d t

is the Fresnel sine integral. According to Lemma 5.1, the smallest solution of (5.9) is equal to the first component of the smallest fixed point of G = (G1, G2), defined by (5.7), with f, c and d given by (5.10) and p ( t ) = t . Calculating the iterations Gnx- it turns out that G4x- = G5x-, whence G 1 4 x - is the smallest solution of (5.9). Similarly, one can show that G 1 3 x + is the greatest solution of (5.9). The exact expressions of these solutions are

u * ( t ) = - 4 5 + 3 2 π 5 - 8 5 t - 3 t t 5 sin 1 t - 6 t 5 cos 1 t (1) - 6 2 π 5 F r e s n e l S 2 π t - tanh 1 1 1 0 ( 2 t - 2 ) H t - 1 2 , (2) u * ( t ) = 3 4 - 1 6 2 π 2 7 + 3 t 2 + 1 6 t t 2 7 sin 1 t + 3 2 t 2 7 cos 1 t (3) + 3 2 2 π 2 7 F r e s n e l S 2 π t + 4 t t 7 + tanh 2 1 2 0 ( 2 t - 2 ) H t - 1 2 . (4) (5)

### 6 Second Order Boundary Value Problems

This section is devoted to the study of the second order boundary value problem (BVP)

- ( p u ) = f ( u , u ) , lim t a + ( p u ) ( t ) = c ( u , u ) , u ( b ) = d ( u , u ) , (6.1)

where f : L 1 [ a , b ] 2 A R [ a , b ] , c, d : L1[a, b]2 → ℝ, and p : [a, b] → ℝ+, -∞ < a < b < ∞. Now we are looking for the smallest and greatest solutions of (6.1) from the set

Z = { u : ( a , b ] : u L 1 [ a , b ] , lim t a + ( p u ) ( t ) e x i s t s , a n d p u R [ a , b ] } . (6.2)

The BVP (6.1) can be transformed into a system of integral equations as follows.

Lemma 6.1. Assume that p : [a, b] → ℝ+, that 1 p L 1 [ a , b ] , and that f ( u , v ) A R [ a , b ] for all u, v L1[a, b]. Then u is a solution of the IVP (6.1) in Z if and only if (u, u') = (u, v), where (u, v) ∈ L1[a, b]2 is a solution of the system

u ( t ) = d ( u , v ) - t b v ( s ) d s , t [ a , b ] , v ( t ) = 1 p ( t ) c ( u , v ) - r a t f ( u , v ) , t [ a , b ] , (6.3)

Proof: Assume that u is a solution of the BVP (6.1) in Z, and denote

v ( t ) = u ( t ) , t [ a , b ] . (6.4)

The differential equation, the boundary conditions of (6.1), the definition (6.2) of Z and the notation (6.4) ensure that

r a t f ( u , v ) = lim r a + r r t f ( u , v ) = lim r a + r r t ( p v ) = lim r a + ( p v ) ( t ) ( p v ) ( r ) ) = ( p v ) ( t ) c ( u , v ) , t [ a , b ] ,

and

u ( t ) - d ( u , v ) = u ( t ) - u ( b ) = - t b u ( s ) d s = - t b v ( s ) d s , t [ a , b ] .

Thus the integral equations of (6.3) hold.

Conversely, let (u, v) be a solution of the system (6.3) in L1[a, b]2. The first equation of (6.3) implies that u is a.e. differentiable and v = u', and that the second boundary condition of (6.1) holds. Since v = u', it follows from the second equation of (6.3) that

( p u ) ( t ) = c ( u , u ) - r a t f ( u , u ) , t [ a , b ] . (6.5)

This equation implies that p · u' belongs to R [ a , b ] , and that the differential equation and first boundary condition of (6.1) are satisfied. Thus u, is a solution of the BVP (6.1) in Z. □

Assume that L1[a, b] is ordered a.e. pointwise, that Z is ordered pointwise. We shall impose the following hypotheses for the functions p, f, c, and d.

(p1) p : [a, b] → ℝ+, and 1 p L 1 [ a , b ] .

(f1) f : L 1 [ a , b ] 2 A R [ a , b ] is order-bounded, and f (u1, v1) ≼ f (u2, v2) whenever ui, vi L1[a, b], i = 1, 2, u1 ≤ u2, and v1 v2.

(c1) c : L1[a, b]2 → ℝ is order-bounded, and c(u2, v2) ≤ c(u1, v1) whenever ui, vi L1[a, b], i = 1, 2, u1 u2, and v1 v2.

(d1) d : L1[a, b]2 → ℝ is order-bounded, and d(u1, v1) ≤ d(u2, v2) whenever ui, vi L1[a, b], i = 1, 2, u1 u2 and v1 v2.

The next theorem is our main existence and comparison result for the BVP (6.1).

Theorem 6.1. Assume that the hypotheses (p1), (f1), (c1), and (d1) hold. Then, the BVP (6.1) has the smallest and greatest solutions in Z, and they are increasing with respect to f and d and decreasing with respect to c.

Proof: Because f, c and d are order-bounded, then the following conditions are valid.

(f0) There exist h A R [ a , b ] such that h- f (u, v) ≼ h+ for all u, v L1[a, b].

(c0) There exist c± ∈ ℝ such that c- c(u, v) ≤ c+ whenever u, v L1[a, b].

(d0) There exist d± ∈ ℝ such that d- d(u, v) ≤ d+ whenever u, v L1[a, b].

Assume that P = L1[a, b]2 is ordered by

( u 1 , v 1 ) ( u 2 , v 2 ) i f a n d o n l y i f u 1 u 2 , a n d v 1 v 2 . (6.6)

We shall first show that the vector-functions x+, x- given by

x - ( t ) = d - - t b 1 p ( s ) ( c + - r a s h - ) d s , 1 p ( t ) ( c + - r a t h - ) , x + ( t ) = d + - t b 1 p ( s ) ( c - - r a s h + ) d s , 1 p ( t ) ( c - - r a t h + ) . (6.7)

belong to P. Since 1/p is Lebesgue integrable and the function t c + - r a t h - belongs to R [ a , b ] , then the second component of x+ is Lebesgue integrable on [a, b]. Similarly one can show that the second component of x- belongs to L1[a, b]. These results ensure that the first components of x± are defined and continuous in t, and hence are in L1[a, b].

Similarly, by applying the given hypotheses one can verify that the relations

G 1 ( u , v ) ( t ) : = d ( u , v ) - K t b v ( s ) d s , t [ a , b ] , G 2 ( u , v ) ( t ) : = 1 p ( t ) c ( u , v r a t f ( u , v ) ) , t [ a , b ] (6.8)

define an increasing mapping G = (G1, G2) : [x-, x+] → [x- , x+].

Let W be a well-ordered chain in the range of G. The set W1 = {u : (u, v) ∈ W} is well ordered, W2 = {v : (u, v) ∈ W } is inversely well-ordered, and both W1 and W2 are order-bounded in L1[a, b]. It then follows from [1, Lemma 9.32] that the supremum of W1 and the infimum of W2 exist in L1[a, b]. Obviously, (sup W1, inf W2) is the supremum of W in (P, ≤). Similarly, one can show that each inversely well-ordered chain of the range of G has the infimum in (P, ≤).

The above proof shows that the operator G = (G1, G2) defined by (6.8) satisfies the hypotheses of Lemma 2.1, whence G has the smallest fixed point x* = (u*, v*) and a greatest fixed point x* = (u*, v*). It follows from (6.8) that (u*, v*) and (u*, v*) are solutions of the system (6.3). According to Lemma 6.1, u* and u* belong to Z and are solutions of the BVP (6.1).

To prove that u* and u* are the smallest and greatest of all solutions of (6.1) in Z, let u Z be any solution of (6.1). In view of Lemma 6.1, (u, v) = (u, u') is a solution of the system (6.3). Applying the properties (f0), (c0), and (d0) it is easy to show that x = (u, v) ∈ [x-, x+], where x± are defined by (6.7). Thus, x = (u, v) is a fixed point of G = (G1, G2) : [x-, x+] → [x- , x+], defined by (6.8). Because x* = (u*, v*) and x* = (u*, v*) are the smallest and greatest fixed points of G, respectively, then (u*, v*) ≤ (u, v) ≤ (u*, v*). In particular, u* u u*, whence u* and u* are the smallest and greatest of all solutions of the BVP (6.1).

The last assertion is an easy consequence of the last conclusion of Lemma 2.1, and the definition (6.8) of G = (G1, G2). □

Consider next a special case of (6.1) where the values of f combined with impulses and Henstock-Kurzweil integrable functions:

d d t ( p ( t ) u ( t ) ) = λ Λ α ( λ ) δ ( t - λ ) + g ( u , u ) ( t ) a . e . o n [ a , b ] , lim t a + ( p ( t ) u ( t ) ) = c ( u , u ) , u ( b ) = d ( u , u ) . (6.9)

Corollary 6.1. Assume that p : [a, b] → ℝ+, 1 p L 1 [ a , b ] , that functions c, d : L1[a, b]2 → ℝ satisfy the hypotheses (ci) and (di), i = 1, 2, that α : Λ → ℝ, λ Λ | α ( λ | < , and that g satisfies the following hypotheses.

(g1) g(u, v) is Henstock-Kurzweil integrable on [a, b] for all u, v L1[a, b].

(g2) There exist Henstock-Kurzweil integrable functions g ¯ , g ¯ : [ a , b ] such that K a t g ¯ K a t g ( u 1 , v 1 ) K a t g ( u 2 , v 2 ) K a t g ¯ , t [ a , b ] , whenever u1 u2 and v1 v2 in L1[a, b].

Then, the impulsive BVP (6.9) has the smallest and greatest solutions that are increasing with respect to g, d and decreasing with respect to c.

Example 6.1. Determine the smallest and greatest solutions of the following singular impulsive BVP.

- d d t ( t u ( t ) ) = δ t - 3 2 + d d t ( t sin 1 t ) tanh 1 1 0 0 0 [ 1 2 ( 3 0 0 0 u ( s ) - 2 0 0 0 u ( s ) ) d s ] a . e . o n [ 0 , 3 ] , lim t 0 + t u ( t ) = [ 1 0 0 0 u ( 1 ) ] 1 + | [ 1 0 0 0 u ( 1 ) ] | , u ( 3 ) = [ 1 0 0 0 u ( 1 ) ] 1 + | [ 1 0 0 0 u ( 1 ) ] | . (6.10)

Solution: System (6.10) is a special case of (6.9) when a = 0, b = 3, p ( t ) = t α 3 2 = 1 , Λ = 3 2 , and g, c, d are given by

g ( u , v ) ( t ) = δ t - 3 2 + d d t ( t sin t t ) tanh 1 1 0 0 0 [ 1 2 ( 3 0 0 0 u ( s ) - 2 0 0 0 v ( s ) ) d s ] , c ( u , v ) = [ 1 0 0 0 v ( 1 ) ] 1 + | [ 1 0 0 0 v ( 1 ) ] | , d ( u , v ) = [ 1 0 0 0 u ( 1 ) ] 1 + | [ 1 0 0 0 u ( 1 ) ] | . (6.11)

It is easy to verify that the hypotheses of Corollary 6.1 are valid. Thus (6.10) has the smallest and greatest solutions. The functions x- and x+ defined by (6.7) can be calculated, and their first components are:

u - ( t ) = - 1 - 6 - 2 3 - 2 3 sin 1 3 - 4 3 3 cos 1 3 - 4 2 π 3 F r e s n e l S 6 3 π + 2 t + 2 t t 3 sin 1 t + 4 t 3 cos 1 t + 4 2 π 3 F r e s n e l S 2 t π + 2 0 t t 3 - ( 2 t - 6 ) H t - 1 3 ,

and

u + ( t ) 1 - 6 + 6 3 + 2 3 sin 1 3 + 4 3 3 cos 1 3 + 4 2 π 3 F r e s n e l S 6 3 π - 2 t - 2 t t 3 sin 1 t - 4 t 3 cos 1 t - 4 2 π 3 F r e s n e l S 2 t π - 2 t t 3 - ( 2 t - 6 ) H t - 1 3 ,

where FresnelS is the Fresnel sine integral.

According to Lemma 6.1 the smallest solution of (6.10) is equal to the first component of the smallest fixed point of G = (G1, G2), defined by (6.3). Calculating the first iterations Gnx- it turns out that G6x- = G7x- . Thus G 1 6 x - is the smallest solution of (6.10). Similarly, one can show that G3x+ = G4x+, whence G 1 3 x + is the greatest solution of (6.10). The exact expressions of these solutions are

u * ( t ) = - 5 6 2 3 5 6 2 4 - 1 1 4 2 1 3 - 2 3 tanh 3 5 5 9 2 0 0 ( 1 + sin 1 3 ) - 6 - 4 3 3 tanh 3 5 5 9 2 0 0 cos 1 3 + 4 2 π 3 tanh 3 5 5 9 2 0 0 F r e s n e l S 6 3 π + 2 8 4 1 t 1 4 2 1 + 2 t t 3 tanh 3 5 5 9 2 0 0 ( sin 1