A one-step 7-stage Hermite-Birkhoff-Taylor ODE solver of order 11

A one-step 7-stage Hermite-Birkhoff-Taylor method of order 11, denoted by HBT(11)7, is constructed for solving nonstiff first-order initial value problems y^′=f(t,y), y(t₀)=y₀. The method adds the derivatives y^′ to y⁽⁶⁾, used in Taylor methods, to a 7-stage Runge-Kutta method of order 6. Forcing an expansion of the numerical solution to agree with a Taylor expansion of the true solution to order 11 leads to Taylor- and Runge-Kutta-type order conditions. These conditions are reorganized into Vandermonde-type linear systems whose solutions are the coefficients of the method. The new method has a larger scaled interval of absolute stability than the Dormand-Prince DP87 and a larger unscaled interval of absolute stability than the Taylor method, T11, of order 11. HBT(11)7 is superior to DP87 and T11 in solving several problems often used to test higher-order ODE solvers on the basis of the number of steps, CPU time, and maximum global error. Numerical results show the benefit of adding high-order derivatives to Runge-Kutta methods.     

One-step 9-stage Hermite–Birkhoff–Taylor ODE solver of order 10

A one-step 9-stage Hermite–Birkhoff–Taylor method of order 10, denoted by HBT(10)9, is constructed for solving nonstiff systems of first-order differential equations of the form y′=f(x,y), y(x ₀)=y ₀. The method uses y′ and higher derivatives y ⁽²⁾ to y ⁽⁴⁾ as in Taylor methods and is combined with a 9-stage Runge–Kutta method. Forcing a Taylor expansion of the numerical solution to agree with an expansion of the true solution leads to Taylor- and Runge–Kutta-type order conditions which are reorganized into Vandermonde-type linear systems whose solutions are the coefficients of the method. The new method has a larger scaled interval of absolute stability than Dormand–Prince DP(8,7)13M. The stepsize is controlled by means of y ⁽²⁾ and y ⁽⁴⁾. HBT(10)9 is superior to DP(8,7)13M and Taylor method of order 10 in solving several problems often used to test high-order ODE solvers on the basis of the number of steps, CPU time, and maximum global error. These numerical results show the benefits of adding high-order derivatives to Runge–Kutta methods.     

Pryce pre-analysis adapted to some DAE solvers

The ODE solver HBT(12)5 of order 12 (T. Nguyen-Ba, H. Hao, H. Yagoub, R. Vaillancourt, One-step 5-stage Hermite-Birkho-Taylor ODE solver of order 12, Appl. Math. Comput. 211 (2009) 313-328. doi:10.1016/j.amc.2009.01.043), which combines a Taylor series method of order 9 with a Runge-Kutta method of order 4, is expanded into the DAE solver HBT(12)5DAE of order 12. Dormand-Prince’s DP(8, 7)13M is also expanded into the DAE solver DP(8, 7)DAE. Pryce structural pre-analysis, extended ODEs and ODE first-order forms are adapted to these DAE solvers with a stepsize control based on local error estimators and a modified Pryce algorithm to advance integration. HBT(12)5DAE uses only the first nine derivatives of the unknown variables as opposed to the first 12 derivatives used by the Taylor series method T12DAE of order 12. Numerical results show the advantage of HBT(12)5DAE over T12DAE, DP(8, 7)DAE and other known DAE solvers.     

Oscillation criteria of a certain class of third order nonlinear delay differential equations with damping

In this paper, we are concerned with the oscillation of third order nonlinear delay differential equations of the form

(r₂(t)(r₁(t)y^′)′)′+p(t)y^′+q(t)f(y(g(t)))=0.     

A new class of efficient one-step contractivity preserving high-order time discretization methods of order 5 to 14

A family of one-step, explicit, contractivity preserving, multi-stage, multi-derivative, Hermite–Birkhoff–Taylor methods of order p =?5,6,…,14, that we denote by CPHBT_RK4(d,s,p), with nonnegative coefficients are constructed by casting s-stage Runge–Kutta methods of order 4 with Taylor methods of order d. The constructed CPHBT_RK4 methods are implemented using efficient variable step control and are compared to other well-known methods on a variety of initial value problems. A selected method: CP 6-stages 9-derivative HBT method of order 12, denoted by CPHBTRK412, has larger region of absolute stability than Dormand–Prince DP(8,7)13M and Taylor method T(12) of order 12. It is superior to DP(8,7)13M and T(12) methods on the basis the number of steps, CPU time, and maximum global error on several problems often used to test higher-order ODE solvers. Also, we show that the contractivity preserving property of CPHBTRK412is very efficient in suppressing the effect of the propagation of discretization errors and the new method compares positively with explicit 17 stages Runge-Kutta-Nyström pair of order 12 by Sharp et al. on a long-term integration of a standard N-body problem. The selected CPHBTRK412is listed in the Appendix.     

The behavior of solutions of second order delay differential equations

In this paper, we study the behavior of solutions of second order delay differential equation

y^″(t)=p₁y^′(t)+p₂y^′(t−τ)+q₁y(t)+q₂y(t−τ),     

