Oscillation of Numerical Solution in the Runge-Kutta Methods for Equation x＇（t） = ax（t） ＋ aox（[t]）

The paper de ls with oscillation of Runge-Kutta methods for equation x＇（t） = ax（t） ＋ aox（[t]）. The conditions of oscillation for the numerical methods are presented by considering the characteristic equation of the corresponding discrete scheme. It is proved that any nodes have the same oscillatory property as the integer nodes. Furthermore, the conditions under which the oscillation of the analytic solution is inherited by the numerical solution are obtained. The relationships between stability and oscillation are considered. Finally, some numerical experiments are given.     

Analysis of a Prey-predator Model with Disease in Prey

In this paper, a system of reaction-diffusion equations arising in ecoepidemiological systems is investigated. The equations model a situation in which a predator species and a prey species inhabit the same bounded region and the predator only eats the prey with transmissible diseases. Local stability of the constant positive solution is considered. A number of existence and non-existence results about the nonconstant steady states of a reaction diffusion system are given. It is proved that if the diffusion coefficient of the prey with disease is treated as a bifurcation parameter, non-constant positive steady-state solutions may bifurcate from the constant steadystate solution under some conditions.     

Two-dimensional wave equations with fractal boundaries

This paper focuses on two cases of two-dimensional wave equations with fractal boundaries. The first case is the equation with classical derivative. The formal solution is obtained. And a definition of the solution is given. Then we prove that under certain conditions, the solution is a kind of fractal function, which is continuous, differentiable nowhere in its domain. Next, for specific given initial position and 3 different initial velocities, the graphs of solutions are sketched. By computing the box dimensions of boundaries of cross-sections for solution surfaces, we evaluate the range of box dimension of the vibrating membrane. The second case is the equation with p-type derivative. The corresponding solution is shown and numerical example is given.     

STABILITY OF GENERAL LINEAR METHODS FOR SYSTEMS OF FUNCTIONAL-DIFFERENTIAL AND FUNCTIONAL EQUATIONS

This paper is concerned with the numerical solution of functional-differential and func-tional equations which include functional-differential equations of neutral type as special cases. The adaptation of general linear methods is considered. It is proved that A-stable general linear methods can inherit the asymptotic stability of underlying linear systems.Some general results of numerical stability are also given.     

带分数阶Robin边界条件的时间-空间分数阶扩散方程的有限差分方法（英文）

In this paper, an efficient numerical method is proposed to solve the Caputo-Riesz fractional diffusion equation with fractional Robin boundary conditions. We approximate the Riesz space fractional derivatives using the fractional central difference scheme with second-order accurate. A priori estimation of the solution of the numerical scheme is given, and the stability and convergence of the numerical scheme are analyzed.Finally, a numerical example is used to verify the accuracy and efficiency...     

THE STABILITY OF θ-METHODS FOR PANTOGRAPH DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for pantograph differential equations with variable delays. We investigate the stability of one leg θ-methods in the numerical solution of these problems. Sufficient conditions for the asymptotic stability of θ-methods are given by Fourier analysis and Ergodic theory.     

具有时滞和毒素系统的动力学分析

In this paper, a toxin producing phytoplankton-zooplankton model with inhibitory substrate and time delay is investigated. A discrete time delay is induced to both of the consume response function and distribution of toxic substance term. Moreover, Tissiet type function is used for zooplankton grazing to account for the effect of toxication by the TPP population. The conditions to guarantee the coexistence of two species and stability of coexistence equilibrium are given. In particular, we show that there exist critical values of the delay parameters below which the coexistence equilibrium is stable and above which it is unstable. Hopf bifurcation occurs when the delay parameters cross their critical values. Some numerical simulations are executed to validate the analytical findings.     

A projected gradient method with nonmonotonic backtracking technique for solving convex constrained monotone variational inequality problem

Based on a differentiable merit function proposed by Taji, et al in “Mathematical Programming, 1993, 58： 369-383”, a projected gradient trust region method for the monotone variational inequality problem with convex constraints is presented. Theoretical analysis is given which proves that the proposed algorithm is globally convergent and has a local quadratic convergence rate under some reasonable conditions. The results of numerical experiments are reported to show the effectiveness of the proposed algorithm.     

STABILITY OF THEORETICAL SOLUTION AND NUMERICAL SOLUTION OF NONLINEAR DIFFERENTIAL EQUATIONS WITH PIECEWISE DELAYS

This paper is concerned with the stability of theoretical solution and numerical solutionof a class of nonlinear differential equations with piecewise delays.At first,a sufficientcondition for the stability of theoretical solution of these problems is given,then numericalstability and asymptotical stability are discussed for a class of multistep methods whenapplied to these problems.     

EXPONENTIAL STABILITY FOR NONLINEAR HYBRID STOCHASTIC PANTOGRAPH EQUATIONS AND NUMERICAL APPROXIMATION

The paper develops exponential stability of the analytic solution and convergence in probability of the numerical method for highly nonlinear hybrid stochastic pantograph equation. The classical linear growth condition is replaced by polynomial growth conditions, under which there exists a unique global solution and the solution is almost surely exponentially stable. On the basis of a series of lemmas, the paper establishes a new criterion on convergence in probability of the Euler-Maruyama approximate solution. The criterion is very general so that many highly nonlinear stochastic pantograph equations can obey these conditions. A highly nonlinear example is provided to illustrate the main theory.     

