Normal Forms for Nonautonomous Differential Equations

We extend Henry Poincarés normal form theory for autonomous differential equations x=f(x) to nonautonomous differential equations x=f(t, x). Poincarés nonresonance condition λ_j−∑ⁿ_i=1 ?_iλ_i≠0 for eigenvalues is generalized to the new nonresonance condition λ_j∩∑ⁿ_i=1 ?_iλ_i=∅ for spectral intervals.     

The improved linear multistep methods for differential equations with piecewise continuous arguments

This paper deals with the convergence of the linear multistep methods for the equation x′(t) = ax(t) + a₀x([t]). Numerical experiments demonstrate that the 2-step Adams-Bashforth method is only of order p = 0 when applied to the given equation. An improved linear multistep methods is constructed. It is proved that these methods preserve their original convergence order for ordinary differential equations (ODEs) and some numerical experiments are given.     

An approximation property of exponential functions

We solve the inhomogeneous linear first order differential equations of the form y′(x) ? λy(x) = Σ _m=0 ^∞ a _m (x ? c) ^m , and prove an approximation property of exponential functions. More precisely, we prove the local Hyers-Ulam stability of linear first order differential equations of the form y′(x) = λy(x) in a special class of analytic functions.     

Linear Differential Equations with Variable Coefficients in Regular and Some Singular Cases

In this paper the asymptotic properties as t → + ∞ for a single linear differential equation of the form x⁽ⁿ⁾ + a₁ (t)x^(n?1)+…. + a_n(t)x = 0, where the coefficients aj (z) are supposed to be of the power order of growth, are considered. The results obtained in the previous publications of the author were related to the so called regular case when a complete set of roots {λ,(t)}, j = 1, 2, …, n of the characteristic polynomial yⁿ + a₁ (t)y^n?1 + … + a_n(t) possesses the property of asymptotic separability. One of the main restrictions of the regular case consists of the demand that the roots of the set {λ,(t)} have not to be equivalent in pairs for t → + ∞. In this paper we consider the some more general case when the set of characteristic roots possesses the property of asymptotic independence which includes the case when the roots may be equivdent in pairs. But some restrictions on the asymptotic behaviour of their differences λ_i(t)→ λ_j(t) are preserved. This case demands more complicated technique of investigation. For this purpose the so called asymptotic spaces were introduced. The theory of asymptotic spaces is used for formal solution of an operator equation of the form x = A(x) and has the analogous meaning as the classical theory of solving this equation in Band spaces. For the considered differential equation, the main asymptotic terms of a fundamental system of solution is given in a simple explicit form and the asymptotic fundamental system is represented in the form of asymptotic Emits for several iterate sequences.     

On linear systems with quasiperiodic coefficients and bounded solutions

For a discrete dynamical system ω _n =ω₀+αn, where a is a constant vector with rationally independent coordinates, on thes-dimensional torus Ω we consider the setL of its linear unitary extensionsx _n+1=A(ω₀+αn)x _n , whereA (Ω) is a continuous function on the torus Ω with values in the space ofm-dimensional unitary matrices. It is proved that linear extensions whose solutions are not almost periodic form a set of the second category inL (representable as an intersection of countably many everywhere dense open subsets). A similar assertion is true for systems of linear differential equations with quasiperiodic skew-symmetric matrices.     

Asymptotical Formulas For Solutions Of Linear Differential Systems Of The First Order

In this paper a system of differential equations y′ ? A(·,λ)y = 0 is considered on the finite interval [a,b] where λ ∈ C, A(·, λ):= λ A₁+ A ₀ +λ _?1A_?1(·,λ) and A ₁,A ₀, A _{? 1} are n × n matrix-functions. The main assumptions: A ₁ is absolutely continuous on the interval [a, b], A ₀ and A _{- 1}(·,λ) are summable on the same interval when ¦λ¦ is sufficiently large; the roots φ₁(x),…,φ_n (x) of the characteristic equation det (φ E — A ₁) = 0 are different for all x ∈ [a,b] and do not vanish; there exists some unlimited set Ω ? C on which the inequalities Re(λφ₁(x)) ≤ … ≤ Re (λφ_n(x)) are fulfilled for all x ∈ [a,b] and for some numeration of the functions φ_j(x). The asymptotic formula of the exponential type for a fundamental matrix of solutions of the system is obtained for sufficiently large ¦λ¦. The remainder term of this formula has a new type dependence on properties of the coefficients A ₁ (x), A _o (x) and A _{- 1} (x).     

Lipschitz lower sigma-exponents of linear differential systems

