One-step 5-stage Hermite-Birkhoff-Taylor ODE solver of order 12

A one-step 5-stage Hermite-Birkhoff-Taylor method, HBT(12)5, of order 12 is constructed for solving nonstiff systems of differential equations y^′=f(t,y), y(t₀)=y₀, where y∈Rⁿ. The method uses derivatives y^′ to y⁽⁹⁾ as in Taylor methods combined with a 5-stage Runge-Kutta method. Forcing an expansion of the numerical solution to agree with a Taylor expansion of the true solution to order 12 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. HBT(12)5 has a larger interval of absolute stability than Dormand-Prince DP(8, 7)13M and Taylor method T12 of order 12. The new method has also a smaller norm of principal error term than T12. It is superior to DP(8, 7)13M and T12 on the basis the number of steps, CPU time and maximum global error on common test problems. The formulae of HBT(12)5 are listed in an appendix.     

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.     

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.     

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.     

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.     

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.     

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−τ),     

A three-stage, VSVO, Hermite-Birkhoff-Taylor, ODE solver

The new variable-step, variable-order, ODE solver, HBT(p) of order p, presented in this paper, combines a three-stage Runge-Kutta method of order 3 with a Taylor series method of order p-2 to solve initial value problems , where y:R→R^d and f:R×R^d→R^d. The order conditions satisfied by HBT(p) are formulated and they lead to Vandermonde-type linear algebraic systems whose solutions are the coefficients in the formulae for HBT(p). A detailed formulation of variable-step HBT(p) in both fixed-order and variable-order modes is presented. The new method and the Taylor series method have similar regions of absolute stability. To obtain high-accuracy results at high order, this method has been implemented in multiple precision.     

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.     

Regions of absolute stability of explicit Runge-Kutta-Nyström methods for y″ = f(x, y, y′)

We examine regions of absolute stability of s-stage explicit Runge-Kutta-Nyström (R-K-N) methods of order s for s = 2, 3, 4 for y″ = f(x, y, y′) by applying these methods to the test equation: y″ + 2αy′ + β²y = O, α, β ? 0, α + β > 0. We also examine the regions of absolute stability of Runge-Kutta (R-K) methods for first order differential equations of respective orders. Interestingly, it turns out that regions of absolute stability of R-K methods and R-K-N methods of respective orders for which the asymptotic relative error does not deteriorate are identical. Our present investigations are in continuation of the recent results of Chawla and Sharma [1].     

On the Attainable Order of Collocation Methods for Delay Differential Equations with Proportional Delay

To analyze the attainable order of m-stage implicit (collocation-based) Runge-Kutta methods for the delay differential equation (DDE) y′(t) = by(qt), 0 < q ≤ 1 with y(0) = 1, and the delay Volterra integral equation (DVIE) y(t) = 1 + $\tfrac{b}{q}\int {_0^{qt} }$ y(s) ds with proportional delay qt, 0 < q ≤ 1, our particular interest lies in the approximations (and their orders) at the first mesh point t = h for the collocation solution v(t) of the DDE and the iterated collocation solution u _it(t) of the DVIE to the solution y(t). Recently, H. Brunner proposed the following open problem: “For m ≤ 3, do there exist collocation points c _i = c _i(q), i = 1, 2,..., m in [0,1] such that the rational approximant v(h)is the (m, m)-Padé approximant to y(h)? If these exist, then |v(h) ? y(h)| = O(h ^2m+1) but what is the collocation polynomial M _m(t; q) = K Π _i=1 ^m (t ? c _i) of v(th), t ∈ [0, 1]?” In this paper, we solve this question affirmatively, and give the related results between the collocation solution v(t) of the DDE and the iterated collocation solution u _it(t) of the DVIE. We also answer to Brunner's second open question in the case that one collocation point is fixed at the right end point of the interval.     

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.     

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),     

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.     

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.     

Absolute stability of explicit Runge-Kutta-Nyström methods for y″=ƒ(x,y, y′)

We examine absolute stability of s-stage explicit Runge-Kutta-Nyström (R-K-N) methods of order s for $\text{s=2, 3, 4}\text{for}\text{y″=?(x, y, y′)}$ by applying these methods to the test equation: y″+2λy′+λ²y=0, λ>0. We show the existence of R-K-N methods of orders two, three and four possessing intervals of absolute stability as large as that of explicit Runge-Kutta (r-K) methods of respective orders.     

Stabilization of unbounded bilinear control systems in Hilbert space

We consider bilinear control systems of the form y^′(t)=Ay(t)+u(t)By(t) where A generates a strongly continuous semigroup of contraction (etA)_t?0 on an infinite-dimensional Hilbert space Y whose scalar product is denoted by 〈.,.〉. We suppose that this system is unbounded in the sense that the linear operator B is unbounded from the state Y into itself. Tacking into account eventual control saturation, we study the problem of stabilization by (possibly nonquadratic) feedback of the form u(t)=−f(〈By(t),y(t)〉). Applications to the heat equation is considered.     

Minimum norm solutions of single stiff linear analytic differential equations

The $\text{l}$₂-norm of the infinite vector of the terms of the Taylor series of an analytic function is used to measure the “unsmoothness” of the function. The sets of solutions to the scalar differential equations y′(t) = λy(t) + f(t) and y′(t) = q(t)y(t) + f(t) are analyzed with respect to this norm. A number of results on the particular solution with minimum norm are given.     

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.     

A convexity-preservingC 2 parametric rational cubic interpolation

Summary AC ² parametric rational cubic interpolantr(t)=x(t) i+y(t) j,t[t ₁,t _n] to data S={(x_j, y_j)|j=1,...,n} is defined in terms of non-negative tension parameters _j ,j=1,...,n–1. LetP be the polygonal line defined by the directed line segments joining the points (x _j ,y _j ),t=1,...,n. Sufficient conditions are derived which ensure thatr(t) is a strictly convex function on strictly left/right winding polygonal line segmentsP. It is then proved that there always exist _j ,j=1,...,n–1 for whichr(t) preserves the local left/righ winding properties of any polygonal lineP. An example application is discussed.This research was supported in part by the natural Sciences and Engineering Research Council of Canada.