Numerical stability analysis of differential equations with piecewise constant arguments with complex coefficients

In this paper we consider the analytical and numerical stability regions of Runge-Kutta methods for differential equations with piecewise continuous arguments with complex coefficients. It is shown that the analytical stability region contained in the numerical one is violated for a∉R by the geometric technique. And we give the conditions under which the analytical stability region is contained in the union of the numerical stability regions of two Runge-Kutta methods. At last, some experiments are given.     

Numerical Stability and Oscillations of Runge-Kutta Methods for Differential Equations with Piecewise Constant Arguments of Advanced Type

For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the numerical stability region contains the analytic stability region are given. The conditions of oscillations for the Runge-Kutta methods are obtained also. We prove that the Runge-Kutta methods preserve the oscillations of the analytic solution. Moreover, the relationship between stability and oscillations is discussed. Several numerical examples which confirm the results of our analysis are presented.     

Runge–Kutta Methods for Systems of Differential Equation with Piecewise Continuous Arguments: Convergence and Stability

This work deals with the convergence and stability of Runge–Kutta methods for systems of differential equation with piecewise continuous arguments x^′(t) = Px(t)+Qx([t+1∕2]) under two cases for coe?cient matrix. First, when P and Q are complex matrices, the su?cient condition under which the analytic solution is asymptotically stable is given. It is proven that the Runge–Kutta methods are convergent with order p. Moreover, the su?cient condition under which the analytical stability region is contained in the numerical stability region is obtained. Second, when P and Q are commutable Hermitian matrices, using the theory of characteristic, the necessary and su?cient conditions under which the analytic solution and the numerical solution are asymptotically stable are presented, respectively. Furthermore, whether the Runge–Kutta methods preserve the stability of analytic solution are investigated by the theory of Padé approximation and order star. To demonstrate the theoretical results, some numerical experiments are adopted.     

向前型分段连续微分方程的数值解

本文讨论了向前型分段连续微分方程Euler-Maclaurin方法的收敛性和稳定性,给出了Euler-Maclaurin方法的稳定条件,证明了方法的收敛阶是2n+2,并且得到了数值解稳定区域包含解析解稳定区域的条件,最后给出了一些数值例子用以验证本文结论的正确性.     

The Numerical Solution of the Fourth Boundary Value Problem for Parabolic Partial Differential Equations

A direct method for the numerical solution of the implicit finitedifference equations derived from a parabolic differential equationwith periodic spatial boundary conditions is presented in algorithmicfrom. Consideration is given to the stability of the roundingerrors involved in the solution process and numerical resultsare derived which compare favourably with those obtained fromthe analytical solution and a matrix spectral resolution methodwhich is closely allied to the method of lines.     

Stability and oscillations of numerical solutions for differential equations with piecewise continuous arguments of alternately advanced and retarded type

This paper is concerned with the numerical properties of θ-methods for the solution of alternately advanced and retarded differential equations with piecewise continuous arguments. Using two θ-methods, namely the one-leg θ-method and the linear θ-method, the necessary and sufficient conditions under which the analytic stability region is contained in the numerical stability region are obtained, and the conditions of oscillations for the θ-methods are also obtained. It is proved that oscillations of the analytic solution are preserved by the θ-methods. Furthermore, the relationships between stability and oscillations are revealed. Some numerical experiments are presented to illustrate our results.     

Stability of Runge-Kutta-Pouzet methods for Volterra integro-differential equations with delays

This paper is concerned with the study of the stability of Runge-Kutta-Pouzet methods for Volterra integro-differential equations with delays. We are interested in the comparison between the analytical and numerical stability regions. First, we focus on scalar equations with real coefficients. It is proved that all Gauss-Pouzet methods can retain the asymptotic stability of the analytical solution. Then, we consider the multidimensional case. A new stability condition for the stability of the analytical solution is given. Under this condition, the asymptotic stability of Gauss-Pouzet methods is investigated.      

一类鞍点型线性半连续动力系统的周期解

讨论了一类鞍点型线性半连续动力系统,给出了阶1周期解存在性、唯一性、稳定性以及阶2周期解不存在性的条件,并指出在一定的条件下,系统可以存在无穷多个阶1周期解.最后给出了相应理论结果的数值模拟.     

滞后型分段连续随机微分方程的稳定性

本文研究了滞后型分段连续随机微分方程的解析稳定性和数值稳定性问题.首先,利用伊藤公式等方法获得了解析解均方稳定的条件,其次,对于包括均方稳定和T-稳定在内的Euler-Maruyama方法的数值稳定性问题,运用不等式技术和随机分析方法获得了一些新的结果,证明了在一定条件下,Euler-Maruyama方法既是均方稳定又是T-稳定的,推广了随机延迟微分方程的数值稳定性结论.     

Backward Euler-Maruyama method applied to nonlinear hybrid stochastic differential equations with time-variable delay

In this paper, we consider strong convergence and almost sure exponential stability of the backward Euler-Maruyama method for nonlinear hybrid stochastic differential equations with time-variable delay. Under the local Lipschitz condition and polynomial growth condition, it is proved that the backward Euler-Maruyama method is strongly convergent. Additionally, the moment estimates and almost sure exponential stability for the analytical solution are proved. Also, under the appropriate condition, we show that the numerical solutions for the backward Euler-Maruyama methods are almost surely exponentially stable. A numerical experiment is given to illustrate the computational effectiveness and the theoretical results of the method.