An inverse problem for computing a leading coefficient in the Sturm-Liouville operator by using the boundary data

We consider an inverse problem for identifying a leading coefficient α(x) in −(α(x)y′(x))′ + q(x)y(x) = H(x), which is known as an inverse coefficient problem for the Sturm-Liouville operator. We transform y(x) to u(x, t) =  (1 + t)y(x) and derive a parabolic type PDE in a fictitious time domain of t. Then we develop a Lie-group adaptive method (LGAM) to find the coefficient function α(x). When α(x) is a continuous function of x, we can identify it very well, by giving boundary data of y, y′ and α. The efficiency of LGAM is confirmed by comparing the numerical results with exact solutions. Although the data used in the identification are limited, we can provide a rather accurate solution of α(x).     

Oscillation theorem for superlinear second order damped differential equations

In this paper, we are concerned with the oscillation of second order superlinear differential equations of the form

(a(t)y^′(t))^′+p(t)y^′(t)+q(t)f(y(t))=0.     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

Duality for a Class of Multiobjective Control Problems

In this paper, a class of multiobjective control problems is considered, where the objective and constraint functions involved are f(t, x(t), ?(t), y(t), z(t)) with x(t) ∈ Rⁿ, y(t) ∈ Rⁿ, and z(t) ∈ R^m, where x(t) and z(t) are the control variables and y(t) is the state variable. Under the assumption of invexity and its generalization, duality theorems are proved through a parametric approach to related properly efficient solutions of the primal and dual problems.     

One-step 9-stage Hermite–Birkhoff–Taylor DAE solver of order 10

The HBT(10)9 method for ODEs is expanded into HBT(10)9DAE for solving nonstiff and moderately stiff systems of fully implicit differential algebraic equations (DAEs) of arbitrarily high fixed index. A scheme to generate first-order derivatives at off-step points is combined with Pryce scheme which generates high order derivatives at step points. The stepsize is controlled by a local error estimator. HBT(10)9DAE uses only the first four derivatives of y instead of the first 10 required by Taylor’s series method T10DAE of order 10. Dormand–Prince’s DP(8,7)13M for ODEs is extended to DP(8,7)DAE for DAEs. HBT(10)9DAE wins over DP(8,7)DAE on several test problems on the basis of CPU time as a function of relative error at the end of the interval of integration. An index-5 problem is equally well solved by HBT(10)9DAE and T10DAE. On this problem, the error in the solution by DP(8,7)DAE increases as time increases.     

Fractal oscillations of self-adjoint and damped linear differential equations of second-order

For a prescribed real number s ∈ [1, 2), we give some sufficient conditions on the coefficients p(x) and q(x) such that every solution y = y(x), y ∈ C²((0, T]) of the linear differential equation (p(x)y′)′ + q(x)y = 0 on (0, T], is bounded and fractal oscillatory near x = 0 with the fractal dimension equal to s. This means that y oscillates near x = 0 and the fractal (box-counting) dimension of the graph Γ(y) of y is equal to s as well as the s dimensional upper Minkowski content (generalized length) of Γ(y) is finite and strictly positive. It verifies that y admits similar kind of the fractal geometric asymptotic behaviour near x = 0 like the chirp function y_ch(x) = a(x)S(φ(x)), which often occurs in the time-frequency analysis and its various applications. Furthermore, this kind of oscillations is established for the Bessel, chirp and other types of damped linear differential equations given in the form y″ + (μ/x)y′ + g(x)y = 0, x ∈ (0, T]. In order to prove the main results, we state a new criterion for fractal oscillations near x = 0 of real continuous functions which essentially improves related one presented in [1].     

New oscillation criteria of second-order nonlinear differential equations

By employing a class of new functions Φ=Φ(t,s,l) and a generalized Riccati technique, some new oscillation and interval oscillation criteria are established for the second-order nonlinear differential equation

(r(t)y^′(t))^′+Q(t,y(t),y^′(t))=0.     

Hermite–Birkhoff–Obrechkoff four-stage four-step ODE solver of order 14 with quantized step size

A four-stage Hermite–Birkhoff–Obrechkoff method of order 14 with four quantized variable steps, denoted by HBOQ(14)4, is constructed for solving non-stiff systems of first-order differential equations of the form y^′=f(t,y)$y^{\prime }=f(t,y)$ with initial conditions y(_t0)=_y0$y\left(t_0\right)=y_0$. Its formula uses y^′$y^{\prime }$, y^″$y^{″}$ and y^?$y^{?}$ as in Obrechkoff methods. Forcing a Taylor expansion of the numerical solution to agree with an expansion of the true solution leads to multistep- and Runge–Kutta-type order conditions which are reorganized into linear Vandermonde-type systems. To reduce overhead, simple formulae are derived only once to obtain the values of Hermite–Birkhoff interpolation polynomials in terms of Lagrange basis functions for 16 quantized step size ratios. The step size is controlled by a local error estimator. When programmed in C ++, HBOQ(14)4 is superior to the Dormand–Prince Runge–Kutta pair DP(8,7)13M of order 8 in solving several problems often used to test higher order ODE solvers at stringent tolerances. When programmed in Matlab, it is superior to ode113 in solving costly problems, on the basis of the number of steps, CPU time, and maximum global error. The code is available on the URL www.site.uottawa.ca/~remi.     

