Quadratically Convergent Method for Simultaneously Approaching the Roots of Polynomial Solutions of a Class of Differential Equations: Application to Orthogonal Polynomials

This paper investigates the application of the method introduced by L. Pasquini (1989) for simultaneously approaching the zeros of polynomial solutions to a class of second-order linear homogeneous ordinary differential equations with polynomial coefficients to a particular case in which these polynomial solutions have zeros symmetrically arranged with respect to the origin. The method is based on a family of nonlinear equations which is associated with a given class of differential equations. The roots of the nonlinear equations are related to the roots of the polynomial solutions of differential equations considered. Newton's method is applied to find the roots of these nonlinear equations. In (Pasquini, 1994) the nonsingularity of the roots of these nonlinear equations is studied. In this paper, following the lines in (Pasquini, 1994), the nonsingularity of the roots of these nonlinear equations is studied. More favourable results than the ones in (Pasquini, 1994) are proven in the particular case of polynomial solutions with symmetrical zeros. The method is applied to approximate the roots of Hermite–Sobolev type polynomials and Freud polynomials. A lower bound for the smallest positive root of Hermite–Sobolev type polynomials is given via the nonlinear equation. The quadratic convergence of the method is proven. A comparison with a classical method that uses the Jacobi matrices is carried out. We show that the algorithm derived by the proposed method is sometimes preferable to the classical QR type algorithms for computing the eigenvalues of the Jacobi matrices even if these matrices are real and symmetric.     

Bifurcations of steady states for the Eguchi–Oki–Matsumura model of phase separation

This article is concerned with bifurcations of steady states for a model system of phase separation, which is introduced by Eguchi–Oki–Matsumura (EOM). The system consists of coupled two evolution equations and admits steady state solutions with different energies. The bifurcation phenomena of these steady states with respect to the principal parameter, which is related to the temperature, are analyzed.     

Applications of Sturm sequences to bifurcation analysis of delay differential equation models

This paper formalizes a method used by several others in the analysis of biological models involving delay differential equations. In such a model, the characteristic equation about a steady state is transcendental. This paper shows that the analysis of the bifurcation due to the introduction of the delay term can be reduced to finding whether a related polynomial equation has simple positive real roots. After this result has been established, we utilize Sturm sequences to determine whether a polynomial equation has positive real roots. This work has extended the stability results found in previous papers and provides a novel theorem about stability switches for low degree characteristic equations.     

Khovanskii–Rolle Continuation for Real Solutions

We present a new continuation algorithm to find all real solutions to a nondegenerate system of polynomial equations. Unlike homotopy methods, the algorithm is not based on a deformation of the system; instead, it traces real curves connecting the solutions to one system of equations to those of another, eventually leading to the desired real solutions. It also differs from homotopy methods in that it follows only real paths and computes no complex solutions to the original equations. The number of curves traced is essentially bounded above by the fewnomial bound for real solutions, and the method takes advantage of any slack in that bound.     

Conditions for polynomial Liénard centers

In 1999, Christopher gave a necessary and sufficient condition for polynomial Li′enard centers, which requires coupled functional equations, where the primitive functions of the damping function and the restoring function are involved, to have polynomial solutions. In order to judge whether the coupled functional equations are solvable, in this paper we give an algorithm to compute a Gr¨obner basis for irreducible decomposition of algebraic varieties so as to find algebraic relations among coefficients of the damping function and the restoring function. We demonstrate the algorithm for polynomial Li′enard systems of degree 5, which are divided into 25 cases. We find all conditions of those coefficients for the polynomial Li′enard center in 13 cases and prove that the origin is not a center in the other 12 cases.     

Solving triangular algebraic systems by means of simultaneous iterations

Our aim in this paper is to extend a variant of the Weierstrass method for the simultaneous computation of the solutions of a triangular algebraic system of equations. The appropriate tools are the symmetric functions of the roots of a polynomial. Using these symmetric functions we give another equivalent formulation for the search of all the roots of a triangular algebraic system. Using the latter formulation our method consists in solving a more simple system (where partial degrees of all the equations do not exceed 1) by Newton’s method. The quadratic convergence of our method is an immediate consequence of Newton’s method and need not be proved explicitly. This revised version was published online in June 2006 with corrections to the Cover Date.     

Steady state bifurcation of a periodically excited system under delayed feedback controls

