Marchuk’s Models of Infection Diseases: New Developments

It is well known that a model of any electrical circuit can be translated into a system of ordinary differential equations. Usually, the simulation is performed in the time-domain using classical discrete approaches convenient for both linear and non-linear systems. An other method consists in the quantization on the variable magnitude which is based on discrete events. The discrete event system specification formalism is adapted to simulate any electrical circuit described with a system of ordinary differential equations using the system state quantification. In this paper, a solution to make more efficient the magnitude quantification for a system of linear ordinary differential equations is given. The proposed technique is used to simulate a simple three phase coupled power systems and it is compared to other solutions given by the explicit integration using formal language and the time-domain discretization using fixed-step fifth order Runge-Kutta integration method.

We consider boundary value problems with periodic boundary conditions for first-order linear systems of impulsive ordinary differential equations with unknown right-hand sides and jumps of solutions at the impulse points entering into the statement of these problems which are assumed to be subjected to some quadratic restrictions. From indirect noisy observations of their solutions on a finite system of intervals, we obtain the optimal, in certain sense, estimates of images of their right-hand sides under linear continuous operators. Under the condition that the unknown correlation functions of noises in observations belong to some special sets, it is established that such estimates and estimation errors are expressed explicitly via solutions of special periodic boundary value problems for linear impulsive systems of ordinary differential equations.

The non-standard Hamiltonian system, also referred to as a partial Hamiltonian system in the literature, of the form q ˙ i = ∂ H ∂ p i , p ˙ i = − ∂ H ∂ q i + Γ i ( t , q i , p i ) ${\dot q^i} = \frac{{\partial H}}{{\partial {p_i}}},{\text{ }}{\dot p^i} = - \frac{{\partial H}}{{\partial {q_i}}} + {\Gamma ^i}(t,{\text{ }}{q^i},{\text{ }}{p_i})$ appears widely in economics, physics, mechanics, and other fields. The non-standard (partial) Hamiltonian systems arise from physical Hamiltonian structures as well as from artificial Hamiltonian structures. We introduce the term ‘artificial Hamiltonian’ for the Hamiltonian of a model having no physical structure. We provide here explicitly the notion of an artificial Hamiltonian for dynamical systems of ordinary differential equations (ODEs). Also, we show that every system of second-order ODEs can be expressed as a non-standard (partial) Hamiltonian system of first-order ODEs by introducing an artificial Hamiltonian. This notion of an artificial Hamiltonian gives a new way to solve dynamical systems of first-order ODEs and systems of second-order ODEs that can be expressed as a non-standard (partial) Hamiltonian system by using the known techniques applicable to the non-standard Hamiltonian systems. We employ the proposed notion to solve dynamical systems of first-order ODEs arising in epidemics.

3 - LINEAR SYSTEMS

