About this Journal Submit a Manuscript Table of Contents 
Abstract and Applied Analysis
VolumeΒ 2011Β (2011), Article IDΒ 631412, 19 pages
doi:10.1155/2011/631412
Research Article

Boundary-Value Problems for Weakly Nonlinear Delay Differential Systems

A. Boichuk,1,2Β J. Diblík,3,4Β D. Khusainov,5Β and M. Růžičková2

1Institute of Mathematics, National Academy of Sciences of Ukraine, Tereshchenkovskaya Street 3, 01601 Kyiv, Ukraine
2Department of Mathematics, University of Žilina, Univerzitná 8215/1, 01026 Žilina, Slovakia
3Department of Mathematics and Descriptive Geometry, Faculty of Civil Engineering, Brno University of Technology, Veveří 331/95, 602 00 Brno, Czech Republic
4Department of Mathematics, Faculty of Electrical Engineering and Communication, Brno University of Technology, Technická 10, 616 00 Brno, Czech Republic
5Department of Complex System Modeling, Faculty of Cybernetics, Taras, Shevchenko National University of Kyiv, Vladimirskaya Street 64, 01033 Kyiv, Ukraine

Received 30 January 2011; Accepted 31 March 2011

Academic Editor: ElenaΒ Braverman

Copyright Β© 2011 A. Boichuk et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Conditions are derived of the existence of solutions of nonlinear boundary-value problems for systems of 𝑛 ordinary differential equations with constant coefficients and single delay (in the linear part) and with a finite number of measurable delays of argument in nonlinearity: Μ‡ 𝑧 ( 𝑑 ) = 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) + 𝑔 ( 𝑑 ) + πœ€ 𝑍 ( 𝑧 ( β„Ž 𝑖 ( 𝑑 ) , 𝑑 , πœ€ ) , 𝑑 ∈ [ π‘Ž , 𝑏 ] , assuming that these solutions satisfy the initial and boundary conditions 𝑧 ( 𝑠 ) ∢ = πœ“ ( 𝑠 ) i f 𝑠 βˆ‰ [ π‘Ž , 𝑏 ] , β„“ 𝑧 ( β‹… ) = 𝛼 ∈ ℝ π‘š . The use of a delayed matrix exponential and a method of pseudoinverse by Moore-Penrose matrices led to an explicit and analytical form of sufficient conditions for the existence of solutions in a given space and, moreover, to the construction of an iterative process for finding the solutions of such problems in a general case when the number of boundary conditions (defined by a linear vector functional β„“ ) does not coincide with the number of unknowns in the differential system with a single delay.

1. Introduction

First, we derive some auxiliary results concerning the theory of differential equations with delay. Consider a system of linear differential equations with concentrated delay Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 ( 𝑑 ) 𝑧 ξ€· β„Ž 0 ( 𝑑 ) ξ€Έ = 𝑔 ( 𝑑 ) i f 𝑑 ∈ [ π‘Ž , 𝑏 ] , ( 1 . 1 ) assuming that 𝑧 ( 𝑠 ) ∢ = πœ“ ( 𝑠 ) i f 𝑠 βˆ‰ [ π‘Ž , 𝑏 ] , ( 1 . 2 ) where 𝐴 is an 𝑛 Γ— 𝑛 real matrix and 𝑔 is an 𝑛 -dimensional real column-vector with components in the space 𝐿 𝑝 [ π‘Ž , 𝑏 ] (where 𝑝 ∈ [ 1 , ∞ ) ) of functions summable on [ π‘Ž , 𝑏 ] ; the delay β„Ž 0 ( 𝑑 ) ≀ 𝑑 is a function β„Ž 0 ∢ [ π‘Ž , 𝑏 ] β†’ ℝ measurable on [ π‘Ž , 𝑏 ] ; πœ“ ∢ ℝ ⧡ [ π‘Ž , 𝑏 ] β†’ ℝ 𝑛 is a given function. Using the denotations ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) ∢ = ⎧ βŽͺ ⎨ βŽͺ ⎩ 𝑧 ξ€· β„Ž 0 ( 𝑑 ) ξ€Έ i f β„Ž 0 ( 𝑑 ) ∈ [ π‘Ž , 𝑏 ] , πœƒ i f β„Ž 0 ( 𝑑 ) βˆ‰ [ π‘Ž , 𝑏 ] , ( 1 . 3 ) πœ“ β„Ž 0 ( 𝑑 ) ∢ = ⎧ βŽͺ ⎨ βŽͺ ⎩ πœƒ i f β„Ž 0 ( 𝑑 ) ∈ [ π‘Ž , 𝑏 ] , πœ“ ξ€· β„Ž 0 ( 𝑑 ) ξ€Έ i f β„Ž 0 ( 𝑑 ) βˆ‰ [ π‘Ž , 𝑏 ] , ( 1 . 4 ) where πœƒ is an 𝑛 -dimensional zero column-vector and assuming 𝑑 ∈ [ π‘Ž , 𝑏 ] , it is possible to rewrite (1.1), (1.2) as ( 𝐿 𝑧 ) ( 𝑑 ) ∢ = Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 ( 𝑑 ) ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) = πœ‘ ( 𝑑 ) , 𝑑 ∈ [ π‘Ž , 𝑏 ] , ( 1 . 5 ) where πœ‘ is an 𝑛 -dimensional column-vector defined by the formula πœ‘ ( 𝑑 ) ∢ = 𝑔 ( 𝑑 ) + 𝐴 ( 𝑑 ) πœ“ β„Ž 0 ( 𝑑 ) ∈ 𝐿 𝑝 [ π‘Ž , 𝑏 ] . ( 1 . 6 ) We will investigate (1.5) assuming that the operator 𝐿 maps a Banach space 𝐷 𝑝 [ π‘Ž , 𝑏 ] of absolutely continuous functions 𝑧 ∢ [ π‘Ž , 𝑏 ] β†’ ℝ 𝑛 into a Banach space 𝐿 𝑝 [ π‘Ž , 𝑏 ] ( 1 ≀ 𝑝 < ∞ ) of functions πœ‘ ∢ [ π‘Ž , 𝑏 ] β†’ ℝ 𝑛 summable on [ π‘Ž , 𝑏 ] ; the operator 𝑆 β„Ž 0 maps the space 𝐷 𝑝 [ π‘Ž , 𝑏 ] into the space 𝐿 𝑝 [ π‘Ž , 𝑏 ] . Transformations (1.3), (1.4) make it possible to add the initial function πœ“ ( 𝑠 ) , 𝑠 < π‘Ž to nonhomogeneity generating an additive and homogeneous operation not depending on πœ“ and without the classical assumption regarding the continuous connection of solution 𝑧 ( 𝑑 ) with the initial function πœ“ ( 𝑑 ) at the point 𝑑 = π‘Ž .

A solution of differential system (1.5) is defined as an 𝑛 -dimensional column vector-function 𝑧 ∈ 𝐷 𝑝 [ π‘Ž , 𝑏 ] , absolutely continuous on [ π‘Ž , 𝑏 ] , with a derivative Μ‡ 𝑧 ∈ 𝐿 𝑝 [ π‘Ž , 𝑏 ] satisfying (1.5) almost everywhere on [ π‘Ž , 𝑏 ] .