For the lower sigma-exponent of the linear differential system ? = A(t)x, x ∈ R ⁿ , t ≥ 0, defined by the formula Δ_σ(A) ≡ inf_λ[Q]≤-σ λ ₁(A + Q), σ > 0, on the basis of the lower characteristic exponents λ ₁(A+Q) of perturbed linear systems with Lyapunov exponents λ[Q] ≤ ?σ < 0 of perturbations Q, we prove the following general form as a function of the parameter σ > 0. For any nondecreasing bounded function f(σ) of the parameter σ ∈ (0,+∞) that coincides with a constant on some infinite interval (σ ₀,+∞), σ ₀ ≥ 0, and satisfies the Lipschitz condition on the complementary interval (0, σ ₀], we prove the existence of a linear system with coefficient matrix A _f (t) bounded on the half-line [0,+∞) whose lower sigma-exponent Δ_σ(A _f ) coincides with the function f(σ) on the entire interval (0,+∞).     

Existence of pseudo-almost automorphic solutions for nonlinear differential equations

In this paper, we discuss the existence of pseudo-almost automorphic solutions to linear differential equation which has an exponential trichotomy~ and the results also hold for some nonlinear equations with the form x＇（t） = f（t,x（t）） ＋ λg（t,x（t））, where f,g are pseudo-almost automorphic functions. We prove our main result by the application of Leray-Schauder fixed point theorem.     

A combinatorial problem involving graphs and matrices