The material presented below amounts roughly to the so-called fundamental theorems of Lie and their implications concerning Lie algebras, Lie subgroups and subalgebras. By straightforward and fairly elementary steps, we shall extend the concept of Lie group to include groups admitting coordinate systems in which the functions defining ab, the group product of elements a, b, possess continuous first order derivatives in the coordinates of a and satisfy a Lipschitz condition in those of b. That such groups are equivalent to classical (i.e. analytical) Lie groups was announced in 1936 by van Kampen, although apparently van Kampen published nothing by way of proof beyond a certain decisive uniqueness theorem concerning systems of ordinary differential equations. (This is contained in our theorem (12.1); the proof given below is essentially that of van Kampen [3]). It is conceivable that the condition relative to coordinates of b could be weakened or dispensed with entirely. For, the only part the condition in question plays in our development is to make it certain that the inverse of a certain mapping-the mapping into canonical coordinates-is single-valued. In this connection mention should be made of the paper [1J of Garrett Birkhoff in which it is shown that the existence of continuous first order derivatives of ab with respect to the coordinates of both the a's and the b's is sufficient for defining a Lie group. We shall consider only groups with a finite number of real parameters (whereas the groups considered by Birkhoff are not necessarily finite dimensional). This makes it possible to use the standard existence theorems for systems of ordinary real differential equations. We make no use whatever of the theory of partial differential equations. We maintain throughout a purely local point of view in the sense that we consider only what happens in the neighborhood of the identity. Historical notes. The proof of the uniqueness of square roots (8.5) is due to Claude Chevalley (unpublished). The proof given below of Lie's theorem that a linear system of vector fields which is closed under commutation defines a Lie group is, we believe, due to van der Waerden [4]. We have supplied a number of preliminary lemmas which make the proof rigorously applicable to the case in which the vector fields are only assumed to possess continuous first order derivatives. The theorem in ?17 that every local subgroup of a Lie group is a Lie group is due essentially to E. Cartan [2]. 481

71 - Integrable Combinations

10 - Applications of Systems of Ordinary Differential Equations

A Hopfield neural network is described by a system of nonlinear ordinary differential equations. We develop a broad range of numerical schemes that are applicable for a wide range of computational problems. We review here our study on an approximate solution of the Fredholm integral equation, and linear and nonlinear singular and hypersingular integral equations, using a continuous method for solving operator equations. This method assumes that the original system is associated with a Cauchy problem for systems of ordinary differential equations on Hopfield neural networks. We present sufficient conditions for the Hopfield networks’ stability defined via coefficients of systems of differential equations.

For each associative, commutative, and unitary algebra over the field of real or complex numbers and an integrable nonlinear ordinary differential equation we can to construct integrable systems of ordinary differential equations and integrable systems of partial differential equations. In this paper we consider in some sense the inverse problem. Determine the conditions under which a given system of ordinary differential equations can be represented as a differential equation in some associative, commutative and unitary algebra. It is also shown that associativity is not a necessary condition.

Analysis and Control of Time Delayed Systems via the Lambert W Function

Using the Lie symmetry approach, we have examined herein the system of partial differential equations corresponding to the Einstein–Maxwell equations for a static axially symmetric spacetime. The method used reduces the system of partial differential equations to a system of ordinary differential equations according to the Lie symmetry admitted. In particular, we found the relevant system of ordinary differential equations is all optimal subgroups. The system of ordinary differential equations is further solved in general to obtain exact solutions. Several new physically important families of exact solutions are derived.

A toxin-producing phytoplankton and zooplankton system is investigated. Considering that zooplankton can be harvested for food in some bodies of water, the harvesting term is introduced to zooplankton population. Firstly, from the ordinary differential equation (ODE) system, we obtain the global asymptotic stability of equilibrium and optimal capture problem. Secondly, based on the ODE system, the diffusion term is introduced and the global asymptotic stability of the steady state solution is obtained. As a result, the diffusion cannot affect the global asymptotic stability of equilibrium, and Turing instability cannot occur. Once again, a delayed differential equation (DDE) system is put forward. The global asymptotic stability of boundary equilibrium and the existence of local Hopf bifurcation at positive equilibrium are discussed. Furthermore, it is proved that there exists at least one positive periodic solution as delay varies in some region by using the global Hopf result of Wu for functional differential equations. Lastly, some numerical simulations are carried out for supporting the theoretical analyses and the positive impacts of harvesting effort, and the release rate of toxin is given. The unstable interval of the positive equilibrium becomes smaller and smaller with the increase of harvesting effort or the release rate of toxin.

<div class="section abstract"><div class="htmlview paragraph">In this paper, a new solution method is presented to study the effect of wave propagation in engine manifolds, which includes solving one-dimensional models for compressible flow of air. Velocity, pressure, and density profiles are found by solving a system of non-linear Partial Differential Equations (PDEs) in space and time derived from Euler’s equations. The 1D model includes frictional losses, area change, and heat transfer. The solution is traditionally found by utilizing the Method of Characteristics and applying finite difference solutions to the resulting system of ordinary differential equations (ODEs) over a discretized grid. In this work, orthogonal collocation is used to solve the system of ODEs that is defined along the characteristic curves. Orthogonal polynomials are utilized to approximate velocity, pressure, sound speed, and the characteristic curves along which the system of PDEs reduce to a system of ODEs. The approximation polynomials are defined over the whole manifold domain, transforming Euler’s equations into a system of ODEs that can be solved using a generic ODE solver. This reduction is done symbolically using a computer algebra system (Maple). The method results in a system of ODEs that has a higher spatial order along the whole space compared to methods based on finite differences, reducing the number of nodes required to find an acceptable solution that captures the state dynamics at different locations inside the manifold in real time. The proposed model is compared against the Method of Characteristics (MOC) that is used as a reference model; this comparison includes the states at the inlet, outlet, and midpoint. In summary, a high order method that can calculate solutions of the 1D manifold model equations is developed by finding the respective polynomial approximations along the 1D space and solving the resultant system using a generic solver in real time.</div></div>

The brief presentation of the method of fast expansions is given to solve nonlinear differential equations. Application rules of the operator of fast expansions are specified for solving differential equations. According to the method of fast expansions, an unknown function can be represented as the sum of the boundary function and Fourier series sines and cosines for one variable. The special construction of the boundary functions leads to reasonably fast convergence of the Fourier series, so that for engineering calculations, it is sufficient to consider only the first three members. The method is applicable both to linear and nonlinear integro-differential systems. By means of applying the method of fast expansions to nonlinear Navier-Stokes equations the problem is reduced to a closed system of ordinary differential equations, which solution doesn't represent special difficulties. We can reapply the method of fast expansions to the resulting system of differential equations and reduce the original problem to a system of algebraic equations. If the problem is n -dimensional, then after n -fold application of the method of fast expansions the problem will be reduced to a closed algebraic system. Finally, we obtain an analytic-form solution of complicated boundary value problem in partial derivatives. The flow of an incompressible viscous fluid of Navier–Stokes is considered in a curvilinear pipe. The problem is reduced to solving a closed system of ordinary differential equations with boundary conditions by the method of fast expansions. The article considers peculiarities of finding the coefficients of boundary functions and Fourier coefficients for the zero-order and first-order operators of fast expansions. Obtaining the analytic-form solution is of great interest, because it allows to analyze and to investigate the influence of various factors on the properties of the viscous fluid in specific cases.

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.