Such approach makes it possible to apply well-developed methods of linear functional analysis to (1.5) with a linear and bounded operator 𝐿 . It is well-known (see: [1, 2]) that a nonhomogeneous operator equation (1.5) with delayed argument is solvable in the space 𝐷 𝑝 [ π‘Ž , 𝑏 ] for an arbitrary right-hand side πœ‘ ∈ 𝐿 𝑝 [ π‘Ž , 𝑏 ] and has an 𝑛 -dimensional family of solutions ( d i m k e r 𝐿 = 𝑛 ) in the form 𝑧 ( 𝑑 ) = 𝑋 ( 𝑑 ) 𝑐 + ξ€œ 𝑏 π‘Ž 𝐾 ( 𝑑 , 𝑠 ) πœ‘ ( 𝑠 ) 𝑑 𝑠 βˆ€ 𝑐 ∈ ℝ 𝑛 , ( 1 . 7 ) where the kernel 𝐾 ( 𝑑 , 𝑠 ) is an 𝑛 Γ— 𝑛 Cauchy matrix defined in the square [ π‘Ž , 𝑏 ] Γ— [ π‘Ž , 𝑏 ] being, for every fixed 𝑠 ≀ 𝑑 , a solution of the matrix Cauchy problem ( 𝐿 𝐾 ( β‹… , 𝑠 ) ) ( 𝑑 ) ∢ = πœ• 𝐾 ( 𝑑 , 𝑠 ) πœ• 𝑑 βˆ’ 𝐴 ( 𝑑 ) ξ€· 𝑆 β„Ž 0 𝐾 ( β‹… , 𝑠 ) ξ€Έ ( 𝑑 ) = Θ , 𝐾 ( 𝑠 , 𝑠 ) = 𝐼 , ( 1 . 8 ) where 𝐾 ( 𝑑 , 𝑠 ) ≑ Θ if π‘Ž ≀ 𝑑 < 𝑠 ≀ 𝑏 , Θ is 𝑛 Γ— 𝑛 null matrix and 𝐼 is 𝑛 Γ— 𝑛 identity matrix. A fundamental 𝑛 Γ— 𝑛 matrix 𝑋 ( 𝑑 ) for the homogeneous ( πœ‘ ≑ πœƒ ) equation (1.5) has the form 𝑋 ( 𝑑 ) = 𝐾 ( 𝑑 , π‘Ž ) , 𝑋 ( π‘Ž ) = 𝐼 [2]. Throughout the paper, we denote by Θ 𝑠 an 𝑠 Γ— 𝑠 null matrix if 𝑠 β‰  𝑛 , by Θ 𝑠 , 𝑝 an 𝑠 Γ— 𝑝 null matrix, by 𝐼 𝑠 an 𝑠 Γ— 𝑠 identity matrix if 𝑠 β‰  𝑛 , and by πœƒ 𝑠 an 𝑠 -dimensional zero column-vector if 𝑠 β‰  𝑛 .

A serious disadvantage of this approach, when investigating the above-formulated problem, is the necessity to find the Cauchy matrix 𝐾 ( 𝑑 , 𝑠 ) [3, 4]. It exists but, as a rule, can only be found numerically. Therefore, it is important to find systems of differential equations with delay such that this problem can be solved directly. Below we consider the case of a system with so-called single delay [5]. In this case, the problem of how to construct the Cauchy matrix is successfully solved analytically due to a delayed matrix exponential defined below.

1.1. A Delayed Matrix Exponential

Consider a Cauchy problem for a linear nonhomogeneous differential system with constant coefficients and with a single delay 𝜏 Μ‡ 𝑧 ( 𝑑 ) = 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) + 𝑔 ( 𝑑 ) , ( 1 . 9 ) 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) , i f 𝑠 ∈ [ βˆ’ 𝜏 , 0 ] , ( 1 . 1 0 ) with an 𝑛 Γ— 𝑛 constant matrix 𝐴 , 𝑔 ∢ [ 0 , ∞ ) β†’ ℝ 𝑛 , πœ“ ∢ [ βˆ’ 𝜏 , 0 ] β†’ ℝ 𝑛 , 𝜏 > 0 and an unknown vector-solution 𝑧 ∢ [ βˆ’ 𝜏 , ∞ ) β†’ ℝ 𝑛 . Together with a nonhomogeneous problem (1.9), (1.10), we consider a related homogeneous problem Μ‡ 𝑧 ( 𝑑 ) = 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) , 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) , i f 𝑠 ∈ [ βˆ’ 𝜏 , 0 ] . ( 1 . 1 1 )

Denote by 𝑒 𝐴 𝑑 𝜏 a matrix function called a delayed matrix exponential (see [5]) and defined as 𝑒 𝐴 𝑑 𝜏 ∢ = ⎧ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ ⎨ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ βŽͺ ⎩ Θ i f βˆ’ ∞ < 𝑑 < βˆ’ 𝜏 , 𝐼 i f βˆ’ 𝜏 ≀ 𝑑 < 0 , 𝐼 + 𝐴 𝑑 1 ! i f 0 ≀ 𝑑 < 𝜏 , 𝐼 + 𝐴 𝑑 1 ! + 𝐴 2 ( 𝑑 βˆ’ 𝜏 ) 2 2 ! i f 𝜏 ≀ 𝑑 < 2 𝜏 , β‹― 𝐼 + 𝐴 𝑑 1 ! + β‹― + 𝐴 π‘˜ ( 𝑑 βˆ’ ( π‘˜ βˆ’ 1 ) 𝜏 ) π‘˜ π‘˜ ! i f ( π‘˜ βˆ’ 1 ) 𝜏 ≀ 𝑑 < π‘˜ 𝜏 , β‹― . ( 1 . 1 2 ) This definition can be reduced to the following expression: 𝑒 𝐴 𝑑 𝜏 = [ 𝑑 / 𝜏 ] + 1  𝑛 = 0 𝐴 𝑛 ( 𝑑 βˆ’ ( 𝑛 βˆ’ 1 ) 𝜏 ) 𝑛 𝑛 ! , ( 1 . 1 3 ) where [ 𝑑 / 𝜏 ] is the greatest integer function. The delayed matrix exponential equals the unit matrix 𝐼 on [ βˆ’ 𝜏 , 0 ] and represents a fundamental matrix of a homogeneous system with single delay. Thus, the delayed matrix exponential solves the Cauchy problem for a homogeneous system (1.11), satisfying the unit initial conditions 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) ≑ 𝑒 𝐴 𝑠 𝜏 = 𝐼 i f βˆ’ 𝜏 ≀ 𝑠 ≀ 0 , ( 1 . 1 4 ) and the following statement holds (see, e.g., [5], [6, Remark 1], [7, Theorem 2.1]).

Lemma 1.1. A solution of a Cauchy problem for a nonhomogeneous system with single delay (1.9), satisfying a constant initial condition 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) = 𝑐 ∈ ℝ 𝑛 i f 𝑠 ∈ [ βˆ’ 𝜏 , 0 ] ( 1 . 1 5 ) has the form 𝑧 ( 𝑑 ) = 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 ) 𝜏 𝑐 + ξ€œ 𝑑 0 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 βˆ’ 𝑠 ) 𝜏 𝑔 ( 𝑠 ) 𝑑 𝑠 . ( 1 . 1 6 )

The delayed matrix exponential was applied, for example, in [6, 7] to investigation of boundary value problems of diffferential systems with a single delay and in [8] to investigation of the stability of linear perturbed systems with a single delay.

1.2. Fredholm Boundary-Value Problem

Without loss of generality, let π‘Ž = 0 and, with a view of the above, the problem (1.9), (1.10) can be transformed ( β„Ž 0 ( 𝑑 ) ∢ = 𝑑 βˆ’ 𝜏 ) to an equation of the type (1.1) (see (1.5)) Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) = πœ‘ ( 𝑑 ) , 𝑑 ∈ [ 0 , 𝑏 ] , ( 1 . 1 7 ) where, in accordance with (1.3),(1.4) ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) = ⎧ βŽͺ ⎨ βŽͺ ⎩ 𝑧 ( 𝑑 βˆ’ 𝜏 ) i f 𝑑 βˆ’ 𝜏 ∈ [ 0 , 𝑏 ] , πœƒ i f 𝑑 βˆ’ 𝜏 βˆ‰ [ 0 , 𝑏 ] , πœ‘ ( 𝑑 ) = 𝑔 ( 𝑑 ) + 𝐴 πœ“ β„Ž 0 ( 𝑑 ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , πœ“ β„Ž 0 ( 𝑑 ) = ⎧ βŽͺ ⎨ βŽͺ ⎩ πœƒ i f 𝑑 βˆ’ 𝜏 ∈ [ 0 , 𝑏 ] , πœ“ ( 𝑑 βˆ’ 𝜏 ) i f 𝑑 βˆ’ 𝜏 βˆ‰ [ 0 , 𝑏 ] . ( 1 . 1 8 )

A general solution of problem (1.17) for a nonhomogeneous system with single delay and zero initial data has the form (1.7) 𝑧 ( 𝑑 ) = 𝑋 ( 𝑑 ) 𝑐 + ξ€œ 𝑏 0 𝐾 ( 𝑑 , 𝑠 ) πœ‘ ( 𝑠 ) 𝑑 𝑠 βˆ€ 𝑐 ∈ ℝ 𝑛 , ( 1 . 1 9 ) where, as can easily be verified (in view of the above-defined delayed matrix exponential) by substituting into (1.17), 𝑋 ( 𝑑 ) = 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 ) 𝜏 , 𝑋 ( 0 ) = 𝑒 βˆ’ 𝐴 𝜏 𝜏 = 𝐼 ( 1 . 2 0 ) is a normal fundamental matrix of the homogeneous system related to (1.17) (or (1.9)) with initial data 𝑋 ( 0 ) = 𝐼 , and the Cauchy matrix 𝐾 ( 𝑑 , 𝑠 ) has the form 𝐾 ( 𝑑 , 𝑠 ) = 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 βˆ’ 𝑠 ) 𝜏 i f 0 ≀ 𝑠 < 𝑑 ≀ 𝑏 , 𝐾 ( 𝑑 , 𝑠 ) ≑ Θ i f 0 ≀ 𝑑 < 𝑠 ≀ 𝑏 . ( 1 . 2 1 ) Obviously 𝐾 ( 𝑑 , 0 ) = 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 ) 𝜏 = 𝑋 ( 𝑑 ) , 𝐾 ( 0 , 0 ) = 𝑒 𝐴 ( βˆ’ 𝜏 ) 𝜏 = 𝑋 ( 0 ) = 𝐼 , ( 1 . 2 2 ) and, therefore, the initial problem (1.17) for systems of ordinary differential equations with constant coefficients and single delay has an 𝑛 -parametric family of linearly independent solutions (1.16).