In this paper we discuss a combinatorial problem involving graphs and matrices. Our problem is a matrix analogue of the classical problem of finding a system of distinct representatives (transversal) of a family of sets and relates closely to an extremal problem involving 1-factors and a long standing conjecture in the dimension theory of partially ordered sets. For an integer n ?1, let n denote the n element set {1,2,3,…, n}. Then let A be a k×t matrix. We say that A satisfies property P(n, k) when the following condition is satisfied: For every k-taple (x₁,x₂,…,x_k?n^k there exist k distinct integers j₁,j₂,…,j_k so that x_i= a_ii for i= 1,2,…,k. The minimum value of t for which there exists a k × t matrix A satisfying property P(n,k) is denoted by f(n,k). For each k?1 and n sufficiently large, we give an explicit formula for f(n, k): for each n?1 and k sufficiently large, we use probabilistic methods to provide inequalities for f(n,k).     

Dichotomies and moving singularities

We consider a linear differential system ε ^σ Φ (t,ε)Y′ =A(t, ε)Y, with ε a small parameter and Φ(t, ε) a function which may vanish in the domain of definition. Under some conditions imposed on the eigenvalues of the matrixA(t, ε), there exists an invertible matrixH(t, ε) which is continuous on ([0,a] × [0, ε₀]). The transformationY=H(t, ε)Z takes then dimensional linear system into two differential systems of orderk andn?k respectively, withk. Thus the investigaton ofn dimensional systems encountered in singular perturbation as well as in stability theory is considerably simplified.     

Note on forced oscillation of nth-order sublinear differential equations

In the case of oscillatory potentials, we establish an oscillation theorem for the forced sublinear differential equation x(n)+q(t)λ|x|sgnx=e(t), t∈[t₀,∞). No restriction is imposed on the forcing term e(t) to be the nth derivative of an oscillatory function. In particular, we show that all solutions of the equation x^″+tαsintλ|x|sgnx=mtβcost, t?0, 0<λ<1 are oscillatory for all m≠0 if β>(α+2)/(1−λ). This provides an analogue of a result of Nasr [Proc. Amer. Math. Soc. 126 (1998) 123] for the forced superlinear equation and answers a question raised in an earlier paper [J.S.W. Wong, SIAM J. Math. Anal. 19 (1988) 673].     

Small deviations of series of independent positive random variables with weights close to exponential

Let ξ, ξ₀, ξ₁, ... be independent identically distributed (i.i.d.) positive random variables. The present paper is a continuation of the article [1] in which the asymptotics of probabilities of small deviations of series S = Σ _j=0 ^∞ a(j)ξ _j was studied under different assumptions on the rate of decrease of the probability ?(ξ < x) as x → 0, as well as of the coefficients a(j) ≥ 0 as j → ∞. We study the asymptotics of ?(S < x) as x → 0 under the condition that the coefficients a(j) are close to exponential. In the case when the coefficients a(j) are exponential and ?(ξ < x) ～ bx ^α as x → 0, b > 0, a > 0, the asymptotics ?(S < x) is obtained in an explicit form up to the factor x ^o(1). Originality of the approach of the present paper consists in employing the theory of delayed differential equations. This approach differs significantly from that in [1].     

A new reduction method in integer programming

The system of two equations in nonnegative integer unknowns x_j and expressed in the form ∑ⁿ_j=1 a_1jx_j=b₁, ∑ⁿ_j=1 a_2jx_j=b₂ has equivalent solutions to the single equation ∑ⁿ_j=1 (t₁a_ij+t₂a_2j)a_j= t₁b₁+t₂b₂, provided t₁ and t₂ are suitably chosen. Glover [2] has n inequalities for determining t₁ and t₂. Elimam and Elmaghraby [1] try to achieve a single inequality for t₁ and t₂. We show, however, that the inequality of [1] may give t₁ and t₂ values that do not produce equivalence. We present a new theorem which leads to a single inequality constraint giving t₁ and t₂ values that consistently produce equivalence.     

Linearized stability analysis of discrete Volterra equations

Various different types of stability are defined, in a unified framework, for discrete Volterra equations of the type x(n)=f(n)+∑ⁿ_j=0K(n,j,x(n)) (n?0). Under appropriate assumptions, stability results are obtainable from those valid in the linear case (K(n,j,x(n))=B(n,j)x(j)), and a linearized stability theory is studied here by using the fundamental and resolvent matrices. Several necessary and sufficient conditions for stability are obtained for solutions of the linear equation by considering the equations in various choices of Banach space , the elements of which are sequences of vectors (, , n,j?0, etc.). We show that the theory, including a number of new results as well as results already known, can be presented in a systematic framework, in which results parallel corresponding results for classical Volterra integral equations of the second kind.     

A RIGOROUS DERIVATION OF THE GROSS-PITAEVSKII HIERARCHY FOR WEAKLY COUPLED TWO-DIMENSIONAL BOSONS

In this article, we consider the dynamics of N two-dimensional boson systems interacting through a pair potential N-1Va(xi-xj) where Va(x) = a-2V (x/a). It is well known that the Gross-Pitaevskii (GP) equation is a nonlinear Schrdinger equation and the GP hierarchy is an infinite BBGKY hierarchy of equations so that if ut solves the GP equation, then the family of k-particle density matrices {k ut, k ≥ 1} solves the GP hierarchy. Denote by ψN,t the solution to the N-particle Schrdinger equation. Under the assumption that a = N-ε for 0 ε 3/4, we prove that as N →∞ the limit points of the k-particle density matrices of ψN,t are solutions of the GP hierarchy with the coupling constant in the nonlinear term of the GP equation given by ∫V (x) dx.     

Finite difference methods for two-point boundary value problems involving high order differential equations

We discuss the construction of finite difference schemes for the two-point nonlinear boundary value problem:y ⁽²ⁿ⁾+f(x,y)=0,y ^(2j)(a)=A _2j ,y ^(2j)(b)=B _2j ,j=0(1)n–1,n2. In the case of linear differential equations, these finite difference schemes lead to (2n+1)-diagonal linear systems. We consider in detail methods of orders two, four and six for two-point boundary value problems involving a fourth order differential equation; convergence of these methods is established and illustrated by numerical examples.     

On the perturbation of the observability equation in linear control systems

The aim of this paper is to investigate stability and sensitivity of the observability variable in linear control systems, (LCS) for short. We first present two results of Hölder continuity in the abstract framework of the ordinary differential equation initial-value problem x′(t) = f(t,x(t)),x(t ₀) = x ₀. Afterwards, we apply our results to automatic systems, providing henceforth the sharpest bounds for the parametric input-output relation in LCS.     

On oscillatory solutions of certain forced Emden-Fowler like equations

We give a constructive proof of existence to oscillatory solutions for the differential equations x^″(t)+a(t)λ|x(t)|sign[x(t)]=e(t), where t?t₀?1 and λ>1, that decay to 0 when t→+∞ as O(t−μ) for μ>0 as close as desired to the “critical quantity” . For this class of equations, we have limt→+∞E(t)=0, where E(t)<0 and E^″(t)=e(t) throughout [t₀,+∞). We also establish that for any μ>μ^? and any negative-valued E(t)=o(t−μ) as t→+∞ the differential equation has a negative-valued solution decaying to 0 at + ∞ as o(t−μ). In this way, we are not in the reach of any of the developments from the recent paper [C.H. Ou, J.S.W. Wong, Forced oscillation of nth-order functional differential equations, J. Math. Anal. Appl. 262 (2001) 722-732].     

Oscillation results for selfadjoint differential systems

We consider the second-order differential system, (1) (R(t) Y′)′ + Q(t) Y = 0, where R, Q, Y are n × n matrices with R(t), Q(t) symmetric and R(t) positive definite for t ϵ [a, + ∞) (R(t) > 0, t ⩾ a). We establish sufficient conditions in order that all prepared solutions Y(t) of (1) are oscillatory; that is, det Y(t) vanishes infinitely often on [a, + ∞). The conditions involve the smallest and largest eigen-values λ_n(R⁻¹(t)) and λ₁(∝_a^t Q(s) ds), respectively. The results obtained can be regarded as generalizing well-known results of Leighton and others in the scalar case.