This paper investigates the steady state bifurcation of a periodically excited system subject to time-delayed feedback controls by the combined method of residue harmonic balance and polynomial homotopy continuation. Three kinds of delayed feedback controls are considered to examine the effects of different delayed feedback controls and delay time on the steady state response. By means of polynomial homotopy continuation, all the possible steady state solutions corresponding the third-order superharmonic and second-subharmonic responses are derived analytically, i.e. without numerical integration. It is found that the delayed feedback changes the bifurcating curves qualitatively and possibly eliminates the saddle-node bifurcation during resonant. The delayed position-velocity coupling and the delayed velocity feedback controls can destabilize the steady state responses. Coexisting periodic solutions, period-doubling bifurcation and even chaos are found in these control systems. The neighborhood of the periodic solutions is verified numerically in the phase portraits. The various effects of time delay on the steady state response are investigated. Many new phenomena are observed.     

稳定性判定与多项式求根算法

本文给出了一种判定多项式根是否全在单位圆内的简便方法．该方法可用于判定离散控制系统的稳定性和求多项式的全部根。     

Qualitative analysis for a ratio-dependent predator-prey model with disease and diffusion

This paper is purported to study a reaction diffusion system arising from a ratio-dependent predator-prey model with disease. We study the dynamical behavior of the predator-prey system. The conditions for the permanent and existence of steady states and their stability are established. We can obtain the bounds for positive steady state of the corresponding elliptic system. The non-existence results of non-constant positive solutions are derived.     

Solution to nonlinear MHDS arising from optimal growth problems

In this paper we propose a method for solving in closed form a general class of nonlinear modified Hamiltonian dynamic systems (MHDS). This method is used to analyze the intertemporal optimization problem from endogenous growth theory, especially the cases with two controls and one state variable. We use the exact solutions to study both uniqueness and indeterminacy of the optimal path when the dynamic system has not a well-defined isolated steady state. With this approach we avoid the linearization process, as well as the reduction of dimension technique usually applied when the dynamic system offers a continuum of steady states or no steady state at all.     

The effects of coupling on pattern formation in a simple autocatalytic system

The spatiotetnporal structures that can arise in two identicalcells, each governed by cubic autocatalator kinetics and coupledvia the diffusive interchange of a reactant, are discussed.The coupling gives rise to five spatially uniform steady states,one of which exists in the uncoupled system. By studying thelinearized equations, it is found that three of these steadystates, including that of the uncoupled system, may give riseto the possibility of bifurcations to spatially nonuniform steadystates. In the case of the steady state corresponding to thatof the uncoupled system, it is seen that the coupling leadsto bifurcations not present in the uncoupled system which giverise to locally stable nonuniform steady states. A weakly nonlinearanalysis is developed for both small and large coupling strengtha, and for parameter values in a neighbourhood of the bifurcationpoints on the new steady states. This clarifies the nature ofthe nonuniform solutions close to bifurcation, which are thenfollowed numerically using a path-following technique. The couplingis found to produce extra nonuniform steady solutions whichare stable close to their bifurcation points.     

Positive solutions of some parabolic system with cross-diffusion and nonlocal initial conditions

The paper is concerned with a system consisting of two coupled nonlinear parabolic equations with a cross-diffusion term, where the solutions at positive times define the initial states. The equations arise as steady state equations of an age-structured predator–prey system with spatial dispersion. Based on unilateral global bifurcation methods for Fredholm operators and on maximal regularity for parabolic equations, global bifurcation of positive solutions is derived.     

Existence of periodic travelling waves solutions in predator prey model with diffusion

This paper deals with the qualitative analysis of the travelling waves solutions of a reaction diffusion model that refers to the competition between the predator and prey with modified Leslie–Gower and Holling type II schemes. The well posedeness of the problem is proved. We establish sufficient conditions for the asymptotic stability of the unique nontrivial positive steady state of the model by analyzing roots of the forth degree exponential polynomial characteristic equation. We also prove the existence of a Hopf bifurcation which leads to periodic oscillating travelling waves by considering the diffusion coefficient as a parameter of bifurcation. Numerical simulations are given to illustrate the analytical study.     

A Gröbner Free Alternative for Polynomial System Solving

Given a system of polynomial equations and inequations with coefficients in the field of rational numbers, we show how to compute a geometric resolution of the set of common roots of the system over the field of complex numbers. A geometric resolution consists of a primitive element of the algebraic extension defined by the set of roots, its minimal polynomial, and the parametrizations of the coordinates. Such a representation of the solutions has a long history which goes back to Leopold Kronecker and has been revisited many times in computer algebra. We introduce a new generation of probabilistic algorithms where all the computations use only univariate or bivariate polynomials. We give a new codification of the set of solutions of a positive dimensional algebraic variety relying on a new global version of Newton's iterator. Roughly speaking the complexity of our algorithm is polynomial in some kind of degree of the system, in its height, and linear in the complexity of evaluation of the system. We present our implementation in the Magma system which is called Kronecker in homage to his method for solving systems of polynomial equations. We show that the theoretical complexity of our algorithm is well reflected in practice and we exhibit some cases for which our program is more efficient than the other available software.     