Now, we will deal with a general boundary-value problem for system (1.17). Using the results [2, 9], it is easy to derive statements for a general boundary-value problem if the number π‘š of boundary conditions does not coincide with the number 𝑛 of unknowns in a differential system with single delay.

We consider a boundary-value problem Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) = 𝑔 ( 𝑑 ) , 𝑑 ∈ [ 0 , 𝑏 ] , 𝑧 ( 𝑠 ) ∢ = πœ“ ( 𝑠 ) , 𝑠 βˆ‰ [ 0 , 𝑏 ] , ( 1 . 2 3 ) assuming that β„“ 𝑧 ( β‹… ) = 𝛼 ∈ ℝ π‘š , ( 1 . 2 4 ) or, using (1.18), its equivalent form Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) = πœ‘ ( 𝑑 ) , 𝑑 ∈ [ 0 , 𝑏 ] , β„“ 𝑧 ( β‹… ) = 𝛼 ∈ ℝ π‘š , ( 1 . 2 5 ) where 𝛼 is an π‘š -dimensional constant vector-column β„“ is an π‘š -dimensional linear vector-functional defined on the space 𝐷 𝑝 [ 0 , 𝑏 ] of an 𝑛 -dimensional vector-functions β„“ = c o l ξ€· β„“ 1 , … , β„“ π‘š ξ€Έ ∢ 𝐷 𝑝 [ 0 , 𝑏 ] ⟢ ℝ π‘š , β„“ 𝑖 ∢ 𝐷 𝑝 [ 0 , 𝑏 ] ⟢ ℝ , 𝑖 = 1 , … , π‘š , ( 1 . 2 6 ) absolutely continuous on [ 0 , 𝑏 ] . Such problems for functional-differential equations are of Fredholm's type (see, e.g., [1, 2]). In order to formulate the following result, we need several auxiliary abbreviations. We set 𝑄 ∢ = β„“ 𝑋 ( β‹… ) = β„“ 𝑒 𝐴 ( β‹… βˆ’ 𝜏 ) 𝜏 . ( 1 . 2 7 ) We define an 𝑛 Γ— 𝑛 -dimensional matrix (orthogonal projection) 𝑃 𝑄 ∢ = 𝐼 βˆ’ 𝑄 + 𝑄 , ( 1 . 2 8 ) projecting space ℝ 𝑛 to k e r 𝑄 of the matrix 𝑄 .

Moreover, we define an π‘š Γ— π‘š -dimensional matrix (orthogonal projection) 𝑃 𝑄 βˆ— ∢ = 𝐼 π‘š βˆ’ 𝑄 𝑄 + , ( 1 . 2 9 ) projecting space ℝ π‘š to k e r 𝑄 βˆ— of the transposed matrix 𝑄 βˆ— = 𝑄 𝑇 , where 𝐼 π‘š is an π‘š Γ— π‘š identity matrix and 𝑄 + is an 𝑛 Γ— π‘š -dimensional matrix pseudoinverse to the π‘š Γ— 𝑛 -dimensional matrix 𝑄 . Denote 𝑑 ∢ = r a n k 𝑃 𝑄 βˆ— and 𝑛 1 ∢ = r a n k 𝑄 = r a n k 𝑄 βˆ— . Since r a n k 𝑃 𝑄 βˆ— = π‘š βˆ’ r a n k 𝑄 βˆ— , ( 1 . 3 0 ) we have 𝑑 = π‘š βˆ’ 𝑛 1 .

We will denote by 𝑃 𝑄 βˆ— 𝑑 an 𝑑 Γ— π‘š -dimensional matrix constructed from 𝑑 linearly independent rows of the matrix 𝑃 𝑄 βˆ— . Denote π‘Ÿ ∢ = r a n k 𝑃 𝑄 . Since r a n k 𝑃 𝑄 = 𝑛 βˆ’ r a n k 𝑄 , ( 1 . 3 1 ) we have π‘Ÿ = 𝑛 βˆ’ 𝑛 1 . By 𝑃 𝑄 π‘Ÿ we will denote an 𝑛 Γ— π‘Ÿ -dimensional matrix constructed from π‘Ÿ linearly independent columns of the matrix 𝑃 𝑄 . Finally, we define 𝑋 π‘Ÿ ( 𝑑 ) ∢ = 𝑋 ( 𝑑 ) 𝑃 𝑄 π‘Ÿ , ( 1 . 3 2 ) and a generalized Green operator ( 𝐺 πœ‘ ) ( 𝑑 ) ∢ = ξ€œ 𝑏 0 𝐺 ( 𝑑 , 𝑠 ) πœ‘ ( 𝑠 ) 𝑑 𝑠 , ( 1 . 3 3 ) where 𝐺 ( 𝑑 , 𝑠 ) ∢ = 𝐾 ( 𝑑 , 𝑠 ) βˆ’ 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 ) 𝜏 𝑄 + β„“ 𝐾 ( β‹… , 𝑠 ) ( 1 . 3 4 ) is a generalized Green matrix corresponding to the boundary-value problem (1.25) (the Cauchy matrix 𝐾 ( 𝑑 , 𝑠 ) has the form (1.21)).

In [6, Theorem 4], the following result (formulating the necessary and sufficient conditions of solvability and giving representations of the solutions 𝑧 ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑧 ∈ 𝐿 𝑝 [ 0 , 𝑏 ] of the boundary-value problem (1.25) in an explicit analytical form) is proved.

Theorem 1.2. If 𝑛 1 ≀ m i n ( π‘š , 𝑛 ) , then: (i)the homogeneous problem Μ‡ 𝑧 ( 𝑑 ) βˆ’ 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) = πœƒ , 𝑑 ∈ [ 0 , 𝑏 ] , β„“ 𝑧 ( β‹… ) = πœƒ π‘š ∈ ℝ π‘š ( 1 . 3 5 ) corresponding to problem (1.25) has exactly π‘Ÿ linearly independent solutions 𝑧 ξ€· 𝑑 , 𝑐 π‘Ÿ ξ€Έ = 𝑋 π‘Ÿ ( 𝑑 ) 𝑐 π‘Ÿ = 𝑒 𝐴 ( 𝑑 βˆ’ 𝜏 ) 𝜏 𝑃 𝑄 π‘Ÿ 𝑐 π‘Ÿ ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , ( 1 . 3 6 ) (ii)nonhomogeneous problem (1.25) is solvable in the space 𝐷 𝑝 [ 0 , 𝑏 ] if and only if πœ‘ ∈ 𝐿 𝑝 [ 0 , 𝑏 ] and 𝛼 ∈ ℝ π‘š satisfy 𝑑 linearly independent conditions 𝑃 𝑄 βˆ— 𝑑 β‹… ξ‚΅ 𝛼 βˆ’ β„“ ξ€œ 𝑏 0 𝐾 ( β‹… , 𝑠 ) πœ‘ ( 𝑠 ) 𝑑 𝑠 ξ‚Ά = πœƒ 𝑑 , ( 1 . 3 7 ) (iii)in that case the nonhomogeneous problem (1.25) has an π‘Ÿ -dimensional family of linearly independent solutions represented in an analytical form 𝑧 ( 𝑑 ) = 𝑧 0 ξ€· 𝑑 , 𝑐 π‘Ÿ ξ€Έ ∢ = 𝑋 π‘Ÿ ( 𝑑 ) 𝑐 π‘Ÿ + ( 𝐺 πœ‘ ) ( 𝑑 ) + 𝑋 ( 𝑑 ) 𝑄 + 𝛼 βˆ€ 𝑐 π‘Ÿ ∈ ℝ π‘Ÿ . ( 1 . 3 8 )