Existence of a positive solution to a system of discrete fractional boundary value problems

We analyze a system of discrete fractional difference equations subject to nonlocal boundary conditions. We consider the system of equations given by -Δ^ν_iy_i(t)=λ_ia_i(t+ν_i-1)f_i(y₁(t+ν₁-1),y₂(t+ν₂-1)), for t∈[0,b]_N0, subject to y_i(ν_i − 2) = ψ_i(y_i) and y_i(ν_i + b) = ?_i(y_i), for i = 1, 2, where ψ_i,?_i:R^b+3→R are given functionals. We also assume that ν_i ∈ (1, 2], for each i. Although we assume that both a_i and f_i(y₁, y₂) are nonnegative for each i, we do not necessarily presume that each ψ_i(y_i) and ?_i(y_i) is nonnegative for each i and each y_i ? 0. This generalizes some recent results both on discrete fractional boundary value problems and on discrete integer-order boundary value problems, and our techniques provide new results in each case.     

Meromorphic Functions That Share Fixed-Points

In this paper, we prove the following result: Let f(z) and g(z) be two nonconstant meromorphic(entire) functions, n ≥ 11(n ≥ 6) a positive integer. If fⁿ(z)f′(z) and gⁿ(z)g′(z) have the same fixed-points, then either f(z) = c₁e^cz2, g(z) = c₂e^− cz2, where c₁, c₂, and c are three constants satisfying 4(c₁c₂)^n + 1c² = −1, or f(z) ≡ tg(z) for a constant t such that t^n + 1 = 1.     

Stability of continuous Runge–Kutta-type methods for nonlinear neutral delay-differential equations

This paper provides results on the correct simulation, when using continuous Runge–Kutta methods, of certain stability properties of nonlinear neutral delay-differential equations (NDDEs) y^′(t)=f(t,y(t),y(t-τ(t)),y^′(t-τ(t)))(t?_t0)$y^{\prime }(t)=f(t\text{,}y(t)\text{,}y(t-\tau (t))\text{,}{y}^{\prime }(t-\tau (t)))(t?t_0)$. In particular, it is shown that certain continuous Runge–Kutta methods based upon the backward Euler method or the 2-stage Lobatto IIIC method, combined with linear interpolation, are GRN$\mathit{GRN}$-stable and asymptotically stable for NDDEs.     

On the existence of periodic solutions to p-Laplacian Rayleigh differential equation with a delay

In this paper, the authors study the existence of periodic solutions to a p-Laplacian Rayleigh differential equation with a delay as follows:

(φp(y^′(t))^)′+f(y^′(t))+g(y(t−τ(t)))=e(t),     

A variational approach to dislocation problems for periodic Schrödinger operators

As a simple model for lattice defects like grain boundaries in solid state physics we consider potentials which are obtained from a periodic potential V=V(x,y) on R² with period lattice Z² by setting Wt(x,y)=V(x+t,y) for x<0 and Wt(x,y)=V(x,y) for x?0, for t∈[0,1]. For Lipschitz-continuous V it is shown that the Schrödinger operators Ht=−Δ+Wt have spectrum (surface states) in the spectral gaps of H₀, for suitable t∈(0,1). We also discuss the density of these surface states as compared to the density of the bulk. Our approach is variational and it is first applied to the well-known dislocation problem (Korotyaev (2000, 2005) [15] and [16]) on the real line. We then proceed to the dislocation problem for an infinite strip and for the plane. In Appendix A, we discuss regularity properties of the eigenvalue branches in the one-dimensional dislocation problem for suitable classes of potentials.     

Two-point Taylor expansions in the asymptotic approximation of double integrals. Application to the second and fourth Appell functions

The main difficulty in Laplace's method of asymptotic expansions of double integrals is originated by a change of variables. We consider a double integral representation of the second Appell function F₂(a,b,b^′,c,c^′;x,y) and illustrate, over this example, a variant of Laplace's method which avoids that change of variables and simplifies the computations. Essentially, the method only requires a Taylor expansion of the integrand at the critical point of the phase function. We obtain in this way an asymptotic expansion of F₂(a,b,b^′,c,c^′;x,y) for large b, b^′, c and c^′. We also consider a double integral representation of the fourth Appell function F₄(a,b,c,d;x,y). We show, in this example, that this variant of Laplace's method is uniform when two or more critical points coalesce or a critical point approaches the boundary of the integration domain. We obtain in this way an asymptotic approximation of F₄(a,b,c,d;x,y) for large values of a,b,c and d. In this second example, the method requires a Taylor expansion of the integrand at two points simultaneously. For this purpose, we also investigate in this paper Taylor expansions of two-variable analytic functions with respect to two points, giving Cauchy-type formulas for the coefficients of the expansion and details about the regions of convergence.