Global stability and travelling waves of a predator-prey model with diffusion and nonlocal maturation delay

In this paper, a diffusive predator-prey system with nonlocal maturation delay is investigated. By analyzing the corresponding characteristic equations, the local stability of each of uniform steady states of the system is discussed. Sufficient conditions are derived for the global stability of the positive steady state and the semi-trivial steady state of the system by using the method of upper–lower solutions and its associated monotone iteration scheme, respectively. The existence of travelling wave solution of the system is established by using the geometric singular perturbation theory. Numerical simulations are carried out to illustrate the theoretical results.     

Global boundedness and stability of solutions for prey-taxis model with handling and searching predators

In this paper, we show the global boundedness and stability of solutions for prey-taxis model with handling and searching predators in a two-dimensional bounded domain with smooth boundary. First, entropy-like equations and boundedness criteria are derived, and it is proved that the system has a unique uniformly bounded global classical solution. In addition, we show that prey-only steady state is globally asymptotically stable if the predator is weak. The convergence rate of solutions to the steady states is derived in the paper.     

The impact of finite precision arithmetic and sensitivity on the numerical solution of partial differential equations

In this paper, we address some fundamental issues concerning “time marching” numerical schemes for computing steady state solutions of boundary value problems for nonlinear partial differential equations. Simple examples are used to illustrate that even theoretically convergent schemes can produce numerical steady state solutions that do not correspond to steady state solutions of the boundary value problem. This phenomenon must be considered in any computational study of nonunique solutions to partial differential equations that govern physical systems such as fluid flows. In particular, numerical calculations have been used to “suggest” that certain Euler equations do not have a unique solution. For Burgers' equation on a finite spatial interval with Neumann boundary conditions the only steady state solutions are constant (in space) functions. Moreover, according to recent theoretical results, for any initial condition the corresponding solution to Burgers' equation must converge to a constant as t → ∞. However, we present a convergent finite difference scheme that produces false nonconstant numerical steady state “solutions.” These erroneous solutions arise out of the necessary finite floating point arithmetic inherent in every digital computer. We suggest the resulting numerical steady state solution may be viewed as a solution to a “nearby” boundary value problem with high sensitivity to changes in the boundary conditions. Finally, we close with some comments on the relevance of this paper to some recent “numerical based proofs” of the existence of nonunique solutions to Euler equations and to aerodynamic design.     

Spatial resonance and Turing–Hopf bifurcations in the Gierer–Meinhardt model

Gierer–Meinhardt system as a molecularly plausible model has been proposed to formalize the observation for pattern formation. In this paper, the Gierer–Meinhardt model without the saturating term is considered. By the linear stability analysis, we not only give out the conditions ensuring the stability and Turing instability of the positive equilibrium but also find the parameter values where possible Turing–Hopf and spatial resonance bifurcation can occur. Then we develop the general algorithm for the calculations of normal form associated with codimension-2 spatial resonance bifurcation to better understand the dynamics neighboring of the bifurcating point. The spatial resonance bifurcation reveals the interaction of two steady state solutions with different modes. Numerical simulations are employed to illustrate the theoretical results for both the Turing–Hopf bifurcation and spatial resonance bifurcation. Some expected solutions including stable spatially inhomogeneous periodic solutions and coexisting stable spatially steady state solutions evolve from Turing–Hopf bifurcation and spatial resonance bifurcation respectively.     

Travelling waves of a predator–prey model with a spatiotemporal delay

This paper is concerned with the existence of traveling waves to a predator–prey model with a spatiotemporal delay. By analyzing the corresponding characteristic equations, the local stability of a positive steady state and each of boundary steady states are established, and the existence of Hopf bifurcation at the positive steady state is also discussed. By constructing a pair of upper–lower solutions and by using the cross‐iteration method as well as the Schauder's fixed‐point theorem, the existence of a traveling wave solution connecting the semi‐trivial steady state and the positive steady state is proved. Numerical simulations are carried out to illustrate the main theoretical results. Copyright © 2013 John Wiley & Sons, Ltd.     

Monotone methods and stability results for nonlocal reaction-diffusion equations with time delay

In this paper, we study the applications of the monotone iteration method for investigating the existence and stability of solutions to nonlocal reaction-diffusion equations with time delay. We emphasize the importance of the idea of monotone iteration schemes for investigating the stability of solutions to such equations. We show that every steady state of such equations obtained by using the monotone iteration method is priori stable and all stable steady states can be obtained by using such method. Finally, we apply our main results to three population models.