2. Perturbed Weakly Nonlinear Boundary Value Problems

As an example of applying Theorem 1.2, we consider a problem of the branching of solutions 𝑧 ∢ [ 0 , 𝑏 ] β†’ ℝ 𝑛 , 𝑏 > 0 of systems of nonlinear ordinary differential equations with a small parameter πœ€ and with a finite number of measurable delays β„Ž 𝑖 ( 𝑑 ) , 𝑖 = 1 , 2 , … , π‘˜ of argument of the form Μ‡ 𝑧 ( 𝑑 ) = 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) + 𝑔 ( 𝑑 ) + πœ€ 𝑍 ξ€· 𝑧 ξ€· β„Ž 𝑖 ( 𝑑 ) ξ€Έ , 𝑑 , πœ€ ξ€Έ , 𝑑 ∈ [ 0 , 𝑏 ] , β„Ž 𝑖 ( 𝑑 ) ≀ 𝑑 , ( 2 . 1 ) satisfying the initial and boundary conditions 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) , i f 𝑠 < 0 , β„“ 𝑧 ( β‹… ) = 𝛼 , 𝛼 ∈ ℝ π‘š , ( 2 . 2 ) and such that its solution 𝑧 = 𝑧 ( 𝑑 , πœ€ ) , satisfying 𝑧 ( β‹… , πœ€ ) ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑧 ( β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑧 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , ( 2 . 3 ) for a sufficiently small πœ€ 0 > 0 , for πœ€ = 0 , turns into one of the generating solutions (1.38); that is, 𝑧 ( 𝑑 , 0 ) = 𝑧 0 ( 𝑑 , 𝑐 π‘Ÿ ) for a 𝑐 π‘Ÿ ∈ ℝ π‘Ÿ . We assume that the 𝑛 Γ— 1 vector-operator 𝑍 satisfies 𝑍 ( β‹… , 𝑑 , πœ€ ) ∈ 𝐢 1 ξ€Ί β€– β€– 𝑧 βˆ’ 𝑧 0 β€– β€– ≀ π‘ž ξ€» , 𝑍 ξ€· 𝑧 ξ€· β„Ž 𝑖 ( 𝑑 ) ξ€Έ , β‹… , πœ€ ξ€Έ ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑍 ξ€· 𝑧 ξ€· β„Ž 𝑖 ( 𝑑 ) ξ€Έ , 𝑑 , β‹… ξ€Έ ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , ( 2 . 4 ) where π‘ž > 0 is sufficiently small. Using denotations (1.3), (1.4), and (1.6), it is easy to show that the perturbed nonlinear boundary value problem (2.1), (2.2) can be rewritten in the form Μ‡ 𝑧 ( 𝑑 ) = 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) + πœ€ 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 ξ€Έ ( 𝑑 ) , 𝑑 , πœ€ ξ€Έ + πœ‘ ( 𝑑 ) , β„“ 𝑧 ( β‹… ) = 𝛼 , 𝑑 ∈ [ 0 , 𝑏 ] . ( 2 . 5 ) In (2.5), 𝐴 is an 𝑛 Γ— 𝑛 constant matrix, β„Ž 0 ∢ [ 0 , 𝑏 ] β†’ ℝ is a single delay defined by β„Ž 0 ( 𝑑 ) ∢ = 𝑑 βˆ’ 𝜏 , 𝜏 > 0 , ξ€· 𝑆 β„Ž 𝑧 ξ€Έ ( 𝑑 ) = c o l ξ€Ί ξ€· 𝑆 β„Ž 1 𝑧 ξ€Έ ( 𝑑 ) , … , ξ€· 𝑆 β„Ž π‘˜ 𝑧 ξ€Έ ( 𝑑 ) ξ€» ( 2 . 6 ) is an 𝑁 -dimensional column vector, where 𝑁 = 𝑛 π‘˜ , and πœ‘ is an 𝑛 -dimensional column vector given by πœ‘ ( 𝑑 ) = 𝑔 ( 𝑑 ) + 𝐴 πœ“ β„Ž 0 ( 𝑑 ) . ( 2 . 7 ) The operator 𝑆 β„Ž maps the space 𝐷 𝑝 into the space 𝐿 𝑁 𝑝 = 𝐿 𝑝 Γ— β‹― Γ— 𝐿 𝑝 ξ„Ώ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…ƒ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…€ ξ…Œ π‘˜ - t i m e s , ( 2 . 8 ) that is, 𝑆 β„Ž ∢ 𝐷 𝑝 β†’ 𝐿 𝑁 𝑝 . Using denotation (1.3) for the operator 𝑆 β„Ž 𝑖 ∢ 𝐷 𝑝 β†’ 𝐿 𝑝 , 𝑖 = 1 , … , π‘˜ , we have the following representation: ξ€· 𝑆 β„Ž 𝑖 𝑧 ξ€Έ ( 𝑑 ) = ξ€œ 𝑏 0 πœ’ β„Ž 𝑖 ( 𝑑 , 𝑠 ) Μ‡ 𝑧 ( 𝑠 ) 𝑑 𝑠 + πœ’ β„Ž 𝑖 ( 𝑑 , 0 ) 𝑧 ( 0 ) , ( 2 . 9 ) where πœ’ β„Ž 𝑖 ( 𝑑 , 𝑠 ) = ⎧ βŽͺ ⎨ βŽͺ ⎩ 1 , i f ( 𝑑 , 𝑠 ) ∈ Ξ© 𝑖 , 0 , i f ( 𝑑 , 𝑠 ) βˆ‰ Ξ© 𝑖 ( 2 . 1 0 ) is the characteristic function of the set Ξ© 𝑖 ∢ = ξ€½ ( 𝑑 , 𝑠 ) ∈ [ 0 , 𝑏 ] Γ— [ 0 , 𝑏 ] ∢ 0 ≀ 𝑠 ≀ β„Ž 𝑖 ( 𝑑 ) ≀ 𝑏 ξ€Ύ . ( 2 . 1 1 ) Assume that the generating boundary value problem Μ‡ 𝑧 ( 𝑑 ) = 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) + πœ‘ ( 𝑑 ) , 𝑙 𝑧 = 𝛼 , ( 2 . 1 2 ) being a particular case of (2.5) for πœ€ = 0 , has solutions for nonhomogeneities πœ‘ ∈ 𝐿 𝑝 [ 0 , 𝑏 ] and 𝛼 ∈ ℝ π‘š that satisfy conditions (1.37). In such a case, by Theorem 1.2, the problem (2.12) possesses an π‘Ÿ -dimensional family of solutions of the form (1.38).

Problem 1. Below, we consider the following problem: derive the necessary and sufficient conditions indicating when solutions of (2.5) turn into solutions (1.38) of the boundary value problem (2.12) for πœ€ = 0 .

Using the theory of generalized inverse operators [2], it is possible to find conditions for the solutions of the boundary value problem (2.5) to be branching from the solutions of (2.5) with πœ€ = 0 . Below, we formulate statements, solving the above problem. As compared with an earlier result [10, page 150], the present result is derived in an explicit analytical form. The progress was possible by using the delayed matrix exponential since, in such a case, all the necessary calculations can be performed to the full.

Theorem 2.1 (necessary condition). Consider the system (2.1); that is, Μ‡ 𝑧 ( 𝑑 ) = 𝐴 𝑧 ( 𝑑 βˆ’ 𝜏 ) + 𝑔 ( 𝑑 ) + πœ€ 𝑍 ξ€· 𝑧 ξ€· β„Ž 𝑖 ( 𝑑 ) ξ€Έ , 𝑑 , πœ€ ξ€Έ , 𝑑 ∈ [ 0 , 𝑏 ] , ( 2 . 1 3 ) where β„Ž 𝑖 ( 𝑑 ) ≀ 𝑑 , 𝑖 = 1 , … , π‘˜ , with the initial and boundary conditions (2.2); that is, 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) , i f 𝑠 < 0 < 𝑏 , β„“ 𝑧 ( β‹… ) = 𝛼 ∈ ℝ π‘š , ( 2 . 1 4 ) and assume that, for nonhomogeneities πœ‘ ( 𝑑 ) = 𝑔 ( 𝑑 ) + 𝐴 πœ“ β„Ž 0 ( 𝑑 ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , ( 2 . 1 5 ) and for 𝛼 ∈ ℝ π‘š , the generating boundary value problem Μ‡ 𝑧 ( 𝑑 ) = 𝐴 ξ€· 𝑆 β„Ž 0 𝑧 ξ€Έ ( 𝑑 ) + πœ‘ ( 𝑑 ) , β„“ 𝑧 ( β‹… ) = 𝛼 , ( 2 . 1 6 ) corresponding to the problem (1.25), has exactly an π‘Ÿ -dimensional family of linearly independent solutions of the form (1.38). Moreover, assume that the boundary value problem (2.13), (2.14) has a solution 𝑧 ( 𝑑 , πœ€ ) which, for πœ€ = 0 , turns into one of solutions 𝑧 0 ( 𝑑 , 𝑐 π‘Ÿ ) in (1.38) with a vector-constant 𝑐 π‘Ÿ ∢ = 𝑐 0 π‘Ÿ ∈ ℝ π‘Ÿ .
Then, the vector 𝑐 0 π‘Ÿ satisfies the equation 𝐹 ξ€· 𝑐 0 π‘Ÿ ξ€Έ ∢ = ξ€œ 𝑏 0 𝐻 ( 𝑠 ) 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 0 ξ€Έ ξ€· 𝑠 , 𝑐 0 π‘Ÿ ξ€Έ , 𝑠 , 0 ξ€Έ 𝑑 𝑠 = πœƒ 𝑑 , ( 2 . 1 7 ) where 𝐻 ( 𝑠 ) ∢ = 𝑃 𝑄 βˆ— 𝑑 β„“ 𝐾 ( β‹… , 𝑠 ) = 𝑃 𝑄 βˆ— 𝑑 β„“ 𝑒 𝐴 ( β‹… βˆ’ 𝜏 βˆ’ 𝑠 ) 𝜏 . ( 2 . 1 8 )

Proof. We consider the nonlinearity in system (2.13), that is, the term πœ€ 𝑍 ( 𝑧 ( β„Ž 𝑖 ( 𝑑 ) ) , 𝑑 , πœ€ ) as an inhomogeneity, and use Theorem 1.2 assuming that condition (1.37) is satisfied. This gives ξ€œ 𝑏 0 𝐻 ( 𝑠 ) 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ 𝑑 𝑠 = πœƒ 𝑑 . ( 2 . 1 9 ) In this integral, letting πœ€ β†’ 0 , we arrive at the required condition (2.17).

Corollary 2.2. For periodic boundary-value problems, the vector-constant 𝑐 π‘Ÿ ∈ ℝ π‘Ÿ has a physical meaning-it is the amplitude of the oscillations generated. For this reason, (2.17) is called an equation generating the amplitude [11]. By analogy with the investigation of periodic problems, it is natural to say (2.17) is an equation for generating the constants of the boundary value problem (2.13), (2.14). If (2.17) is solvable, then the vector constant 𝑐 0 π‘Ÿ ∈ ℝ π‘Ÿ specifies the generating solution 𝑧 0 ( 𝑑 , 𝑐 0 π‘Ÿ ) corresponding to the solution 𝑧 = 𝑧 ( 𝑑 , πœ€ ) of the original problem such that 𝑧 ( β‹… , πœ€ ) ∢ [ 0 , 𝑏 ] ⟢ ℝ 𝑛 , 𝑧 ( β‹… , πœ€ ) ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑧 ( β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑧 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , 𝑧 ( 𝑑 , 0 ) = 𝑧 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ . ( 2 . 2 0 ) Also, if (2.17) is unsolvable, the problem (2.13), (2.14) has no solution in the analyzed space. Note that, here and in what follows, all expressions are obtained in the real form and hence, we are interested in real solutions of (2.17), which can be algebraic or transcendental.

Sufficient conditions for the existence of solutions of the boundary-value problem (2.13), (2.14) can be derived using results in [10, page 155] and [2]. By changing the variables in system (2.13), (2.14) 𝑧 ( 𝑑 , πœ€ ) = 𝑧 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ + 𝑦 ( 𝑑 , πœ€ ) , ( 2 . 2 1 ) we arrive at a problem of finding sufficient conditions for the existence of solutions of the problem Μ‡ 𝑦 ( 𝑑 ) = 𝐴 ξ€· 𝑆 β„Ž 0 𝑦 ξ€Έ ( 𝑑 ) + πœ€ 𝑍 ξ€· 𝑆 β„Ž ξ€· 𝑧 0 + 𝑦 ξ€Έ ( 𝑑 ) , 𝑑 , πœ€ ξ€Έ , β„“ 𝑦 = πœƒ π‘š , 𝑑 ∈ [ 0 , 𝑏 ] , ( 2 . 2 2 ) and such that 𝑦 ( β‹… , πœ€ ) ∢ [ 0 , 𝑏 ] ⟢ ℝ 𝑛 , 𝑦 ( β‹… , πœ€ ) ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑦 ( β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑦 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , 𝑦 ( 𝑑 , 0 ) = πœƒ . ( 2 . 2 3 ) Since the vector function 𝑍 ( ( 𝑆 β„Ž 𝑧 ) ( 𝑑 ) , 𝑑 , πœ€ ) is continuously differentiable with respect to 𝑧 and continuous in πœ€ in the neighborhood of the point ( 𝑧 , πœ€ ) = ξ€· 𝑧 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ , 0 ξ€Έ , ( 2 . 2 4 ) we can separate its linear term as a function depending on 𝑦 and terms of order zero with respect to πœ€ 𝑍 ξ€· 𝑆 β„Ž ξ€· 𝑧 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ + 𝑦 ξ€Έ , 𝑑 , πœ€ ξ€Έ = 𝑓 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ + 𝐴 1 ( 𝑑 ) ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑑 ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑑 ) , 𝑑 , πœ€ ξ€Έ , ( 2 . 2 5 ) where 𝑓 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ ∢ = 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 0 ξ€Έ ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ , 𝑑 , 0 ξ€Έ , 𝑓 0 ξ€· β‹… , 𝑐 0 π‘Ÿ ξ€Έ ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝐴 1 ( 𝑑 ) = 𝐴 1 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ = πœ• 𝑍 ( 𝑆 β„Ž π‘₯ , 𝑑 , 0 ) πœ• 𝑆 β„Ž π‘₯ | | | | π‘₯ = 𝑧 0 ( 𝑑 , 𝑐 0 π‘Ÿ ) , 𝐴 1 ( β‹… ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑅 ( πœƒ , 𝑑 , 0 ) = πœƒ , πœ• 𝑅 ( πœƒ , 𝑑 , 0 ) πœ• 𝑦 = Θ , 𝑅 ( 𝑦 , β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] . ( 2 . 2 6 ) We now consider the vector function 𝑍 ( ( 𝑆 β„Ž ( 𝑧 0 + 𝑦 ) ) ( 𝑑 ) , 𝑑 , πœ€ ) in (2.22) as an inhomogeneity and we apply Theorem 1.2 to this system. As the result, we obtain the following representation for the solution of (2.22): 𝑦 ( 𝑑 , πœ€ ) = 𝑋 π‘Ÿ ( 𝑑 ) 𝑐 + 𝑦 ( 1 ) ( 𝑑 , πœ€ ) . ( 2 . 2 7 ) In this expression, the unknown vector of constants 𝑐 = 𝑐 ( πœ€ ) ∈ 𝐢 [ 0 , πœ€ 0 ] is determined from a condition similar to condition (1.37) for the existence of solution of problem (2.22): 𝐡 0 𝑐 = ξ€œ 𝑏 0 𝐻 ( 𝑠 ) ξ€Ί 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑦 ( 1 ) ξ€Έ ( 𝑠 , πœ€ ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ ξ€» 𝑑 𝑠 , ( 2 . 2 8 ) where 𝐡 0 = ξ€œ 𝑏 0 𝐻 ( 𝑠 ) 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑋 π‘Ÿ ξ€Έ ( 𝑠 ) 𝑑 𝑠 ( 2 . 2 9 ) is a 𝑑 Γ— π‘Ÿ matrix, and 𝐻 ( 𝑠 ) ∢ = 𝑃 𝑄 βˆ— 𝑑 β„“ 𝐾 ( β‹… , 𝑠 ) = 𝑃 𝑄 βˆ— 𝑑 β„“ 𝑒 𝐴 ( β‹… βˆ’ 𝜏 βˆ’ 𝑠 ) 𝜏 . ( 2 . 3 0 ) The unknown vector function 𝑦 ( 1 ) ( 𝑑 , πœ€ ) is determined by using the generalized Green operator as follows: 𝑦 ( 1 ) ( 𝑑 , πœ€ ) = πœ€ ξ€· 𝐺 ξ€Ί 𝑍 ξ€· 𝑆 β„Ž ξ€· 𝑧 0 ξ€· 𝑠 , 𝑐 0 π‘Ÿ ξ€Έ + 𝑦 ξ€Έ , 𝑠 , πœ€ ξ€Έ ξ€» ξ€Έ ( 𝑑 ) . ( 2 . 3 1 ) Let 𝑃 𝑁 ( 𝐡 0 ) be an π‘Ÿ Γ— π‘Ÿ matrix orthoprojector ℝ π‘Ÿ β†’ 𝑁 ( 𝐡 0 ) , and let 𝑃 𝑁 ( 𝐡 βˆ— 0 ) be a 𝑑 Γ— 𝑑 matrix-orthoprojector ℝ 𝑑 β†’ 𝑁 ( 𝐡 βˆ— 0 ) . Equation (2.28) is solvable with respect to 𝑐 ∈ ℝ π‘Ÿ if and only if 𝑃 𝑁 ( 𝐡 βˆ— 0 ) ξ€œ 𝑏 0 𝐻 ( 𝑠 ) ξ€Ί 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑦 ( 1 ) ξ€Έ ( 𝑠 , πœ€ ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ ξ€» 𝑑 𝑠 = πœƒ 𝑑 . ( 2 . 3 2 ) For 𝑃 𝑁 ( 𝐡 βˆ— 0 ) = Θ 𝑑 , ( 2 . 3 3 ) the last condition is always satisfied and (2.28) is solvable with respect to 𝑐 ∈ ℝ π‘Ÿ up to an arbitrary vector constant 𝑃 𝑁 ( 𝐡 0 ) 𝑐 ∈ ℝ π‘Ÿ from the null space of the matrix 𝐡 0 𝑐 = 𝐡 + 0 ξ€œ 𝑏 0 𝐻 ( 𝑠 ) ξ€Ί 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑦 ( 1 ) ξ€Έ ( 𝑠 , πœ€ ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ ξ€» 𝑑 𝑠 + 𝑃 𝑁 ( 𝐡 0 ) 𝑐 . ( 2 . 3 4 ) To find a solution 𝑦 = 𝑦 ( 𝑑 , πœ€ ) of (2.28) such that 𝑦 ( β‹… , πœ€ ) ∢ [ 0 , 𝑏 ] ⟢ 𝑅 𝑛 , 𝑦 ( β‹… , πœ€ ) ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑦 ( β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑦 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , 𝑦 ( 𝑑 , 0 ) = πœƒ , ( 2 . 3 5 ) it is necessary to solve the following operator system: 𝑦 ( 𝑑 , πœ€ ) = 𝑋 π‘Ÿ ( 𝑑 ) 𝑐 + 𝑦 ( 1 ) ( 𝑑 , πœ€ ) , 𝑐 = 𝐡 + 0 ξ€œ 𝑏 0 𝐻 ( 𝑠 ) ξ€Ί 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑦 ( 1 ) ξ€Έ ( 𝑠 , πœ€ ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ ξ€» 𝑑 𝑠 , 𝑦 ( 1 ) ( 𝑑 , πœ€ ) = πœ€ 𝐺 ξ€Ί 𝑍 ξ€· 𝑆 β„Ž ξ€· 𝑧 0 ξ€· 𝑠 , 𝑐 0 π‘Ÿ ξ€Έ + 𝑦 ξ€Έ , 𝑠 , πœ€ ξ€Έ ξ€» ( 𝑑 ) . ( 2 . 3 6 ) The operator system (2.36) belongs to the class of systems solvable by the method of simple iterations, convergent for sufficiently small πœ€ ∈ [ 0 , πœ€ 0 ] (see [10, page 188]). Indeed, system (2.36) can be rewritten in the form 𝑒 = 𝐿 ( 1 ) 𝑒 + 𝐹 𝑒 , ( 2 . 3 7 ) where 𝑒 = c o l ( 𝑦 ( 𝑑 , πœ€ ) , 𝑐 ( πœ€ ) , 𝑦 ( 1 ) ( 𝑑 , πœ€ ) ) is a ( 2 𝑛 + π‘Ÿ ) -dimensional column vector, 𝐿 ( 1 ) is a linear operator 𝐿 ( 1 ) ∢ = βŽ› ⎜ ⎜ ⎜ ⎝ Θ 𝑋 π‘Ÿ ( 𝑑 ) 𝐼 Θ π‘Ÿ , 𝑛 Θ π‘Ÿ , π‘Ÿ 𝐿 1 Θ Θ 𝑛 , π‘Ÿ Θ ⎞ ⎟ ⎟ ⎟ ⎠ , ( 2 . 3 8 ) where 𝐿 1 ( βˆ— ) = 𝐡 + 0 ξ€œ 𝑏 0 𝐻 ( 𝑠 ) 𝐴 1 ( 𝑠 ) ( βˆ— ) 𝑑 𝑠 , ( 2 . 3 9 ) and 𝐹 is a nonlinear operator 𝐹 𝑒 ∢ = βŽ› ⎜ ⎜ ⎜ ⎜ ⎝ πœƒ 𝐡 + 0 ξ€œ 𝑏 0 𝐻 ( 𝑠 ) 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ 𝑑 𝑠 πœ€ ξ€· 𝐺 ξ€Ί 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 0 ξ€Έ ξ€· 𝑠 , 𝑐 0 π‘Ÿ ξ€Έ , 𝑠 , 0 ξ€Έ + 𝐴 1 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) + 𝑅 ξ€· ξ€· 𝑆 β„Ž 𝑦 ξ€Έ ( 𝑠 , πœ€ ) , 𝑠 , πœ€ ξ€Έ ξ€» ξ€Έ ( 𝑑 ) ⎞ ⎟ ⎟ ⎟ ⎟ ⎠ . ( 2 . 4 0 ) In view of the structure of the operator 𝐿 ( 1 ) containing zero blocks on and below the main diagonal, the inverse operator ξ€· 𝐼 2 𝑛 + π‘Ÿ βˆ’ 𝐿 ( 1 ) ξ€Έ βˆ’ 1 ( 2 . 4 1 ) exists. System (2.37) can be transformed into 𝑒 = 𝑆 𝑒 , ( 2 . 4 2 ) where 𝑆 ∢ = ξ€· 𝐼 2 𝑛 + π‘Ÿ βˆ’ 𝐿 ( 1 ) ξ€Έ βˆ’ 1 𝐹 ( 2 . 4 3 ) is a contraction operator in a sufficiently small neighborhood of the point ( 𝑧 , πœ€ ) = ξ€· 𝑧 0 ξ€· 𝑑 , 𝑐 0 π‘Ÿ ξ€Έ , 0 ξ€Έ . ( 2 . 4 4 ) Thus, the solvability of the last operator system can be established by using one of the existing versions of the fixed-point principles [12] applicable to the system for sufficiently small πœ€ ∈ [ 0 , πœ€ 0 ] . It is easy to prove that the sufficient condition 𝑃 𝑁 ( 𝐡 βˆ— 0 ) = Θ 𝑑 for the existence of solutions of the boundary value problem (2.13), (2.14) means that the constant 𝑐 0 π‘Ÿ ∈ ℝ π‘Ÿ of the equation for generating constant (2.17) is a simple root of equation (2.17) [2]. By using the method of simple iterations, we can find the solution of the operator system and hence the solution of the original boundary value problem (2.13), (2.14). Now, we arrive at the following theorem.

Theorem 2.3 (sufficient condition). Assume that the boundary value problem (2.13), (2.14) satisfies the conditions listed above and the corresponding linear boundary value problem (1.25) has an π‘Ÿ -dimensional family of linearly independent solutions of the form (1.38). Then, for any simple root 𝑐 π‘Ÿ = 𝑐 0 π‘Ÿ ∈ ℝ π‘Ÿ of the equation for generating the constants (2.17), there exist at least one solution of the boundary value problem (2.13), (2.14). The indicated solution 𝑧 ( 𝑑 , πœ€ ) is such that 𝑧 ( β‹… , πœ€ ) ∈ 𝐷 𝑝 [ 0 , 𝑏 ] , Μ‡ 𝑧 ( β‹… , πœ€ ) ∈ 𝐿 𝑝 [ 0 , 𝑏 ] , 𝑧 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , ( 2 . 4 5 ) and, for πœ€ = 0 , turns into one of the generating solutions (1.38) with a constant 𝑐 0 π‘Ÿ ∈ ℝ π‘Ÿ ; that is, 𝑧 ( 𝑑 , 0 ) = 𝑧 0 ( 𝑑 , 𝑐 0 π‘Ÿ ) . This solution can be found by the method of simple iterations, which is convergent for a sufficiently small πœ€ ∈ [ 0 , πœ€ 0 ] .

Corollary 2.4. If the number 𝑛 of unknown variables is equal to the number π‘š of boundary conditions (and hence π‘Ÿ = 𝑑 ), the boundary value problem (2.13), (2.14) has a unique solution. In such a case, the problems considered for functional-differential equations are of Fredholm's type with a zero index. By using the procedure proposed in [2] with some simplifying assumptions, we can generalize the proposed method to the case of multiple roots of equation (2.17) to determine sufficient conditions for the existence of solutions of the boundary-value problem (2.13), (2.14).

3. Example

We will illustrate the above proved theorems on the example of a weakly perturbed linear boundary value problem. Consider the following simplest boundary value problem-a periodic problem for the delayed differential equation: Μ‡ 𝑧 ( 𝑑 ) = 𝑧 ( 𝑑 βˆ’ 𝜏 ) + πœ€ π‘˜  𝑖 = 1 𝐡 𝑖 ( 𝑑 ) 𝑧 ξ€· β„Ž 𝑖 ( 𝑑 ) ξ€Έ + 𝑔 ( 𝑑 ) , 𝑑 ∈ ( 0 , 𝑇 ] , 𝑧 ( 𝑠 ) = πœ“ ( 𝑠 ) , i f 𝑠 < 0 , 𝑧 ( 0 ) = 𝑧 ( 𝑇 ) , ( 3 . 1 ) where 0 < 𝜏 , 𝑇 = c o n s t , 𝐡 𝑖 are 𝑛 Γ— 𝑛 matrices, 𝐡 𝑖 , 𝑔 ∈ 𝐿 𝑝 [ 0 , 𝑇 ] , πœ“ ∢ ℝ 1 ⧡ ( 0 , 𝑇 ] β†’ ℝ 𝑛 , β„Ž 𝑖 ( 𝑑 ) ≀ 𝑑 are measurable functions. Using the symbols 𝑆 β„Ž 𝑖 and πœ“ β„Ž 𝑖 (see (1.3), (1.4), (2.9)), we arrive at the following operator system: Μ‡ 𝑧 ( 𝑑 ) = 𝑧 ( 𝑑 βˆ’ 𝜏 ) + πœ€ 𝐡 ( 𝑑 ) ξ€· 𝑆 β„Ž 𝑧 ξ€Έ ( 𝑑 ) + πœ‘ ( 𝑑 , πœ€ ) , β„“ 𝑧 ∢ = 𝑧 ( 0 ) βˆ’ 𝑧 ( 𝑇 ) = πœƒ 𝑛 , ( 3 . 2 ) where 𝐡 ( 𝑑 ) ∢ = ( 𝐡 1 ( 𝑑 ) , … , 𝐡 π‘˜ ( 𝑑 ) ) is an 𝑛 Γ— 𝑁 matrix ( 𝑁 = 𝑛 π‘˜ ) , and πœ‘ ( 𝑑 , πœ€ ) ∢ = 𝑔 ( 𝑑 ) + πœ“ β„Ž 0 ( 𝑑 ) + πœ€ π‘˜  𝑖 = 1 𝐡 𝑖 ( 𝑑 ) πœ“ β„Ž 𝑖 ( 𝑑 ) ∈ 𝐿 𝑝 [ 0 , 𝑇 ] . ( 3 . 3 ) We will consider the simplest case with 𝑇 ≀ 𝜏 . Utilizing the delayed matrix exponential, it can be easily verified that in this case, the matrix 𝑋 ( 𝑑 ) = 𝑒 𝐼 ( 𝑑 βˆ’ 𝜏 ) 𝜏 = 𝐼 ( 3 . 4 ) is a normal fundamental matrix for the homogeneous generating system Μ‡ 𝑧 ( 𝑑 ) = 𝑧 ( 𝑑 βˆ’ 𝜏 ) . ( 3 . 5 ) Then, 𝑄 ∢ = β„“ 𝑋 ( β‹… ) = 𝑒 βˆ’ 𝐼 𝜏 𝜏 βˆ’ 𝑒 𝐼 ( 𝑇 βˆ’ 𝜏 ) 𝜏 = πœƒ 𝑛 , 𝑃 𝑄 = 𝑃 𝑄 βˆ— = 𝐼 , ( π‘Ÿ = 𝑛 , 𝑑 = π‘š = 𝑛 ) , 𝐾 ( 𝑑 , 𝑠 ) = ⎧ βŽͺ ⎨ βŽͺ ⎩ 𝑒 𝐼 ( 𝑑 βˆ’ 𝜏 βˆ’ 𝑠 ) 𝜏 = 𝐼 , 0 ≀ 𝑠 ≀ 𝑑 ≀ 𝑇 , Θ , 𝑠 > 𝑑 , β„“ 𝐾 ( β‹… , 𝑠 ) = 𝐾 ( 0 , 𝑠 ) βˆ’ 𝐾 ( 𝑇 , 𝑠 ) = βˆ’ 𝐼 , 𝐻 ( 𝜏 ) = 𝑃 𝑄 βˆ— β„“ 𝐾 ( β‹… , 𝑠 ) = βˆ’ 𝐼 , ξ€· 𝑆 β„Ž 𝑖 𝐼 ξ€Έ ( 𝑑 ) = πœ’ β„Ž 𝑖 ( 𝑑 , 0 ) β‹… 𝐼 = 𝐼 β‹… ⎧ βŽͺ ⎨ βŽͺ ⎩ 1 , i f 0 ≀ β„Ž 𝑖 ( 𝑑 ) ≀ 𝑇 , 0 , i f β„Ž 𝑖 ( 𝑑 ) < 0 . ( 3 . 6 ) To illustrate the theorems proved above, we will find the conditions for which the boundary value problem (3.1) has a solution 𝑧 ( 𝑑 , πœ€ ) that, for πœ€ = 0 , turns into one of solutions (1.38) 𝑧 0 ( 𝑑 , 𝑐 π‘Ÿ ) of the generating problem. In contrast to the previous works [7, 9], we consider the case when the unperturbed boundary-value problem Μ‡ 𝑧 ( 𝑑 ) = 𝑧 ( 𝑑 βˆ’ 𝜏 ) + πœ‘ ( 𝑑 , 0 ) , 𝑧 ( 0 ) = 𝑧 ( 𝑇 ) ( 3 . 7 ) has an 𝑛 -parametric family of linear-independent solutions of the form(1.38) 𝑧 ∢ = 𝑧 0 ξ€· 𝑑 , 𝑐 𝑛 ξ€Έ = 𝑐 𝑛 + ( 𝐺 πœ‘ ) ( 𝑑 ) , βˆ€ 𝑐 𝑛 ∈ ℝ 𝑛 . ( 3 . 8 ) For this purpose, it is necessary and sufficient for the vector function πœ‘ ( 𝑑 ) = 𝑔 ( 𝑑 ) + πœ“ β„Ž 0 ( 𝑑 ) ( 3 . 9 ) to satisfy the condition of type (1.37) ξ€œ 𝑇 0 𝐻 ( 𝑠 ) πœ‘ ( 𝑠 ) 𝑑 𝑠 = βˆ’ ξ€œ 𝑇 0 πœ‘ ( 𝑠 ) 𝑑 𝑠 = πœƒ 𝑛 . ( 3 . 1 0 ) Then, according to the Theorem 2.1, the constant 𝑐 𝑛 = 𝑐 0 𝑛 ∈ ℝ 𝑛 must satisfy (2.17), that is, the equation 𝐹 ξ€· 𝑐 0 𝑛 ξ€Έ ∢ = ξ€œ 𝑇 0 𝐻 ( 𝑠 ) 𝑍 ξ€· ξ€· 𝑆 β„Ž 𝑧 0 ξ€Έ ξ€· 𝑠 , 𝑐 0 𝑛 ξ€Έ , 𝑠 , 0 ξ€Έ 𝑑 𝑠 = πœƒ 𝑛 , ( 3 . 1 1 ) which in our case is a linear algebraic system 𝐡 0 𝑐 0 𝑛 = βˆ’ ξ€œ 𝑇 0 𝐡 ( 𝑠 ) ξ€· 𝑆 β„Ž ( 𝐺 πœ‘ ) ξ€Έ ( 𝑠 ) 𝑑 𝑠 , ( 3 . 1 2 ) with the 𝑛 Γ— 𝑛 matrix 𝐡 0 in the form 𝐡 0 = ξ€œ 𝑇 0 𝐻 ( 𝑠 ) 𝐡 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝐼 ξ€Έ ( 𝑠 ) 𝑑 𝑠 = βˆ’ ξ€œ 𝑇 0 π‘˜  𝑖 = 1 𝐡 𝑖 ( 𝑠 ) ξ€· 𝑆 β„Ž 𝑖 𝐼 ξ€Έ ( 𝑠 ) 𝑑 𝑠 = βˆ’ π‘˜  𝑖 = 1 ξ€œ 𝑇 0 𝐡 𝑖 ( 𝑠 ) πœ’ β„Ž 𝑖 ( 𝑠 , 0 ) 𝑑 𝑠 . ( 3 . 1 3 ) According to Corollary 2.4, if d e t 𝐡 0 β‰  0 , the problem (3.1) for the case 𝑇 ≀ 𝜏 has a unique solution 𝑧 ( 𝑑 , πœ€ ) with the properties 𝑧 ( β‹… , πœ€ ) ∈ 𝐷 𝑛 𝑝 [ 0 , 𝑇 ] , Μ‡ 𝑧 ( β‹… , πœ€ ) ∈ 𝐿 𝑛 𝑝 [ 0 , 𝑇 ] , 𝑧 ( 𝑑 , β‹… ) ∈ 𝐢 ξ€Ί 0 , πœ€ 0 ξ€» , 𝑧 ( 𝑑 , 0 ) = 𝑧 0 ξ€· 𝑑 , 𝑐 0 𝑛 ξ€Έ , ( 3 . 1 4 ) for 𝑔 ∈ 𝐿 𝑝 [ 0 , 𝑇 ] , πœ“ ( 𝑑 ) ∈ 𝐿 𝑝 [ 0 , 𝑇 ] , and for measurable delays β„Ž 𝑖 that which satisfy the criterion (3.10) of the existence of a generating solution where 𝑐 0 𝑛 = βˆ’ 𝐡 + 0 ξ€œ 𝑇 0 𝐡 ( 𝑠 ) ξ€· 𝑆 β„Ž ( 𝐺 πœ‘ ) ξ€Έ ( 𝑠 ) 𝑑 𝑠 . ( 3 . 1 5 ) A solution 𝑧 ( 𝑑 , πœ€ ) of the boundary value problem (3.1) can be found by the convergent method of simple iterations (see Theorem 2.3).

If, for example, β„Ž 𝑖 ( 𝑑 ) = 𝑑 βˆ’ Ξ” 𝑖 , where 0 < Ξ” 𝑖 = c o n s t < 𝑇 , 𝑖 = 1 , … , π‘˜ , then πœ’ β„Ž 𝑖 ( 𝑑 , 0 ) = ⎧ βŽͺ ⎨ βŽͺ ⎩ 1 i f 0 ≀ β„Ž 𝑖 ( 𝑑 ) = 𝑑 βˆ’ Ξ” 𝑖 ≀ 𝑇 , 0 i f β„Ž 𝑖 ( 𝑑 ) = 𝑑 βˆ’ Ξ” 𝑖 < 0 , = ⎧ βŽͺ ⎨ βŽͺ ⎩ 1 i f Ξ” 𝑖 ≀ 𝑑 ≀ 𝑇 + Ξ” 𝑖 , 0 , i f 𝑑 < Ξ” 𝑖 . ( 3 . 1 6 ) The 𝑛 Γ— 𝑛 matrix 𝐡 0 can be rewritten in the form 𝐡 0 = ξ€œ 𝑇 0 𝐻 ( 𝑠 ) π‘˜  𝑖 = 1 𝐡 𝑖 ( 𝑠 ) πœ’ β„Ž 𝑖 ( 𝑠 , 0 ) 𝑑 𝜏 = βˆ’ π‘˜  𝑖 = 1 ξ€œ 𝑇 0 𝐡 𝑖 ( 𝑠 ) πœ’ β„Ž 𝑖 ( 𝑠 , 0 ) 𝑑 𝑠 = βˆ’ π‘˜  𝑖 = 1 ξ€œ 𝑇 Ξ” 𝑖 𝐡 𝑖 ( 𝑠 ) 𝑑 𝑠 , ( 3 . 1 7 ) and the unique solvability condition of the boundary value problem (3.1) takes the form d e t ⎑ ⎒ ⎣ π‘˜  𝑖 = 1 ξ€œ 𝑇 Ξ” 𝑖 𝐡 𝑖 ( 𝑠 ) 𝑑 𝑠 ⎀ βŽ₯ ⎦ β‰  0 . ( 3 . 1 8 ) It is easy to see that if the vector function 𝑍 ( 𝑧 ( β„Ž 𝑖 ( 𝑑 ) ) , 𝑑 , πœ€ ) is nonlinear in 𝑧 , for example as a square, then (3.11) generating the constants will be a square-algebraic system and, in this case, the boundary value problem (3.1) can have two solutions branching from the point πœ€ = 0 .

Acknowledgments

The first and the fourth authors were supported by the Grant no. 1/0090/09 of the Grant Agency of Slovak Republic (VEGA) and by the project APVV-0700-07 of Slovak Research and Development Agency. The second author was supported by the Grant no. P201/11/0768 of Czech Grant Agency, by the Council of Czech Government MSM 0021630503 and by the Project FEKT/FSI-S-11-1-1159. The third author was supported by the Project no. M/34-2008 of Ukrainian Ministry of Education, Ukraine.

References

  1. N. V. Azbelev and V. P. Maksimov, β€œEquations with delayed arguments,” Differential Equations, vol. 18, pp. 1419–1441, 1983, translated in Differentsial'nye Uravneniya, Vol. 18, no. 12, pp. 2027-2050.
  2. A. A. Boichuk and A. M. Samoilenko, Generalized Inverse Operators and Fredholm Boundary Value Problems, VSP, Boston, Mass, USA, 2004.
  3. J. Hale, Theory of Functional Differential Equations, Applied Mathematical Sciences, Vol. 3, Springer, New York, NY, USA, 2nd edition, 1977.
  4. J. Mallet-Paret, β€œThe fredholm alternative for functional-differential equations of mixed type,” Journal of Dynamics and Differential Equations, vol. 11, no. 1, pp. 1–47, 1999. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  5. D. Ya. Khusainov and G. V. Shuklin, β€œRelative controllability in systems with pure delay,” International Applied Mechanics, vol. 41, no. 2, pp. 210–221, 2005, translated in Prikladnaya Mekhanika, Vol. 41, no. 2, pp. 118-130. View at Publisher Β· View at Google Scholar Β· View at MathSciNet
  6. A. A. Boichuk, J. Diblík, D. Khusainov, and M. Růžičková, β€œFredholm's boundary-value problems for differential systems with a single delay,” Nonlinear Analysis, Theory, Methods and Applications, vol. 72, no. 5, pp. 2251–2258, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  7. A. A. Boichuk, J. Diblík, D. Y. A. Khusainov, and M. Růžičková, β€œBoundary value problems for delay differential systems,” Advances in Difference Equations, vol. 2010, Article ID 593834, 10 pages, 2010. View at Zentralblatt MATH
  8. M. Medvĕd, M. Pospíšil, and L. Škripková, β€œStability and the nonexistence of blowing-up solutions of nonlinear delay systems with linear parts defined by permutable matrices,” Nonlinear Analysis, Theory, Methods and Applications, vol. 74, no. 12, pp. 3903–3911, 2011. View at Publisher Β· View at Google Scholar
  9. A. A. Boichuk and M. K. Grammatikopoulos, β€œPerturbed fredholm boundary value problems for delay differential systems,” Abstract and Applied Analysis, no. 15, pp. 843–864, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  10. A. A. Boichuk, V. F. Zhuravlev, and A. M. Samoilenko, β€œGeneralized inverse operators and noether boundary-value problems,” in Proceedings of the Institute of Mathematics of the National Academy of Sciences of the Ukraine, vol. 13, p. 320, Kyiv, Ukrania, 1995.
  11. I. G. Malkin, Some Problems in the Theory of Nonlinear Oscillations, Gostekhizdat, Moscow, Russia, 1956.
  12. M. A. Krasnosel'skii, G. M. Vainikko, P. P. Zabreiko, Y. A. B. Rutitskii, and V. Ya. Stetsenko, Approximate Solution of Operator Equations, Nauka, Moscow, Russia, 1968, Translation: Noordhoff, Groningen, 1972.