Sample records for ordinary differential equation
-
Solving Ordinary Differential Equations
Krogh, F. T.
1987-01-01
Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.
-
Ordinary differential equations
Greenberg, Michael D
2014-01-01
Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps
-
Introduction to ordinary differential equations
Rabenstein, Albert L
1966-01-01
Introduction to Ordinary Differential Equations is a 12-chapter text that describes useful elementary methods of finding solutions using ordinary differential equations. This book starts with an introduction to the properties and complex variable of linear differential equations. Considerable chapters covered topics that are of particular interest in applications, including Laplace transforms, eigenvalue problems, special functions, Fourier series, and boundary-value problems of mathematical physics. Other chapters are devoted to some topics that are not directly concerned with finding solutio
-
Ordinary differential equations
Miller, Richard K
1982-01-01
Ordinary Differential Equations is an outgrowth of courses taught for a number of years at Iowa State University in the mathematics and the electrical engineering departments. It is intended as a text for a first graduate course in differential equations for students in mathematics, engineering, and the sciences. Although differential equations is an old, traditional, and well-established subject, the diverse backgrounds and interests of the students in a typical modern-day course cause problems in the selection and method of presentation of material. In order to compensate for this diversity,
-
Dichotomies for generalized ordinary differential equations and applications
Bonotto, E. M.; Federson, M.; Santos, F. L.
2018-03-01
In this work we establish the theory of dichotomies for generalized ordinary differential equations, introducing the concepts of dichotomies for these equations, investigating their properties and proposing new results. We establish conditions for the existence of exponential dichotomies and bounded solutions. Using the correspondences between generalized ordinary differential equations and other equations, we translate our results to measure differential equations and impulsive differential equations. The fact that we work in the framework of generalized ordinary differential equations allows us to manage functions with many discontinuities and of unbounded variation.
-
Schwarz maps of algebraic linear ordinary differential equations
Sanabria Malagón, Camilo
2017-12-01
A linear ordinary differential equation is called algebraic if all its solution are algebraic over its field of definition. In this paper we solve the problem of finding closed form solution to algebraic linear ordinary differential equations in terms of standard equations. Furthermore, we obtain a method to compute all algebraic linear ordinary differential equations with rational coefficients by studying their associated Schwarz map through the Picard-Vessiot Theory.
-
Ordinary differential equations with applications in molecular biology.
Ilea, M; Turnea, M; Rotariu, M
2012-01-01
Differential equations are of basic importance in molecular biology mathematics because many biological laws and relations appear mathematically in the form of a differential equation. In this article we presented some applications of mathematical models represented by ordinary differential equations in molecular biology. The vast majority of quantitative models in cell and molecular biology are formulated in terms of ordinary differential equations for the time evolution of concentrations of molecular species. Assuming that the diffusion in the cell is high enough to make the spatial distribution of molecules homogenous, these equations describe systems with many participating molecules of each kind. We propose an original mathematical model with small parameter for biological phospholipid pathway. All the equations system includes small parameter epsilon. The smallness of epsilon is relative to the size of the solution domain. If we reduce the size of the solution region the same small epsilon will result in a different condition number. It is clear that the solution for a smaller region is less difficult. We introduce the mathematical technique known as boundary function method for singular perturbation system. In this system, the small parameter is an asymptotic variable, different from the independent variable. In general, the solutions of such equations exhibit multiscale phenomena. Singularly perturbed problems form a special class of problems containing a small parameter which may tend to zero. Many molecular biology processes can be quantitatively characterized by ordinary differential equations. Mathematical cell biology is a very active and fast growing interdisciplinary area in which mathematical concepts, techniques, and models are applied to a variety of problems in developmental medicine and bioengineering. Among the different modeling approaches, ordinary differential equations (ODE) are particularly important and have led to significant advances
-
Differential equations a dynamical systems approach ordinary differential equations
Hubbard, John H
1991-01-01
This is a corrected third printing of the first part of the text Differential Equations: A Dynamical Systems Approach written by John Hubbard and Beverly West. The authors' main emphasis in this book is on ordinary differential equations. The book is most appropriate for upper level undergraduate and graduate students in the fields of mathematics, engineering, and applied mathematics, as well as the life sciences, physics and economics. Traditional courses on differential equations focus on techniques leading to solutions. Yet most differential equations do not admit solutions which can be written in elementary terms. The authors have taken the view that a differential equations defines functions; the object of the theory is to understand the behavior of these functions. The tools the authors use include qualitative and numerical methods besides the traditional analytic methods. The companion software, MacMath, is designed to bring these notions to life.
-
Monograph - The Numerical Integration of Ordinary Differential Equations.
Hull, T. E.
The materials presented in this monograph are intended to be included in a course on ordinary differential equations at the upper division level in a college mathematics program. These materials provide an introduction to the numerical integration of ordinary differential equations, and they can be used to supplement a regular text on this…
-
Schiesser, William E
2014-01-01
Features a solid foundation of mathematical and computational tools to formulate and solve real-world ODE problems across various fields With a step-by-step approach to solving ordinary differential equations (ODEs), Differential Equation Analysis in Biomedical Science and Engineering: Ordinary Differential Equation Applications with R successfully applies computational techniques for solving real-worldODE problems that are found in a variety of fields, including chemistry, physics, biology,and physiology. The book provides readers with the necessary knowledge to reproduce andextend the comp
-
Ordinary differential equation for local accumulation time.
Berezhkovskii, Alexander M
2011-08-21
Cell differentiation in a developing tissue is controlled by the concentration fields of signaling molecules called morphogens. Formation of these concentration fields can be described by the reaction-diffusion mechanism in which locally produced molecules diffuse through the patterned tissue and are degraded. The formation kinetics at a given point of the patterned tissue can be characterized by the local accumulation time, defined in terms of the local relaxation function. Here, we show that this time satisfies an ordinary differential equation. Using this equation one can straightforwardly determine the local accumulation time, i.e., without preliminary calculation of the relaxation function by solving the partial differential equation, as was done in previous studies. We derive this ordinary differential equation together with the accompanying boundary conditions and demonstrate that the earlier obtained results for the local accumulation time can be recovered by solving this equation. © 2011 American Institute of Physics
-
Inverse problems in ordinary differential equations and applications
Llibre, Jaume
2016-01-01
This book is dedicated to study the inverse problem of ordinary differential equations, that is it focuses in finding all ordinary differential equations that satisfy a given set of properties. The Nambu bracket is the central tool in developing this approach. The authors start characterizing the ordinary differential equations in R^N which have a given set of partial integrals or first integrals. The results obtained are applied first to planar polynomial differential systems with a given set of such integrals, second to solve the 16th Hilbert problem restricted to generic algebraic limit cycles, third for solving the inverse problem for constrained Lagrangian and Hamiltonian mechanical systems, fourth for studying the integrability of a constrained rigid body. Finally the authors conclude with an analysis on nonholonomic mechanics, a generalization of the Hamiltonian principle, and the statement an solution of the inverse problem in vakonomic mechanics.
-
From ordinary to partial differential equations
Esposito, Giampiero
2017-01-01
This book is addressed to mathematics and physics students who want to develop an interdisciplinary view of mathematics, from the age of Riemann, Poincaré and Darboux to basic tools of modern mathematics. It enables them to acquire the sensibility necessary for the formulation and solution of difficult problems, with an emphasis on concepts, rigour and creativity. It consists of eight self-contained parts: ordinary differential equations; linear elliptic equations; calculus of variations; linear and non-linear hyperbolic equations; parabolic equations; Fuchsian functions and non-linear equations; the functional equations of number theory; pseudo-differential operators and pseudo-differential equations. The author leads readers through the original papers and introduces new concepts, with a selection of topics and examples that are of high pedagogical value.
-
A textbook on ordinary differential equations
Ahmad, Shair
2015-01-01
This book offers readers a primer on the theory and applications of Ordinary Differential Equations. The style used is simple, yet thorough and rigorous. Each chapter ends with a broad set of exercises that range from the routine to the more challenging and thought-provoking. Solutions to selected exercises can be found at the end of the book. The book contains many interesting examples on topics such as electric circuits, the pendulum equation, the logistic equation, the Lotka-Volterra system, the Laplace Transform, etc., which introduce students to a number of interesting aspects of the theory and applications. The work is mainly intended for students of Mathematics, Physics, Engineering, Computer Science and other areas of the natural and social sciences that use ordinary differential equations, and who have a firm grasp of Calculus and a minimal understanding of the basic concepts used in Linear Algebra. It also studies a few more advanced topics, such as Stability Theory and Boundary Value Problems, whic...
-
Numerical analysis of systems of ordinary and stochastic differential equations
Artemiev, S S
1997-01-01
This text deals with numerical analysis of systems of both ordinary and stochastic differential equations. It covers numerical solution problems of the Cauchy problem for stiff ordinary differential equations (ODE) systems by Rosenbrock-type methods (RTMs).
-
Ordinary differential equations
Cox, William
1995-01-01
Building on introductory calculus courses, this text provides a sound foundation in the underlying principles of ordinary differential equations. Important concepts, including uniqueness and existence theorems, are worked through in detail and the student is encouraged to develop much of the routine material themselves, thus helping to ensure a solid understanding of the fundamentals required.The wide use of exercises, problems and self-assessment questions helps to promote a deeper understanding of the material and it is developed in such a way that it lays the groundwork for further
-
Kurzweil, J
1986-01-01
The author, Professor Kurzweil, is one of the world's top experts in the area of ordinary differential equations - a fact fully reflected in this book. Unlike many classical texts which concentrate primarily on methods of integration of differential equations, this book pursues a modern approach: the topic is discussed in full generality which, at the same time, permits us to gain a deep insight into the theory and to develop a fruitful intuition. The basic framework of the theory is expanded by considering further important topics like stability, dependence of a solution on a parameter, Car
-
Institute of Scientific and Technical Information of China (English)
WANG; Shunjin; ZHANG; Hua
2006-01-01
The problem of preserving fidelity in numerical computation of nonlinear ordinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure of the dynamical system.The ordinary differential equations are lifted to the corresponding partial differential equations in the framework of algebraic dynamics,and a new algorithm-algebraic dynamics algorithm is proposed based on the exact analytical solutions of the ordinary differential equations by the algebraic dynamics method.In the new algorithm,the time evolution of the ordinary differential system is described locally by the time translation operator and globally by the time evolution operator.The exact analytical piece-like solution of the ordinary differential equations is expressd in terms of Taylor series with a local convergent radius,and its finite order truncation leads to the new numerical algorithm with a controllable precision better than Runge Kutta Algorithm and Symplectic Geometric Algorithm.
-
FORSIM, Solution of Ordinary or Partial Differential Equation with Initial Conditions
International Nuclear Information System (INIS)
Carver, M.B.
1985-01-01
1 - Description of problem or function: FORSIM is a FORTRAN oriented simulation program which automates the continuous transient solution of systems of ordinary and/or partial differential equations. The user writes his equations in a FORTRAN subroutine, following prescribed rules, and loads this routine along with the executive routines. The executive routines then read in initial data supplied by the user and proceed with the integration. 2 - Method of solution: Partial differential equations are converted to coupled ordinary differential equations by suitable discretization formulae. Integration is done by variable order, variable step-size error controlled algorithms. 3 - Restrictions on the complexity of the problem - Maximum of: 1000 ordinary differential equations
-
Solutions manual to accompany Ordinary differential equations
Greenberg, Michael D
2014-01-01
Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps
-
A course in ordinary differential equations
Swift, Randall J
2014-01-01
Praise for the First Edition:"A Course in Ordinary Differential Equations deserves to be on the MAA's Basic Library List … the book with its layout, is very student friendly-it is easy to read and understand; every chapter and explanations flow smoothly and coherently … the reviewer would recommend this book highly for undergraduate introductory differential equation courses." -Srabasti Dutta, College of Saint Elizabeth, MAA Online, July 2008"An important feature is that the exposition is richly accompanied by computer algebra code (equally distributed between MATLAB, Mathematica, and Maple). The major part of the book is devoted to classical theory (both for systems and higher order equations). The necessary material from linear algebra is also covered. More advanced topics include numerical methods, stability of equilibria, bifurcations, Laplace transforms, and the power series method."-EMS Newsletter, June 2007"This is a delightful textbook for a standard one-semester undergraduate course in ordinary d...
-
A Unified Introduction to Ordinary Differential Equations
Lutzer, Carl V.
2006-01-01
This article describes how a presentation from the point of view of differential operators can be used to (partially) unify the myriad techniques in an introductory course in ordinary differential equations by providing students with a powerful, flexible paradigm that extends into (or from) linear algebra. (Contains 1 footnote.)
-
Time-course window estimator for ordinary differential equations linear in the parameters
Vujacic, Ivan; Dattner, Itai; Gonzalez, Javier; Wit, Ernst
In many applications obtaining ordinary differential equation descriptions of dynamic processes is scientifically important. In both, Bayesian and likelihood approaches for estimating parameters of ordinary differential equations, the speed and the convergence of the estimation procedure may
-
International Nuclear Information System (INIS)
Ka-Lin, Su; Yuan-Xi, Xie
2010-01-01
By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ordinary differential equation method is presented for solving the (2 + 1)-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the (2 + 1)-dimensional sine-Poisson equation are derived in a simple manner by this technique. (general)
-
Network Reconstruction From High-Dimensional Ordinary Differential Equations.
Chen, Shizhe; Shojaie, Ali; Witten, Daniela M
2017-01-01
We consider the task of learning a dynamical system from high-dimensional time-course data. For instance, we might wish to estimate a gene regulatory network from gene expression data measured at discrete time points. We model the dynamical system nonparametrically as a system of additive ordinary differential equations. Most existing methods for parameter estimation in ordinary differential equations estimate the derivatives from noisy observations. This is known to be challenging and inefficient. We propose a novel approach that does not involve derivative estimation. We show that the proposed method can consistently recover the true network structure even in high dimensions, and we demonstrate empirical improvement over competing approaches. Supplementary materials for this article are available online.
-
Ordinary differential equations principles and applications
Nandakumaran, A K; George, Raju K
2017-01-01
Written in a clear, logical and concise manner, this comprehensive resource allows students to quickly understand the key principles, techniques and applications of ordinary differential equations. Important topics including first and second order linear equations, initial value problems and qualitative theory are presented in separate chapters. The concepts of two point boundary value problems, physical models and first order partial differential equations are discussed in detail. The text uses tools of calculus and real analysis to get solutions in explicit form. While discussing first order linear systems, linear algebra techniques are used. The real-life applications are interspersed throughout the book to invoke reader's interest. The methods and tricks to solve numerous mathematical problems with sufficient derivations and explanation are provided. The proofs of theorems are explained for the benefit of the readers.
-
Hojman's theorem of the third-order ordinary differential equation
International Nuclear Information System (INIS)
Hong-Sheng, Lü; Hong-Bin, Zhang; Shu-Long, Gu
2009-01-01
This paper extends Hojman's conservation law to the third-order differential equation. A new conserved quantity is constructed based on the Lie group of transformation generators of the equations of motion. The generators contain variations of the time and generalized coordinates. Two independent non-trivial conserved quantities of the third-order ordinary differential equation are obtained. A simple example is presented to illustrate the applications of the results. (general)
-
Lie group classification of first-order delay ordinary differential equations
Dorodnitsyn, Vladimir A.; Kozlov, Roman; Meleshko, Sergey V.; Winternitz, Pavel
2018-05-01
A group classification of first-order delay ordinary differential equations (DODEs) accompanied by an equation for the delay parameter (delay relation) is presented. A subset of such systems (delay ordinary differential systems or DODSs), which consists of linear DODEs and solution-independent delay relations, have infinite-dimensional symmetry algebras—as do nonlinear ones that are linearizable by an invertible transformation of variables. Genuinely nonlinear DODSs have symmetry algebras of dimension n, . It is shown how exact analytical solutions of invariant DODSs can be obtained using symmetry reduction.
-
A new numerical approximation of the fractal ordinary differential equation
Atangana, Abdon; Jain, Sonal
2018-02-01
The concept of fractal medium is present in several real-world problems, for instance, in the geological formation that constitutes the well-known subsurface water called aquifers. However, attention has not been quite devoted to modeling for instance, the flow of a fluid within these media. We deem it important to remind the reader that the concept of fractal derivative is not to represent the fractal sharps but to describe the movement of the fluid within these media. Since this class of ordinary differential equations is highly complex to solve analytically, we present a novel numerical scheme that allows to solve fractal ordinary differential equations. Error analysis of the method is also presented. Application of the method and numerical approximation are presented for fractal order differential equation. The stability and the convergence of the numerical schemes are investigated in detail. Also some exact solutions of fractal order differential equations are presented and finally some numerical simulations are presented.
-
Algorithmic Verification of Linearizability for Ordinary Differential Equations
Lyakhov, Dmitry A.
2017-07-19
For a nonlinear ordinary differential equation solved with respect to the highest order derivative and rational in the other derivatives and in the independent variable, we devise two algorithms to check if the equation can be reduced to a linear one by a point transformation of the dependent and independent variables. The first algorithm is based on a construction of the Lie point symmetry algebra and on the computation of its derived algebra. The second algorithm exploits the differential Thomas decomposition and allows not only to test the linearizability, but also to generate a system of nonlinear partial differential equations that determines the point transformation and the coefficients of the linearized equation. The implementation of both algorithms is discussed and their application is illustrated using several examples.
-
Ordinary differential equations a graduate text
Bhamra, K S
2015-01-01
ORDINARY DIFFERENTIAL EQUATIONS: A Graduate Text presents a systematic and comprehensive introduction to ODEs for graduate and postgraduate students. The systematic organized text on differential inequalities, Gronwall's inequality, Nagumo's theorems, Osgood's criteria and applications of different equations of first order is dealt with in a greater depth. The book discusses qualitative and quantitative aspects of the Strum - Liouville problems, Green's function, integral equations, Laplace transform and is supported by a number of worked-out examples in each lesson to make the concepts clear. A lot of stress on stability theory is laid down, especially on Lyapunov and Poincare stability theory. A numerous figures in various lessons (in particular lessons dealing with stability theory) have been added to clarify the key concepts in DE theory. Nonlinear oscillation in conservative systems and Hamiltonian systems highlights basic nature of the systems considered. Perturbation techniques lesson deals in fairly d...
-
Generation and Identification of Ordinary Differential Equations of Maximal Symmetry Algebra
Directory of Open Access Journals (Sweden)
J. C. Ndogmo
2016-01-01
Full Text Available An effective method for generating linear ordinary differential equations of maximal symmetry in their most general form is found, and an explicit expression for the point transformation reducing the equation to its canonical form is obtained. New expressions for the general solution are also found, as well as several identification and other results and a direct proof of the fact that a linear ordinary differential equation is iterative if and only if it is reducible to the canonical form by a point transformation. New classes of solvable equations parameterized by an arbitrary function are also found, together with simple algebraic expressions for the corresponding general solution.
-
Low Dimensional Vessiot-Guldberg-Lie Algebras of Second-Order Ordinary Differential Equations
Directory of Open Access Journals (Sweden)
Rutwig Campoamor-Stursberg
2016-03-01
Full Text Available A direct approach to non-linear second-order ordinary differential equations admitting a superposition principle is developed by means of Vessiot-Guldberg-Lie algebras of a dimension not exceeding three. This procedure allows us to describe generic types of second-order ordinary differential equations subjected to some constraints and admitting a given Lie algebra as Vessiot-Guldberg-Lie algebra. In particular, well-known types, such as the Milne-Pinney or Kummer-Schwarz equations, are recovered as special cases of this classification. The analogous problem for systems of second-order differential equations in the real plane is considered for a special case that enlarges the generalized Ermakov systems.
-
Numerical solutions of ordinary and partial differential equations in the frequency domain
International Nuclear Information System (INIS)
Hazi, G.; Por, G.
1997-01-01
Numerical problems during the noise simulation in a nuclear power plant are discussed. The solutions of ordinary and partial differential equations are studied in the frequency domain. Numerical methods by the transfer function method are applied. It is shown that the correctness of the numerical methods is limited for ordinary differential equations in the frequency domain. To overcome the difficulties, step-size selection is suggested. (author)
-
Symmetries of th-Order Approximate Stochastic Ordinary Differential Equations
Fredericks, E.; Mahomed, F. M.
2012-01-01
Symmetries of $n$ th-order approximate stochastic ordinary differential equations (SODEs) are studied. The determining equations of these SODEs are derived in an Itô calculus context. These determining equations are not stochastic in nature. SODEs are normally used to model nature (e.g., earthquakes) or for testing the safety and reliability of models in construction engineering when looking at the impact of random perturbations.
-
Generalized Ordinary Differential Equation Models.
Miao, Hongyu; Wu, Hulin; Xue, Hongqi
2014-10-01
Existing estimation methods for ordinary differential equation (ODE) models are not applicable to discrete data. The generalized ODE (GODE) model is therefore proposed and investigated for the first time. We develop the likelihood-based parameter estimation and inference methods for GODE models. We propose robust computing algorithms and rigorously investigate the asymptotic properties of the proposed estimator by considering both measurement errors and numerical errors in solving ODEs. The simulation study and application of our methods to an influenza viral dynamics study suggest that the proposed methods have a superior performance in terms of accuracy over the existing ODE model estimation approach and the extended smoothing-based (ESB) method.
-
A textbook on ordinary differential equations
Ahmad, Shair
2014-01-01
The book is a primer of the theory of Ordinary Differential Equations. Each chapter is completed by a broad set of exercises; the reader will also find a set of solutions of selected exercises. The book contains many interesting examples as well (like the equations for the electric circuits, the pendium equation, the logistic equation, the Lotka-Volterra system, and many other) which introduce the reader to some interesting aspects of the theory and its applications. The work is mainly addressed to students of Mathematics, Physics, Engineering, Statistics, Computer Sciences, with knowledge of Calculus and Linear Algebra, and contains more advanced topics for further developments, such as Laplace transform; Stability theory and existence of solutions to Boundary Value problems. The authors are preparing a complete solutions manual, containing solutions to all the exercises published in the book. The manual will be available Summer 2014. Instructors who wish to adopt the book may request the manual by writing...
-
Design of TIR collimating lens for ordinary differential equation of extended light source
Zhan, Qianjing; Liu, Xiaoqin; Hou, Zaihong; Wu, Yi
2017-10-01
The source of LED has been widely used in our daily life. The intensity angle distribution of single LED is lambert distribution, which does not satisfy the requirement of people. Therefore, we need to distribute light and change the LED's intensity angle distribution. The most commonly method to change its intensity angle distribution is the free surface. Generally, using ordinary differential equations to calculate free surface can only be applied in a point source, but it will lead to a big error for the expand light. This paper proposes a LED collimating lens based on the ordinary differential equation, combined with the LED's light distribution curve, and adopt the method of calculating the center gravity of the extended light to get the normal vector. According to the law of Snell, the ordinary differential equations are constructed. Using the runge-kutta method for solution of ordinary differential equation solution, the curve point coordinates are gotten. Meanwhile, the edge point data of lens are imported into the optical simulation software TracePro. Based on 1mm×1mm single lambert body for light conditions, The degrees of collimating light can be close to +/-3. Furthermore, the energy utilization rate is higher than 85%. In this paper, the point light source is used to calculate partial differential equation method and compared with the simulation of the lens, which improve the effect of 1 degree of collimation.
-
On oscillation of second-order linear ordinary differential equations
Czech Academy of Sciences Publication Activity Database
Lomtatidze, A.; Šremr, Jiří
2011-01-01
Roč. 54, - (2011), s. 69-81 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z10190503 Keywords : linear second-order ordinary differential equation * Kamenev theorem * oscillation Subject RIV: BA - General Mathematics http://www.rmi.ge/jeomj/memoirs/vol54/abs54-4.htm
-
Hunt, L. R.; Villarreal, Ramiro
1987-01-01
System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.
-
Antiperiodic Boundary Value Problems for Second-Order Impulsive Ordinary Differential Equations
Directory of Open Access Journals (Sweden)
2009-02-01
Full Text Available We consider a second-order ordinary differential equation with antiperiodic boundary conditions and impulses. By using Schaefer's fixed-point theorem, some existence results are obtained.
-
Numerical integration of asymptotic solutions of ordinary differential equations
Thurston, Gaylen A.
1989-01-01
Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.
-
Rosenbaum, J. S.
1971-01-01
Systems of ordinary differential equations in which the magnitudes of the eigenvalues (or time constants) vary greatly are commonly called stiff. Such systems of equations arise in nuclear reactor kinetics, the flow of chemically reacting gas, dynamics, control theory, circuit analysis and other fields. The research reported develops an A-stable numerical integration technique for solving stiff systems of ordinary differential equations. The method, which is called the generalized trapezoidal rule, is a modification of the trapezoidal rule. However, the method is computationally more efficient than the trapezoidal rule when the solution of the almost-discontinuous segments is being calculated.
-
Ding, A Adam; Wu, Hulin
2014-10-01
We propose a new method to use a constrained local polynomial regression to estimate the unknown parameters in ordinary differential equation models with a goal of improving the smoothing-based two-stage pseudo-least squares estimate. The equation constraints are derived from the differential equation model and are incorporated into the local polynomial regression in order to estimate the unknown parameters in the differential equation model. We also derive the asymptotic bias and variance of the proposed estimator. Our simulation studies show that our new estimator is clearly better than the pseudo-least squares estimator in estimation accuracy with a small price of computational cost. An application example on immune cell kinetics and trafficking for influenza infection further illustrates the benefits of the proposed new method.
-
The relationship among the solutions of two auxiliary ordinary differential equations
International Nuclear Information System (INIS)
Liu Xiaoping; Liu Chunping
2009-01-01
In a recent article [Phys. Lett. A 356 (2006) 124], Sirendaoreji extended their auxiliary equation method by introducing a new auxiliary ordinary differential equation (NAODE) and its 14 solutions. Then the author studied some nonlinear evolution equations (NLEEs) and got more exact travelling wave solutions. In this paper, we will show that the 14 solutions of the NAODE are actually the same as the solutions obtained by original auxiliary equation method, and they are only different in the form.
-
Ordinary Differential Equation Models for Adoptive Immunotherapy.
Talkington, Anne; Dantoin, Claudia; Durrett, Rick
2018-05-01
Modified T cells that have been engineered to recognize the CD19 surface marker have recently been shown to be very successful at treating acute lymphocytic leukemias. Here, we explore four previous approaches that have used ordinary differential equations to model this type of therapy, compare their properties, and modify the models to address their deficiencies. Although the four models treat the workings of the immune system in slightly different ways, they all predict that adoptive immunotherapy can be successful to move a patient from the large tumor fixed point to an equilibrium with little or no tumor.
-
VODE, Variable Coefficient Ordinary Differential Equations (ODE) Solver
International Nuclear Information System (INIS)
Brown, P.N.; Hindmarsh, A.C.; Byrne, G.D.
2002-01-01
1 - Description of program or function: VODE is a package of subroutines for the numerical solution of the initial-value problem for systems of first-order ordinary differential equations. The package can be used for either stiff or non-stiff systems. In the stiff case, the Jacobian matrix is treated as full or banded. An algorithm is included for saving and reusing the Jacobian matrix under certain conditions. If storage is limited, this option may be suppressed. 2 - Method of solution - VODE uses the variable-order, variable- coefficient Adams-Moulton method for non-stiff systems and the variable-order, fixed-leading-coefficient Backward Differentiation Formula (BDF) method for stiff systems
-
A Simple Method to Find out when an Ordinary Differential Equation Is Separable
Cid, Jose Angel
2009-01-01
We present an alternative method to that of Scott (D. Scott, "When is an ordinary differential equation separable?", "Amer. Math. Monthly" 92 (1985), pp. 422-423) to teach the students how to discover whether a differential equation y[prime] = f(x,y) is separable or not when the nonlinearity f(x, y) is not explicitly factorized. Our approach is…
-
Error estimation in the neural network solution of ordinary differential equations.
Filici, Cristian
2010-06-01
In this article a method of error estimation for the neural approximation of the solution of an Ordinary Differential Equation is presented. Some examples of the application of the method support the theory presented. Copyright 2010. Published by Elsevier Ltd.
-
Approximate analytical methods for solving ordinary differential equations
Radhika, TSL; Rani, T Raja
2015-01-01
Approximate Analytical Methods for Solving Ordinary Differential Equations (ODEs) is the first book to present all of the available approximate methods for solving ODEs, eliminating the need to wade through multiple books and articles. It covers both well-established techniques and recently developed procedures, including the classical series solution method, diverse perturbation methods, pioneering asymptotic methods, and the latest homotopy methods.The book is suitable not only for mathematicians and engineers but also for biologists, physicists, and economists. It gives a complete descripti
-
Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations
Radhakrishnan, Krishnan; Hindmarsh, Alan C.
1993-01-01
LSODE, the Livermore Solver for Ordinary Differential Equations, is a package of FORTRAN subroutines designed for the numerical solution of the initial value problem for a system of ordinary differential equations. It is particularly well suited for 'stiff' differential systems, for which the backward differentiation formula method of orders 1 to 5 is provided. The code includes the Adams-Moulton method of orders 1 to 12, so it can be used for nonstiff problems as well. In addition, the user can easily switch methods to increase computational efficiency for problems that change character. For both methods a variety of corrector iteration techniques is included in the code. Also, to minimize computational work, both the step size and method order are varied dynamically. This report presents complete descriptions of the code and integration methods, including their implementation. It also provides a detailed guide to the use of the code, as well as an illustrative example problem.
-
Mallet, D. G.; McCue, S. W.
2009-01-01
The solution of linear ordinary differential equations (ODEs) is commonly taught in first-year undergraduate mathematics classrooms, but the understanding of the concept of a solution is not always grasped by students until much later. Recognizing what it is to be a solution of a linear ODE and how to postulate such solutions, without resorting to…
-
From Ordinary Differential Equations to Structural Causal Models: the deterministic case
Mooij, J.M.; Janzing, D.; Schölkopf, B.; Nicholson, A.; Smyth, P.
2013-01-01
We show how, and under which conditions, the equilibrium states of a first-order Ordinary Differential Equation (ODE) system can be described with a deterministic Structural Causal Model (SCM). Our exposition sheds more light on the concept of causality as expressed within the framework of
-
Solving Second-Order Ordinary Differential Equations without Using Complex Numbers
Kougias, Ioannis E.
2009-01-01
Ordinary differential equations (ODEs) is a subject with a wide range of applications and the need of introducing it to students often arises in the last year of high school, as well as in the early stages of tertiary education. The usual methods of solving second-order ODEs with constant coefficients, among others, rely upon the use of complex…
-
Runge-Kutta Methods for Linear Ordinary Differential Equations
Zingg, David W.; Chisholm, Todd T.
1997-01-01
Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.
-
Random ordinary differential equations and their numerical solution
Han, Xiaoying
2017-01-01
This book is intended to make recent results on the derivation of higher order numerical schemes for random ordinary differential equations (RODEs) available to a broader readership, and to familiarize readers with RODEs themselves as well as the closely associated theory of random dynamical systems. In addition, it demonstrates how RODEs are being used in the biological sciences, where non-Gaussian and bounded noise are often more realistic than the Gaussian white noise in stochastic differential equations (SODEs). RODEs are used in many important applications and play a fundamental role in the theory of random dynamical systems. They can be analyzed pathwise with deterministic calculus, but require further treatment beyond that of classical ODE theory due to the lack of smoothness in their time variable. Although classical numerical schemes for ODEs can be used pathwise for RODEs, they rarely attain their traditional order since the solutions of RODEs do not have sufficient smoothness to have Taylor ...
-
Institute of Scientific and Technical Information of China (English)
无
2009-01-01
In this paper, we consider a two-point boundary value problem for a system of second order ordinary differential equations. Under some conditions, we show the existence of positive solution to the system of second order ordinary differential equa-tions.
-
Khataybeh, S. N.; Hashim, I.
2018-04-01
In this paper, we propose for the first time a method based on Bernstein polynomials for solving directly a class of third-order ordinary differential equations (ODEs). This method gives a numerical solution by converting the equation into a system of algebraic equations which is solved directly. Some numerical examples are given to show the applicability of the method.
-
Stochastic Computational Approach for Complex Nonlinear Ordinary Differential Equations
International Nuclear Information System (INIS)
Khan, Junaid Ali; Raja, Muhammad Asif Zahoor; Qureshi, Ijaz Mansoor
2011-01-01
We present an evolutionary computational approach for the solution of nonlinear ordinary differential equations (NLODEs). The mathematical modeling is performed by a feed-forward artificial neural network that defines an unsupervised error. The training of these networks is achieved by a hybrid intelligent algorithm, a combination of global search with genetic algorithm and local search by pattern search technique. The applicability of this approach ranges from single order NLODEs, to systems of coupled differential equations. We illustrate the method by solving a variety of model problems and present comparisons with solutions obtained by exact methods and classical numerical methods. The solution is provided on a continuous finite time interval unlike the other numerical techniques with comparable accuracy. With the advent of neuroprocessors and digital signal processors the method becomes particularly interesting due to the expected essential gains in the execution speed. (general)
-
Camporesi, Roberto
2016-01-01
This book presents a method for solving linear ordinary differential equations based on the factorization of the differential operator. The approach for the case of constant coefficients is elementary, and only requires a basic knowledge of calculus and linear algebra. In particular, the book avoids the use of distribution theory, as well as the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The case of variable coefficients is addressed using Mammana’s result for the factorization of a real linear ordinary differential operator into a product of first-order (complex) factors, as well as a recent generalization of this result to the case of complex-valued coefficients.
-
Martini, Ruud; Kersten, P.H.M.
1983-01-01
Using 1-1 mappings, the complete symmetry groups of contact transformations of general linear second-order ordinary differential equations are determined from two independent solutions of those equations, and applied to the harmonic oscillator with and without damping.
-
Mohammed, Ahmed; Zeleke, Aklilu
2015-01-01
We introduce a class of second-order ordinary differential equations (ODEs) with variable coefficients whose closed-form solutions can be obtained by the same method used to solve ODEs with constant coefficients. General solutions for the homogeneous case are discussed.
-
Camacho-Machín, M.; Guerrero-Ortiz, C.
2015-01-01
The aim of this paper is to present and discuss some of the evidence regarding the resources that students use when they establish relationships between a contextual situation and an ordinary differential equation (ODE). We present research results obtained from work by seven students in a graduate level course in mathematics education, where they…
-
Blow up of solutions to ordinary differential equations arising in nonlinear dispersive problems
Directory of Open Access Journals (Sweden)
Milena Dimova
2018-03-01
Full Text Available We study a new class of ordinary differential equations with blow up solutions. Necessary and sufficient conditions for finite blow up time are proved. Based on the new differential equation, a revised version of the concavity method of Levine is proposed. As an application we investigate the non-existence of global solutions to the Cauchy problem of Klein-Gordon, and to the double dispersive equations. We obtain necessary and sufficient condition for finite time blow up with arbitrary positive energy. A very general sufficient condition for blow up is also given.
-
KRYSI, Ordinary Differential Equations Solver with Sdirk Krylov Method
International Nuclear Information System (INIS)
Hindmarsh, A.C.; Norsett, S.P.
2001-01-01
1 - Description of program or function: KRYSI is a set of FORTRAN subroutines for solving ordinary differential equations initial value problems. It is suitable for both stiff and non-stiff systems. When solving the implicit stage equations in the stiff case, KRYSI uses a Krylov subspace iteration method called the SPIGMR (Scaled Preconditioned Incomplete Generalized Minimum Residual) method. No explicit Jacobian storage is required, except where used in pre- conditioning. A demonstration problem is included with a description of two pre-conditioners that are natural for its solution by KRYSI. 2 - Method of solution: KRYSI uses a three-stage, third-order singly diagonally implicit Runge-Kutta (SDIRK) method. In the stiff case, a preconditioned Krylov subspace iteration within a (so-called) inexact Newton iteration is used to solve the system of nonlinear algebraic equations
-
Ordinary differential equations basics and beyond
Schaeffer, David G
2016-01-01
This book develops the theory of ordinary differential equations (ODEs), starting from an introductory level (with no prior experience in ODEs assumed) through to a graduate-level treatment of the qualitative theory, including bifurcation theory (but not chaos). While proofs are rigorous, the exposition is reader-friendly, aiming for the informality of face-to-face interactions. A unique feature of this book is the integration of rigorous theory with numerous applications of scientific interest. Besides providing motivation, this synthesis clarifies the theory and enhances scientific literacy. Other features include: (i) a wealth of exercises at various levels, along with commentary that explains why they matter; (ii) figures with consistent color conventions to identify nullclines, periodic orbits, stable and unstable manifolds; and (iii) a dedicated website with software templates, problem solutions, and other resources supporting the text. Given its many applications, the book may be used comfortably in sc...
-
Vujačić, Ivan; Dattner, Itai
In this paper we use the sieve framework to prove consistency of the ‘direct integral estimator’ of parameters for partially observed systems of ordinary differential equations, which are commonly used for modeling dynamic processes.
-
A variational approach to parameter estimation in ordinary differential equations
Directory of Open Access Journals (Sweden)
Kaschek Daniel
2012-08-01
Full Text Available Abstract Background Ordinary differential equations are widely-used in the field of systems biology and chemical engineering to model chemical reaction networks. Numerous techniques have been developed to estimate parameters like rate constants, initial conditions or steady state concentrations from time-resolved data. In contrast to this countable set of parameters, the estimation of entire courses of network components corresponds to an innumerable set of parameters. Results The approach presented in this work is able to deal with course estimation for extrinsic system inputs or intrinsic reactants, both not being constrained by the reaction network itself. Our method is based on variational calculus which is carried out analytically to derive an augmented system of differential equations including the unconstrained components as ordinary state variables. Finally, conventional parameter estimation is applied to the augmented system resulting in a combined estimation of courses and parameters. Conclusions The combined estimation approach takes the uncertainty in input courses correctly into account. This leads to precise parameter estimates and correct confidence intervals. In particular this implies that small motifs of large reaction networks can be analysed independently of the rest. By the use of variational methods, elements from control theory and statistics are combined allowing for future transfer of methods between the two fields.
-
A variational approach to parameter estimation in ordinary differential equations.
Kaschek, Daniel; Timmer, Jens
2012-08-14
Ordinary differential equations are widely-used in the field of systems biology and chemical engineering to model chemical reaction networks. Numerous techniques have been developed to estimate parameters like rate constants, initial conditions or steady state concentrations from time-resolved data. In contrast to this countable set of parameters, the estimation of entire courses of network components corresponds to an innumerable set of parameters. The approach presented in this work is able to deal with course estimation for extrinsic system inputs or intrinsic reactants, both not being constrained by the reaction network itself. Our method is based on variational calculus which is carried out analytically to derive an augmented system of differential equations including the unconstrained components as ordinary state variables. Finally, conventional parameter estimation is applied to the augmented system resulting in a combined estimation of courses and parameters. The combined estimation approach takes the uncertainty in input courses correctly into account. This leads to precise parameter estimates and correct confidence intervals. In particular this implies that small motifs of large reaction networks can be analysed independently of the rest. By the use of variational methods, elements from control theory and statistics are combined allowing for future transfer of methods between the two fields.
-
Using trees to compute approximate solutions to ordinary differential equations exactly
Grossman, Robert
1991-01-01
Some recent work is reviewed which relates families of trees to symbolic algorithms for the exact computation of series which approximate solutions of ordinary differential equations. It turns out that the vector space whose basis is the set of finite, rooted trees carries a natural multiplication related to the composition of differential operators, making the space of trees an algebra. This algebraic structure can be exploited to yield a variety of algorithms for manipulating vector fields and the series and algebras they generate.
-
Numerical methods for the solution of ordinary differential equations
International Nuclear Information System (INIS)
Azeem, M.
1999-01-01
The ode 113 code solves non-stiff differential equations and is a fully variable step, variable order, PECE implementation in terms of modified divided differences of Adams-Bashforth-Moulton family of formulas of order 1-12. The main objectives of this project were to modify PECE mode of ode 113 into PEC mode, study the variable step size and variable order strategy of both the modes and finally, develop the switching strategy between both PECE and PEC modes to minimize the cost of solving the ordinary differential equations. Using some test problems (including stiff, mild stiff and non-stiff), it was found that the PEC mode was more efficient for non-stiff problems at crude and intermediate tolerances and the PECE mode for all problems at the stringent tolerance. An automatic switching strategy was developed using the results observed from the step size and order plots of all the test problems for both the modes and gave the optimum results. (author)
-
Numerical solution of ordinary differential equations
Fox, L
1987-01-01
Nearly 20 years ago we produced a treatise (of about the same length as this book) entitled Computing methods for scientists and engineers. It was stated that most computation is performed by workers whose mathematical training stopped somewhere short of the 'professional' level, and that some books are therefore needed which use quite simple mathematics but which nevertheless communicate the essence of the 'numerical sense' which is exhibited by the real computing experts and which is surely needed, at least to some extent, by all who use modern computers and modern numerical software. In that book we treated, at no great length, a variety of computational problems in which the material on ordinary differential equations occupied about 50 pages. At that time it was quite common to find books on numerical analysis, with a little on each topic ofthat field, whereas today we are more likely to see similarly-sized books on each major topic: for example on numerical linear algebra, numerical approximation, numeri...
-
The qualitative theory of ordinary differential equations an introduction
Brauer, Fred
1989-01-01
""This is a very good book ... with many well-chosen examples and illustrations."" - American Mathematical MonthlyThis highly regarded text presents a self-contained introduction to some important aspects of modern qualitative theory for ordinary differential equations. It is accessible to any student of physical sciences, mathematics or engineering who has a good knowledge of calculus and of the elements of linear algebra. In addition, algebraic results are stated as needed; the less familiar ones are proved either in the text or in appendixes.The topics covered in the first three chapters a
-
Maat, Siti Mistima; Zakaria, Effandi
2011-01-01
Ordinary differential equations (ODEs) are one of the important topics in engineering mathematics that lead to the understanding of technical concepts among students. This study was conducted to explore the students' understanding of ODEs when they solve ODE questions using a traditional method as well as a computer algebraic system, particularly…
-
Camporesi, Roberto
2011-01-01
We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of…
-
Rosenbaum, J. S.
1976-01-01
If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.
-
Sánchez Pérez, J F; Conesa, M; Alhama, I; Alhama, F; Cánovas, M
2017-01-01
Classical dimensional analysis and nondimensionalization are assumed to be two similar approaches in the search for dimensionless groups. Both techniques, simplify the study of many problems. The first approach does not need to know the mathematical model, being sufficient a deep understanding of the physical phenomenon involved, while the second one begins with the governing equations and reduces them to their dimensionless form by simple mathematical manipulations. In this work, a formal protocol is proposed for applying the nondimensionalization process to ordinary differential equations, linear or not, leading to dimensionless normalized equations from which the resulting dimensionless groups have two inherent properties: In one hand, they are physically interpreted as balances between counteracting quantities in the problem, and on the other hand, they are of the order of magnitude unity. The solutions provided by nondimensionalization are more precise in every case than those from dimensional analysis, as it is illustrated by the applications studied in this work.
-
SAHYB-2, Solution of Ordinary Differential Equation with User-Supplied Subroutine
International Nuclear Information System (INIS)
Hoop, H. d'; Monterosso, R.
1967-01-01
1 - Nature of physical problem solved: SAHYB-2 is a general purpose programme for the solution of ordinary differential equations. These are written in a user-supplied subroutine called DER, which uses notations very close to mathematical formulas. Special mathematical functions are included in the programme, as: Function generation, delay generation, steps, ramps and pulses, as well as a simplified standard output procedure - boundary value problems or parametric optimisation may be handled by iterations adding a subroutine called REPEAT. 2 - Method of solution: Integration is carried out by constant step fourth-order Runge-Kutta method, or by a fixed or variable step Adams-Moulton predictor corrector method. 3 - Restrictions on the complexity of the problem: Maximum 150 differential equations of the first order. Maximum 30 tables for function generator or delay lines
-
Calculus & ordinary differential equations
Pearson, David
1995-01-01
Professor Pearson's book starts with an introduction to the area and an explanation of the most commonly used functions. It then moves on through differentiation, special functions, derivatives, integrals and onto full differential equations. As with other books in the series the emphasis is on using worked examples and tutorial-based problem solving to gain the confidence of students.
-
Gainetdinova, A A; Gazizov, R K
2017-01-01
We suggest an algorithm for integrating systems of two second-order ordinary differential equations with four symmetries. In particular, if the admitted transformation group has two second-order differential invariants, the corresponding system can be integrated by quadratures using invariant representation and the operator of invariant differentiation. Otherwise, the systems reduce to partially uncoupled forms and can also be integrated by quadratures.
-
Directory of Open Access Journals (Sweden)
Yoshitsugu Takei
2015-01-01
Full Text Available Using two concrete examples, we discuss the multisummability of WKB solutions of singularly perturbed linear ordinary differential equations. Integral representations of solutions and a criterion for the multisummability based on the Cauchy-Heine transform play an important role in the proof.
-
International Nuclear Information System (INIS)
Basak, K C; Ray, P C; Bera, R K
2009-01-01
The aim of the present analysis is to apply the Adomian decomposition method and He's variational method for the approximate analytical solution of a nonlinear ordinary fractional differential equation. The solutions obtained by the above two methods have been numerically evaluated and presented in the form of tables and also compared with the exact solution. It was found that the results obtained by the above two methods are in excellent agreement with the exact solution. Finally, a surface plot of the approximate solutions of the fractional differential equation by the above two methods is drawn for 0≤t≤2 and 1<α≤2.
-
Solving ordinary differential equations by electrical analogy: a multidisciplinary teaching tool
Sanchez Perez, J. F.; Conesa, M.; Alhama, I.
2016-11-01
Ordinary differential equations are the mathematical formulation for a great variety of problems in science and engineering, and frequently, two different problems are equivalent from a mathematical point of view when they are formulated by the same equations. Students acquire the knowledge of how to solve these equations (at least some types of them) using protocols and strict algorithms of mathematical calculation without thinking about the meaning of the equation. The aim of this work is that students learn to design network models or circuits in this way; with simple knowledge of them, students can establish the association of electric circuits and differential equations and their equivalences, from a formal point of view, that allows them to associate knowledge of two disciplines and promote the use of this interdisciplinary approach to address complex problems. Therefore, they learn to use a multidisciplinary tool that allows them to solve these kinds of equations, even students of first course of engineering, whatever the order, grade or type of non-linearity. This methodology has been implemented in numerous final degree projects in engineering and science, e.g., chemical engineering, building engineering, industrial engineering, mechanical engineering, architecture, etc. Applications are presented to illustrate the subject of this manuscript.
-
Directory of Open Access Journals (Sweden)
SURE KÖME
2014-12-01
Full Text Available In this paper, we investigated the effect of Magnus Series Expansion Method on homogeneous stiff ordinary differential equations with different stiffness ratios. A Magnus type integrator is used to obtain numerical solutions of two different examples of stiff problems and exact and approximate results are tabulated. Furthermore, absolute error graphics are demonstrated in detail.
-
Quantifying uncertainty, variability and likelihood for ordinary differential equation models
LENUS (Irish Health Repository)
Weisse, Andrea Y
2010-10-28
Abstract Background In many applications, ordinary differential equation (ODE) models are subject to uncertainty or variability in initial conditions and parameters. Both, uncertainty and variability can be quantified in terms of a probability density function on the state and parameter space. Results The partial differential equation that describes the evolution of this probability density function has a form that is particularly amenable to application of the well-known method of characteristics. The value of the density at some point in time is directly accessible by the solution of the original ODE extended by a single extra dimension (for the value of the density). This leads to simple methods for studying uncertainty, variability and likelihood, with significant advantages over more traditional Monte Carlo and related approaches especially when studying regions with low probability. Conclusions While such approaches based on the method of characteristics are common practice in other disciplines, their advantages for the study of biological systems have so far remained unrecognized. Several examples illustrate performance and accuracy of the approach and its limitations.
-
On one two-point BVP for the fourth order linear ordinary differential equation
Czech Academy of Sciences Publication Activity Database
Mukhigulashvili, Sulkhan; Manjikashvili, M.
2017-01-01
Roč. 24, č. 2 (2017), s. 265-275 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : fourth order linear ordinary differential equations * two-point boundary value problems Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0077/gmj-2016-0077. xml
-
Czech Academy of Sciences Publication Activity Database
Lomtatidze, Alexander
2017-01-01
Roč. 24, č. 2 (2017), s. 241-263 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : nonlinear ordinary differential equations * periodic boundary value problem * solvability Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2017-0009/gmj-2017-0009. xml
-
Czech Academy of Sciences Publication Activity Database
Lomtatidze, Alexander
2017-01-01
Roč. 24, č. 2 (2017), s. 241-263 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : nonlinear ordinary differential equations * periodic boundary value problem * solvability Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2017-0009/gmj-2017-0009.xml
-
On one two-point BVP for the fourth order linear ordinary differential equation
Czech Academy of Sciences Publication Activity Database
Mukhigulashvili, Sulkhan; Manjikashvili, M.
2017-01-01
Roč. 24, č. 2 (2017), s. 265-275 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : fourth order linear ordinary differential equations * two-point boundary value problems Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0077/gmj-2016-0077.xml
-
Boyko, Vyacheslav M; Popovych, Roman O; Shapoval, Nataliya M
2013-01-01
Lie symmetries of systems of second-order linear ordinary differential equations with constant coefficients are exhaustively described over both the complex and real fields. The exact lower and upper bounds for the dimensions of the maximal Lie invariance algebras possessed by such systems are obtained using an effective algebraic approach.
-
A block Krylov subspace time-exact solution method for linear ordinary differential equation systems
Bochev, Mikhail A.
2013-01-01
We propose a time-exact Krylov-subspace-based method for solving linear ordinary differential equation systems of the form $y'=-Ay+g(t)$ and $y"=-Ay+g(t)$, where $y(t)$ is the unknown function. The method consists of two stages. The first stage is an accurate piecewise polynomial approximation of
-
International Nuclear Information System (INIS)
Lima, M.L.; Mignaco, J.A.
1985-01-01
It is shown that the rational power law potentials in the two-body radial Schoedinger equation admit a systematic treatment available from the classical theory of ordinary linear differential equations of the second order. The admissible potentials come into families evolved from equations having a fixed number of elementary singularities. As a consequence, relations are found and discussed among the several potentials in a family. (Author) [pt
-
The geometric approach to sets of ordinary differential equations and Hamiltonian dynamics
Estabrook, F. B.; Wahlquist, H. D.
1975-01-01
The calculus of differential forms is used to discuss the local integration theory of a general set of autonomous first order ordinary differential equations. Geometrically, such a set is a vector field V in the space of dependent variables. Integration consists of seeking associated geometric structures invariant along V: scalar fields, forms, vectors, and integrals over subspaces. It is shown that to any field V can be associated a Hamiltonian structure of forms if, when dealing with an odd number of dependent variables, an arbitrary equation of constraint is also added. Families of integral invariants are an immediate consequence. Poisson brackets are isomorphic to Lie products of associated CT-generating vector fields. Hamilton's variational principle follows from the fact that the maximal regular integral manifolds of a closed set of forms must include the characteristics of the set.
-
Efficient solution of ordinary differential equations modeling electrical activity in cardiac cells.
Sundnes, J; Lines, G T; Tveito, A
2001-08-01
The contraction of the heart is preceded and caused by a cellular electro-chemical reaction, causing an electrical field to be generated. Performing realistic computer simulations of this process involves solving a set of partial differential equations, as well as a large number of ordinary differential equations (ODEs) characterizing the reactive behavior of the cardiac tissue. Experiments have shown that the solution of the ODEs contribute significantly to the total work of a simulation, and there is thus a strong need to utilize efficient solution methods for this part of the problem. This paper presents how an efficient implicit Runge-Kutta method may be adapted to solve a complicated cardiac cell model consisting of 31 ODEs, and how this solver may be coupled to a set of PDE solvers to provide complete simulations of the electrical activity.
-
On the coupling of systems of hyperbolic conservation laws with ordinary differential equations
International Nuclear Information System (INIS)
Borsche, Raul; Colombo, Rinaldo M; Garavello, Mauro
2010-01-01
Motivated by applications to the piston problem, to a manhole model, to blood flow and to supply chain dynamics, this paper deals with a system of conservation laws coupled with a system of ordinary differential equations. The former is defined on a domain with boundary and the coupling is provided by the boundary condition. For each of the examples considered, numerical integrations are provided
-
Hou, Jie; Dong, Jianji; Zhang, Xinliang
2017-06-15
Systems of ordinary differential equations (SODEs) are crucial for describing the dynamic behaviors in various systems such as modern control systems which require observability and controllability. In this Letter, we propose and experimentally demonstrate an all-optical SODE solver based on the silicon-on-insulator platform. We use an add/drop microring resonator to construct two different ordinary differential equations (ODEs) and then introduce two external feedback waveguides to realize the coupling between these ODEs, thus forming the SODE solver. A temporal coupled mode theory is used to deduce the expression of the SODE. A system experiment is carried out for further demonstration. For the input 10 GHz NRZ-like pulses, the measured output waveforms of the SODE solver agree well with the calculated results.
-
Sun, Mengyang; Cheng, Xianrui; Socolar, Joshua E S
2013-06-01
A common approach to the modeling of gene regulatory networks is to represent activating or repressing interactions using ordinary differential equations for target gene concentrations that include Hill function dependences on regulator gene concentrations. An alternative formulation represents the same interactions using Boolean logic with time delays associated with each network link. We consider the attractors that emerge from the two types of models in the case of a simple but nontrivial network: a figure-8 network with one positive and one negative feedback loop. We show that the different modeling approaches give rise to the same qualitative set of attractors with the exception of a possible fixed point in the ordinary differential equation model in which concentrations sit at intermediate values. The properties of the attractors are most easily understood from the Boolean perspective, suggesting that time-delay Boolean modeling is a useful tool for understanding the logic of regulatory networks.
-
International Nuclear Information System (INIS)
Lima, M.L.; Mignaco, J.A.
1985-01-01
It is shown that the rational power law potentials in the two-body radial Schrodinger equations admit a systematic treatment available from the classical theory of ordinary linear differential equations of the second order. The resulting potentials come into families evolved from equations having a fixed number of elementary regular singularities. As a consequence, relations are found and discussed among the several potentials in a family. (Author) [pt
-
Robust estimation for ordinary differential equation models.
Cao, J; Wang, L; Xu, J
2011-12-01
Applied scientists often like to use ordinary differential equations (ODEs) to model complex dynamic processes that arise in biology, engineering, medicine, and many other areas. It is interesting but challenging to estimate ODE parameters from noisy data, especially when the data have some outliers. We propose a robust method to address this problem. The dynamic process is represented with a nonparametric function, which is a linear combination of basis functions. The nonparametric function is estimated by a robust penalized smoothing method. The penalty term is defined with the parametric ODE model, which controls the roughness of the nonparametric function and maintains the fidelity of the nonparametric function to the ODE model. The basis coefficients and ODE parameters are estimated in two nested levels of optimization. The coefficient estimates are treated as an implicit function of ODE parameters, which enables one to derive the analytic gradients for optimization using the implicit function theorem. Simulation studies show that the robust method gives satisfactory estimates for the ODE parameters from noisy data with outliers. The robust method is demonstrated by estimating a predator-prey ODE model from real ecological data. © 2011, The International Biometric Society.
-
Nonlinear ordinary differential equations analytical approximation and numerical methods
Hermann, Martin
2016-01-01
The book discusses the solutions to nonlinear ordinary differential equations (ODEs) using analytical and numerical approximation methods. Recently, analytical approximation methods have been largely used in solving linear and nonlinear lower-order ODEs. It also discusses using these methods to solve some strong nonlinear ODEs. There are two chapters devoted to solving nonlinear ODEs using numerical methods, as in practice high-dimensional systems of nonlinear ODEs that cannot be solved by analytical approximate methods are common. Moreover, it studies analytical and numerical techniques for the treatment of parameter-depending ODEs. The book explains various methods for solving nonlinear-oscillator and structural-system problems, including the energy balance method, harmonic balance method, amplitude frequency formulation, variational iteration method, homotopy perturbation method, iteration perturbation method, homotopy analysis method, simple and multiple shooting method, and the nonlinear stabilized march...
-
Camporesi, Roberto
2016-01-01
We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations of any order based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The approach presented here can be used in a first course on differential equations for science and engineering majors.
-
Cause and cure of sloppiness in ordinary differential equation models.
Tönsing, Christian; Timmer, Jens; Kreutz, Clemens
2014-08-01
Data-based mathematical modeling of biochemical reaction networks, e.g., by nonlinear ordinary differential equation (ODE) models, has been successfully applied. In this context, parameter estimation and uncertainty analysis is a major task in order to assess the quality of the description of the system by the model. Recently, a broadened eigenvalue spectrum of the Hessian matrix of the objective function covering orders of magnitudes was observed and has been termed as sloppiness. In this work, we investigate the origin of sloppiness from structures in the sensitivity matrix arising from the properties of the model topology and the experimental design. Furthermore, we present strategies using optimal experimental design methods in order to circumvent the sloppiness issue and present nonsloppy designs for a benchmark model.
-
Cause and cure of sloppiness in ordinary differential equation models
Tönsing, Christian; Timmer, Jens; Kreutz, Clemens
2014-08-01
Data-based mathematical modeling of biochemical reaction networks, e.g., by nonlinear ordinary differential equation (ODE) models, has been successfully applied. In this context, parameter estimation and uncertainty analysis is a major task in order to assess the quality of the description of the system by the model. Recently, a broadened eigenvalue spectrum of the Hessian matrix of the objective function covering orders of magnitudes was observed and has been termed as sloppiness. In this work, we investigate the origin of sloppiness from structures in the sensitivity matrix arising from the properties of the model topology and the experimental design. Furthermore, we present strategies using optimal experimental design methods in order to circumvent the sloppiness issue and present nonsloppy designs for a benchmark model.
-
From differential to difference equations for first order ODEs
Freed, Alan D.; Walker, Kevin P.
1991-01-01
When constructing an algorithm for the numerical integration of a differential equation, one should first convert the known ordinary differential equation (ODE) into an ordinary difference equation. Given this difference equation, one can develop an appropriate numerical algorithm. This technical note describes the derivation of two such ordinary difference equations applicable to a first order ODE. The implicit ordinary difference equation has the same asymptotic expansion as the ODE itself, whereas the explicit ordinary difference equation has an asymptotic that is similar in structure but different in value when compared with that of the ODE.
-
Barbu, Viorel
2016-01-01
This textbook is a comprehensive treatment of ordinary differential equations, concisely presenting basic and essential results in a rigorous manner. Including various examples from physics, mechanics, natural sciences, engineering and automatic theory, Differential Equations is a bridge between the abstract theory of differential equations and applied systems theory. Particular attention is given to the existence and uniqueness of the Cauchy problem, linear differential systems, stability theory and applications to first-order partial differential equations. Upper undergraduate students and researchers in applied mathematics and systems theory with a background in advanced calculus will find this book particularly useful. Supplementary topics are covered in an appendix enabling the book to be completely self-contained.
-
Camporesi, Roberto
2016-01-01
We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations of any order based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as…
-
International Nuclear Information System (INIS)
Datta, Dhurjati Prasad; Bose, Manoj Kumar
2004-01-01
We present a new one parameter family of second derivative discontinuous solutions to the simplest scale invariant linear ordinary differential equation. We also point out how the construction could be extended to generate families of higher derivative discontinuous solutions as well. The discontinuity can occur only for a subset of even order derivatives, viz., 2nd, 4th, 8th, 16th,.... The solutions are shown to break the discrete parity (reflection) symmetry of the underlying equation. These results are expected to gain significance in the contemporary search of a new dynamical principle for understanding complex phenomena in nature
-
Avellar, J.; Claudino, A. L. G. C.; Duarte, L. G. S.; da Mota, L. A. C. P.
2015-10-01
For the Darbouxian methods we are studying here, in order to solve first order rational ordinary differential equations (1ODEs), the most costly (computationally) step is the finding of the needed Darboux polynomials. This can be so grave that it can render the whole approach unpractical. Hereby we introduce a simple heuristics to speed up this process for a class of 1ODEs.
-
VCODE, Ordinary Differential Equation Solver for Stiff and Non-Stiff Problems
International Nuclear Information System (INIS)
Cohen, Scott D.; Hindmarsh, Alan C.
2001-01-01
1 - Description of program or function: CVODE is a package written in ANSI standard C for solving initial value problems for ordinary differential equations. It solves both stiff and non stiff systems. In the stiff case, it includes a variety of options for treating the Jacobian of the system, including dense and band matrix solvers, and a preconditioned Krylov (iterative) solver. 2 - Method of solution: Integration is by Adams or BDF (Backward Differentiation Formula) methods, at user option. Corrector iteration is by functional iteration or Newton iteration. For the solution of linear systems within Newton iteration, users can select a dense solver, a band solver, a diagonal approximation, or a preconditioned Generalized Minimal Residual (GMRES) solver. In the dense and band cases, the user can supply a Jacobian approximation or let CVODE generate it internally. In the GMRES case, the pre-conditioner is user-supplied
-
Directory of Open Access Journals (Sweden)
Hasan Bulut
2013-01-01
Full Text Available We introduce the rudiments of fractional calculus and the consequent applications of the Sumudu transform on fractional derivatives. Once this connection is firmly established in the general setting, we turn to the application of the Sumudu transform method (STM to some interesting nonhomogeneous fractional ordinary differential equations (FODEs. Finally, we use the solutions to form two-dimensional (2D graphs, by using the symbolic algebra package Mathematica Program 7.
-
Assessment of available integration algorithms for initial value ordinary differential equations
International Nuclear Information System (INIS)
Carver, M.B.; Stewart, D.G.
1979-11-01
There exists an extremely large number of algorithms designed for the ordinary differential equation initial value problem. The integration is normally done by a finite sum at time intervals which are chosen dynamically to satisfy an imposed error tolerance. This report describes the basic logistics of the integration process, identifies common areas of difficulty, and establishes a comprehensive test profile for integration algorithms. A number of algorithms are described, and selected published subroutines are evaluated using the test profile. It concludes that an effective library for general use need have only two such routines. The two selected are versions of the well-known Gear and Runge-Kutta-Fehlberg algorithms. Full documentation and listings are included. (auth)
-
Directory of Open Access Journals (Sweden)
Winfried Auzinger
2006-01-01
Full Text Available We demonstrate that eigenvalue problems for ordinary differential equations can be recast in a formulation suitable for the solution by polynomial collocation. It is shown that the well-posedness of the two formulations is equivalent in the regular as well as in the singular case. Thus, a collocation code equipped with asymptotically correct error estimation and adaptive mesh selection can be successfully applied to compute the eigenvalues and eigenfunctions efficiently and with reliable control of the accuracy. Numerical examples illustrate this claim.
-
On the relation between elementary partial difference equations and partial differential equations
van den Berg, I.P.
1998-01-01
The nonstandard stroboscopy method links discrete-time ordinary difference equations of first-order and continuous-time, ordinary differential equations of first order. We extend this method to the second order, and also to an elementary, yet general class of partial difference/differential
-
International Nuclear Information System (INIS)
Li Pengfei; Hu Gang; Chen Runsheng
2008-01-01
Gene transcriptional regulation (TR) processes are often described by coupled nonlinear ordinary differential equations (ODEs). When the dimension of TR circuits is high (e.g. n ≥ 3) the motions of the corresponding ODEs may, very probably, show self-sustained oscillations and chaos. On the other hand, chaoticity may be harmful for the normal biological functions of TR processes. In this letter we numerically study the dynamics of 3-gene TR ODEs in great detail, and investigate many 4-, 5-, and 10-gene TR systems by randomly choosing figures and parameters in the conventionally accepted ranges. And we find that oscillations are very seldom and no chaotic motion is observed, even if the dimension of systems is sufficiently high (n ≥ 3). It is argued that the observation of nonchaoticity of these ODEs agrees with normal functions of actual TR processes
-
Directory of Open Access Journals (Sweden)
M. P. Markakis
2010-01-01
Full Text Available Certain nonlinear autonomous ordinary differential equations of the second order are reduced to Abel equations of the first kind ((Ab-1 equations. Based on the results of a previous work, concerning a closed-form solution of a general (Ab-1 equation, and introducing an arbitrary function, exact one-parameter families of solutions are derived for the original autonomous equations, for the most of which only first integrals (in closed or parametric form have been obtained so far. Two-dimensional autonomous systems of differential equations of the first order, equivalent to the considered herein autonomous forms, are constructed and solved by means of the developed analysis.
-
Calculation of similarity solutions of partial differential equations
International Nuclear Information System (INIS)
Dresner, L.
1980-08-01
When a partial differential equation in two independent variables is invariant to a group G of stretching transformations, it has similarity solutions that can be found by solving an ordinary differential equation. Under broad conditions, this ordinary differential equation is also invariant to another stretching group G', related to G. The invariance of the ordinary differential equation to G' can be used to simplify its solution, particularly if it is of second order. Then a method of Lie's can be used to reduce it to a first-order equation, the study of which is greatly facilitated by analysis of its direction field. The method developed here is applied to three examples: Blasius's equation for boundary layer flow over a flat plate and two nonlinear diffusion equations, cc/sub t/ = c/sub zz/ and c/sub t/ = (cc/sub z/)/sub z/
-
Differential equations methods and applications
Said-Houari, Belkacem
2015-01-01
This book presents a variety of techniques for solving ordinary differential equations analytically and features a wealth of examples. Focusing on the modeling of real-world phenomena, it begins with a basic introduction to differential equations, followed by linear and nonlinear first order equations and a detailed treatment of the second order linear equations. After presenting solution methods for the Laplace transform and power series, it lastly presents systems of equations and offers an introduction to the stability theory. To help readers practice the theory covered, two types of exercises are provided: those that illustrate the general theory, and others designed to expand on the text material. Detailed solutions to all the exercises are included. The book is excellently suited for use as a textbook for an undergraduate class (of all disciplines) in ordinary differential equations. .
-
Lavrova, Anastasia I; Postnikov, Eugene B; Zyubin, Andrey Yu; Babak, Svetlana V
2017-04-01
We consider two approaches to modelling the cell metabolism of 6-mercaptopurine, one of the important chemotherapy drugs used for treating acute lymphocytic leukaemia: kinetic ordinary differential equations, and Boolean networks supplied with one controlling node, which takes continual values. We analyse their interplay with respect to taking into account ATP concentration as a key parameter of switching between different pathways. It is shown that the Boolean networks, which allow avoiding the complexity of general kinetic modelling, preserve the possibility of reproducing the principal switching mechanism.
-
International Nuclear Information System (INIS)
Pradeep, R Gladwin; Chandrasekar, V K; Senthilvelan, M; Lakshmanan, M
2011-01-01
In this paper, we devise a systematic procedure to obtain nonlocal symmetries of a class of scalar nonlinear ordinary differential equations (ODEs) of arbitrary order related to linear ODEs through nonlocal relations. The procedure makes use of the Lie point symmetries of the linear ODEs and the nonlocal connection to deduce the nonlocal symmetries of the corresponding nonlinear ODEs. Using these nonlocal symmetries, we obtain reduction transformations and reduced equations to specific examples. We find that the reduced equations can be explicitly integrated to deduce the general solutions for these cases. We also extend this procedure to coupled higher order nonlinear ODEs with specific reference to second-order nonlinear ODEs. (paper)
-
Sparse Additive Ordinary Differential Equations for Dynamic Gene Regulatory Network Modeling.
Wu, Hulin; Lu, Tao; Xue, Hongqi; Liang, Hua
2014-04-02
The gene regulation network (GRN) is a high-dimensional complex system, which can be represented by various mathematical or statistical models. The ordinary differential equation (ODE) model is one of the popular dynamic GRN models. High-dimensional linear ODE models have been proposed to identify GRNs, but with a limitation of the linear regulation effect assumption. In this article, we propose a sparse additive ODE (SA-ODE) model, coupled with ODE estimation methods and adaptive group LASSO techniques, to model dynamic GRNs that could flexibly deal with nonlinear regulation effects. The asymptotic properties of the proposed method are established and simulation studies are performed to validate the proposed approach. An application example for identifying the nonlinear dynamic GRN of T-cell activation is used to illustrate the usefulness of the proposed method.
-
Pure soliton solutions of some nonlinear partial differential equations
International Nuclear Information System (INIS)
Fuchssteiner, B.
1977-01-01
A general approach is given to obtain the system of ordinary differential equations which determines the pure soliton solutions for the class of generalized Korteweg-de Vries equations. This approach also leads to a system of ordinary differential equations for the pure soliton solutions of the sine-Gordon equation. (orig.) [de
-
Nahay, John Michael
2008-01-01
We present a new approach to solving polynomial ordinary differential equations by transforming them to linear functional equations and then solving the linear functional equations. We will focus most of our attention upon the first-order Abel differential equation with two nonlinear terms in order to demonstrate in as much detail as possible the computations necessary for a complete solution. We mention in our section on further developments that the basic transformation idea can be generali...
-
Directory of Open Access Journals (Sweden)
Mehmet Tarik Atay
2013-01-01
Full Text Available The Variational Iteration Method (VIM and Modified Variational Iteration Method (MVIM are used to find solutions of systems of stiff ordinary differential equations for both linear and nonlinear problems. Some examples are given to illustrate the accuracy and effectiveness of these methods. We compare our results with exact results. In some studies related to stiff ordinary differential equations, problems were solved by Adomian Decomposition Method and VIM and Homotopy Perturbation Method. Comparisons with exact solutions reveal that the Variational Iteration Method (VIM and the Modified Variational Iteration Method (MVIM are easier to implement. In fact, these methods are promising methods for various systems of linear and nonlinear stiff ordinary differential equations. Furthermore, VIM, or in some cases MVIM, is giving exact solutions in linear cases and very satisfactory solutions when compared to exact solutions for nonlinear cases depending on the stiffness ratio of the stiff system to be solved.
-
Directory of Open Access Journals (Sweden)
Decio Levi
2013-10-01
Full Text Available We briefly review two different methods of applying Lie group theory in the numerical solution of ordinary differential equations. On specific examples we show how the symmetry preserving discretization provides difference schemes for which the “first differential approximation” is invariant under the same Lie group as the original ordinary differential equation.
-
Elements of partial differential equations
Sneddon, Ian Naismith
1957-01-01
Geared toward students of applied rather than pure mathematics, this volume introduces elements of partial differential equations. Its focus is primarily upon finding solutions to particular equations rather than general theory.Topics include ordinary differential equations in more than two variables, partial differential equations of the first and second orders, Laplace's equation, the wave equation, and the diffusion equation. A helpful Appendix offers information on systems of surfaces, and solutions to the odd-numbered problems appear at the end of the book. Readers pursuing independent st
-
Non-instantaneous impulses in differential equations
Agarwal, Ravi; O'Regan, Donal
2017-01-01
This monograph is the first published book devoted to the theory of differential equations with non-instantaneous impulses. It aims to equip the reader with mathematical models and theory behind real life processes in physics, biology, population dynamics, ecology and pharmacokinetics. The authors examine a wide scope of differential equations with non-instantaneous impulses through three comprehensive chapters, providing an all-rounded and unique presentation on the topic, including: - Ordinary differential equations with non-instantaneous impulses (scalar and n-dimensional case) - Fractional differential equa tions with non-instantaneous impulses (with Caputo fractional derivatives of order q ϵ (0, 1)) - Ordinary differential equations with non-instantaneous impulses occurring at random moments (with exponential, Erlang, or Gamma distribution) Each chapter focuses on theory, proofs and examples, and contains numerous graphs to enrich the reader’s understanding. Additionally, a carefully selected bibliogr...
-
Chow, Sy-Miin; Bendezú, Jason J; Cole, Pamela M; Ram, Nilam
2016-01-01
Several approaches exist for estimating the derivatives of observed data for model exploration purposes, including functional data analysis (FDA; Ramsay & Silverman, 2005 ), generalized local linear approximation (GLLA; Boker, Deboeck, Edler, & Peel, 2010 ), and generalized orthogonal local derivative approximation (GOLD; Deboeck, 2010 ). These derivative estimation procedures can be used in a two-stage process to fit mixed effects ordinary differential equation (ODE) models. While the performance and utility of these routines for estimating linear ODEs have been established, they have not yet been evaluated in the context of nonlinear ODEs with mixed effects. We compared properties of the GLLA and GOLD to an FDA-based two-stage approach denoted herein as functional ordinary differential equation with mixed effects (FODEmixed) in a Monte Carlo (MC) study using a nonlinear coupled oscillators model with mixed effects. Simulation results showed that overall, the FODEmixed outperformed both the GLLA and GOLD across all the embedding dimensions considered, but a novel use of a fourth-order GLLA approach combined with very high embedding dimensions yielded estimation results that almost paralleled those from the FODEmixed. We discuss the strengths and limitations of each approach and demonstrate how output from each stage of FODEmixed may be used to inform empirical modeling of young children's self-regulation.
-
Auxiliary equation method for solving nonlinear partial differential equations
International Nuclear Information System (INIS)
Sirendaoreji,; Jiong, Sun
2003-01-01
By using the solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained with the aid of symbolic computation
-
Symmetries, Integrals and Solutions of Ordinary Differential ...
Indian Academy of Sciences (India)
Second-and third-order scalar ordinary differential equations of maximal symmetry in the traditional sense of point, respectively contact, symmetry are examined for the mappings they produce in solutions and fundamental first integrals. The properties of the `exceptional symmetries', i.e. those not considered to be generic to ...
-
Nonlinear differential equations
Energy Technology Data Exchange (ETDEWEB)
Dresner, L.
1988-01-01
This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics.
-
Nonlinear differential equations
International Nuclear Information System (INIS)
Dresner, L.
1988-01-01
This report is the text of a graduate course on nonlinear differential equations given by the author at the University of Wisconsin-Madison during the summer of 1987. The topics covered are: direction fields of first-order differential equations; the Lie (group) theory of ordinary differential equations; similarity solutions of second-order partial differential equations; maximum principles and differential inequalities; monotone operators and iteration; complementary variational principles; and stability of numerical methods. The report should be of interest to graduate students, faculty, and practicing scientists and engineers. No prior knowledge is required beyond a good working knowledge of the calculus. The emphasis is on practical results. Most of the illustrative examples are taken from the fields of nonlinear diffusion, heat and mass transfer, applied superconductivity, and helium cryogenics
-
Linares, Oscar A; Schiesser, William E; Fudin, Jeffrey; Pham, Thien C; Bettinger, Jeffrey J; Mathew, Roy O; Daly, Annemarie L
2015-01-01
Background There is a need to have a model to study methadone’s losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. Aim To build a one-dimensional (1-D), hollow-fiber geometry, ordinary differential equation (ODE) and partial differential equation (PDE) countercurrent hemodialyzer model (ODE/PDE model). Methodology We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone’s overall intradialytic mass transfer rate coefficient, km; and, methadone’s removal rate, jME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. Results The ODE/PDE model revealed a significant increase in the change of methadone’s mass transfer with increased dialysate flow rate, %Δkm=18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11). This was accompanied by a small significant increase in methadone’s mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11). The ODE/PDE model accurately predicted methadone’s removal during dialysis. The absolute value of the prediction errors for methadone’s extraction and throughput were less than 2%. Conclusion ODE/PDE modeling of methadone’s hemodialysis is a new approach to study methadone’s removal, in particular, and opioid removal, in general, in patients with end-stage renal disease on hemodialysis. ODE/PDE modeling accurately quantified the fundamental phenomena of methadone’s mass transfer during hemodialysis. This methodology may lead to development of optimally designed intradialytic opioid treatment protocols, and allow dynamic monitoring of outflow plasma opioid concentrations for model
-
Linares, Oscar A; Schiesser, William E; Fudin, Jeffrey; Pham, Thien C; Bettinger, Jeffrey J; Mathew, Roy O; Daly, Annemarie L
2015-01-01
There is a need to have a model to study methadone's losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. To build a one-dimensional (1-D), hollow-fiber geometry, ordinary differential equation (ODE) and partial differential equation (PDE) countercurrent hemodialyzer model (ODE/PDE model). We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone's overall intradialytic mass transfer rate coefficient, km ; and, methadone's removal rate, j ME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. The ODE/PDE model revealed a significant increase in the change of methadone's mass transfer with increased dialysate flow rate, %Δkm =18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11). This was accompanied by a small significant increase in methadone's mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11). The ODE/PDE model accurately predicted methadone's removal during dialysis. The absolute value of the prediction errors for methadone's extraction and throughput were less than 2%. ODE/PDE modeling of methadone's hemodialysis is a new approach to study methadone's removal, in particular, and opioid removal, in general, in patients with end-stage renal disease on hemodialysis. ODE/PDE modeling accurately quantified the fundamental phenomena of methadone's mass transfer during hemodialysis. This methodology may lead to development of optimally designed intradialytic opioid treatment protocols, and allow dynamic monitoring of outflow plasma opioid concentrations for model predictive control during dialysis in humans.
-
International Nuclear Information System (INIS)
Srivastava, A.C.; Hazarika, G.C.
1990-01-01
An algorithm based on the shooting method has been developed for the solution of a two-point boundary value problem consisting of a system of third order simultaneous ordinary differential equations. The Falkner-Skan equations for electrically conducting viscous fluid with applied magnetic field has been solved by using this algorithm for various values of the wedge angle and magnetic parameters. The shooting method seems to be well convergent for a system as the results are in good agreement with those obtained by other methods. It is observed that both viscous boundary layer and magnetic boundary layer decrease while velocity as well as magnetic field increase with the increase of the wedge angle. (author). 6 tabs., 7 refs
-
FORSIM-6, Automatic Solution of Coupled Differential Equation System
International Nuclear Information System (INIS)
Carver, M.B.; Stewart, D.G.; Blair, J.M.; Selander, W.N.
1983-01-01
1 - Description of problem or function: The FORSIM program is a versatile package which automates the solution of coupled differential equation systems. The independent variables are time, and up to three space coordinates, and the equations may be any mixture of partial and/or ordinary differential equations. The philosophy of the program is to provide a tool which will solve a system of differential equations for a user who has basic but unspecialized knowledge of numerical analysis and FORTRAN. The equations to be solved, together with the initial conditions and any special instructions, may be specified by the user in a single FORTRAN subroutine, although he may write a number of routines if this is more suitable. These are then loaded with the control routines, which perform the solution and any requested input and output. 2 - Method of solution: Partial differential equations are automatically converted into sets of coupled ordinary differential equations by variable order discretization in the spatial dimensions. These and other ordinary differential equations are integrated continuously in time using efficient variable order, variable step, error-controlled algorithms
-
Differential-algebraic solutions of the heat equation
Buchstaber, Victor M.; Netay, Elena Yu.
2014-01-01
In this work we introduce the notion of differential-algebraic ansatz for the heat equation and explicitly construct heat equation and Burgers equation solutions given a solution of a homogeneous non-linear ordinary differential equation of a special form. The ansatz for such solutions is called the $n$-ansatz, where $n+1$ is the order of the differential equation.
-
A perturbative solution to metadynamics ordinary differential equation
Tiwary, Pratyush; Dama, James F.; Parrinello, Michele
2015-12-01
Metadynamics is a popular enhanced sampling scheme wherein by periodic application of a repulsive bias, one can surmount high free energy barriers and explore complex landscapes. Recently, metadynamics was shown to be mathematically well founded, in the sense that the biasing procedure is guaranteed to converge to the true free energy surface in the long time limit irrespective of the precise choice of biasing parameters. A differential equation governing the post-transient convergence behavior of metadynamics was also derived. In this short communication, we revisit this differential equation, expressing it in a convenient and elegant Riccati-like form. A perturbative solution scheme is then developed for solving this differential equation, which is valid for any generic biasing kernel. The solution clearly demonstrates the robustness of metadynamics to choice of biasing parameters and gives further confidence in the widely used method.
-
A perturbative solution to metadynamics ordinary differential equation.
Tiwary, Pratyush; Dama, James F; Parrinello, Michele
2015-12-21
Metadynamics is a popular enhanced sampling scheme wherein by periodic application of a repulsive bias, one can surmount high free energy barriers and explore complex landscapes. Recently, metadynamics was shown to be mathematically well founded, in the sense that the biasing procedure is guaranteed to converge to the true free energy surface in the long time limit irrespective of the precise choice of biasing parameters. A differential equation governing the post-transient convergence behavior of metadynamics was also derived. In this short communication, we revisit this differential equation, expressing it in a convenient and elegant Riccati-like form. A perturbative solution scheme is then developed for solving this differential equation, which is valid for any generic biasing kernel. The solution clearly demonstrates the robustness of metadynamics to choice of biasing parameters and gives further confidence in the widely used method.
-
Tao, Youshan; Guo, Qian; Aihara, Kazuyuki
2014-10-01
Hormonal therapy with androgen suppression is a common treatment for advanced prostate tumors. The emergence of androgen-independent cells, however, leads to a tumor relapse under a condition of long-term androgen deprivation. Clinical trials suggest that intermittent androgen suppression (IAS) with alternating on- and off-treatment periods can delay the relapse when compared with continuous androgen suppression (CAS). In this paper, we propose a mathematical model for prostate tumor growth under IAS therapy. The model elucidates initial hormone sensitivity, an eventual relapse of a tumor under CAS therapy, and a delay of a relapse under IAS therapy, which are due to the coexistence of androgen-dependent cells, androgen-independent cells resulting from reversible changes by adaptation, and androgen-independent cells resulting from irreversible changes by genetic mutations. The model is formulated as a free boundary problem of partial differential equations that describe the evolution of populations of the abovementioned three types of cells during on-treatment periods and off-treatment periods. Moreover, the model can be transformed into a piecewise linear ordinary differential equation model by introducing three new volume variables, and the study of the resulting model may help to devise optimal IAS schedules.
-
SIVA/DIVA- INITIAL VALUE ORDINARY DIFFERENTIAL EQUATION SOLUTION VIA A VARIABLE ORDER ADAMS METHOD
Krogh, F. T.
1994-01-01
The SIVA/DIVA package is a collection of subroutines for the solution of ordinary differential equations. There are versions for single precision and double precision arithmetic. These solutions are applicable to stiff or nonstiff differential equations of first or second order. SIVA/DIVA requires fewer evaluations of derivatives than other variable order Adams predictor-corrector methods. There is an option for the direct integration of second order equations which can make integration of trajectory problems significantly more efficient. Other capabilities of SIVA/DIVA include: monitoring a user supplied function which can be separate from the derivative; dynamically controlling the step size; displaying or not displaying output at initial, final, and step size change points; saving the estimated local error; and reverse communication where subroutines return to the user for output or computation of derivatives instead of automatically performing calculations. The user must supply SIVA/DIVA with: 1) the number of equations; 2) initial values for the dependent and independent variables, integration stepsize, error tolerance, etc.; and 3) the driver program and operational parameters necessary for subroutine execution. SIVA/DIVA contains an extensive diagnostic message library should errors occur during execution. SIVA/DIVA is written in FORTRAN 77 for batch execution and is machine independent. It has a central memory requirement of approximately 120K of 8 bit bytes. This program was developed in 1983 and last updated in 1987.
-
Iterative Splitting Methods for Differential Equations
Geiser, Juergen
2011-01-01
Iterative Splitting Methods for Differential Equations explains how to solve evolution equations via novel iterative-based splitting methods that efficiently use computational and memory resources. It focuses on systems of parabolic and hyperbolic equations, including convection-diffusion-reaction equations, heat equations, and wave equations. In the theoretical part of the book, the author discusses the main theorems and results of the stability and consistency analysis for ordinary differential equations. He then presents extensions of the iterative splitting methods to partial differential
-
Fast integration-based prediction bands for ordinary differential equation models.
Hass, Helge; Kreutz, Clemens; Timmer, Jens; Kaschek, Daniel
2016-04-15
To gain a deeper understanding of biological processes and their relevance in disease, mathematical models are built upon experimental data. Uncertainty in the data leads to uncertainties of the model's parameters and in turn to uncertainties of predictions. Mechanistic dynamic models of biochemical networks are frequently based on nonlinear differential equation systems and feature a large number of parameters, sparse observations of the model components and lack of information in the available data. Due to the curse of dimensionality, classical and sampling approaches propagating parameter uncertainties to predictions are hardly feasible and insufficient. However, for experimental design and to discriminate between competing models, prediction and confidence bands are essential. To circumvent the hurdles of the former methods, an approach to calculate a profile likelihood on arbitrary observations for a specific time point has been introduced, which provides accurate confidence and prediction intervals for nonlinear models and is computationally feasible for high-dimensional models. In this article, reliable and smooth point-wise prediction and confidence bands to assess the model's uncertainty on the whole time-course are achieved via explicit integration with elaborate correction mechanisms. The corresponding system of ordinary differential equations is derived and tested on three established models for cellular signalling. An efficiency analysis is performed to illustrate the computational benefit compared with repeated profile likelihood calculations at multiple time points. The integration framework and the examples used in this article are provided with the software package Data2Dynamics, which is based on MATLAB and freely available at http://www.data2dynamics.org helge.hass@fdm.uni-freiburg.de Supplementary data are available at Bioinformatics online. © The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e
-
Particle Swarm Optimization Based on Local Attractors of Ordinary Differential Equation System
Directory of Open Access Journals (Sweden)
Wenyu Yang
2014-01-01
Full Text Available Particle swarm optimization (PSO is inspired by sociological behavior. In this paper, we interpret PSO as a finite difference scheme for solving a system of stochastic ordinary differential equations (SODE. In this framework, the position points of the swarm converge to an equilibrium point of the SODE and the local attractors, which are easily defined by the present position points, also converge to the global attractor. Inspired by this observation, we propose a class of modified PSO iteration methods (MPSO based on local attractors of the SODE. The idea of MPSO is to choose the next update state near the present local attractor, rather than the present position point as in the original PSO, according to a given probability density function. In particular, the quantum-behaved particle swarm optimization method turns out to be a special case of MPSO by taking a special probability density function. The MPSO methods with six different probability density functions are tested on a few benchmark problems. These MPSO methods behave differently for different problems. Thus, our framework not only gives an interpretation for the ordinary PSO but also, more importantly, provides a warehouse of PSO-like methods to choose from for solving different practical problems.
-
Mathematical modelling of tissue formation on the basis of ordinary differential equations
Directory of Open Access Journals (Sweden)
Maxim N. Nazarov
2017-10-01
Full Text Available A mathematical model is proposed for describing the population dynamics of cellular clusters on the basis of systems of the first-order ordinary differential equations. The main requirement for the construction of model equations was to obtain a formal biological justification for their derivation, as well as proof of their correctness. In addition, for all the parameters involved in the model equations, the presence of biological meaning was guaranteed, as well as the possibility of evaluating them either during the experiment or by using models of intracellular biochemistry. In the desired model the intercellular exchange of a special signal molecules was chosen as the main mechanism for coordination of the tissue growth and new types selection during cell division. For simplicity, all signalling molecules that can create cells of the same type were not considered separately in the model, but were instead combined in a single complex of molecules: a ‘generalized signal’. Such an approach allows us to eventually assign signals as a functions of cell types and introduce their effects in the form of matrices in the models, where the rows are responsible for the types of cells receiving the signals, and the columns for the types of cells emitting signals.
-
International Nuclear Information System (INIS)
Diaz Sanchidrian, C.
1989-01-01
The present paper applies dimensional analysis with spatial discrimination to transform the differential equations in partial derivatives developed in the theory of heat transmission into ordinary ones. The effectivity of the method is comparable to that methods based in transformations of uni or multiparametric groups, with the advantage of being more direct and simple. (Author)
-
A first course in ordinary differential equations analytical and numerical methods
Hermann, Martin
2014-01-01
This book presents a modern introduction to analytical and numerical techniques for solving ordinary differential equations (ODEs). Contrary to the traditional format—the theorem-and-proof format—the book is focusing on analytical and numerical methods. The book supplies a variety of problems and examples, ranging from the elementary to the advanced level, to introduce and study the mathematics of ODEs. The analytical part of the book deals with solution techniques for scalar first-order and second-order linear ODEs, and systems of linear ODEs—with a special focus on the Laplace transform, operator techniques and power series solutions. In the numerical part, theoretical and practical aspects of Runge-Kutta methods for solving initial-value problems and shooting methods for linear two-point boundary-value problems are considered. The book is intended as a primary text for courses on the theory of ODEs and numerical treatment of ODEs for advanced undergraduate and early graduate students. It is assumed t...
-
Energy Technology Data Exchange (ETDEWEB)
Moryakov, A. V., E-mail: sailor@orc.ru [National Research Centre Kurchatov Institute (Russian Federation)
2016-12-15
An algorithm for solving the linear Cauchy problem for large systems of ordinary differential equations is presented. The algorithm for systems of first-order differential equations is implemented in the EDELWEISS code with the possibility of parallel computations on supercomputers employing the MPI (Message Passing Interface) standard for the data exchange between parallel processes. The solution is represented by a series of orthogonal polynomials on the interval [0, 1]. The algorithm is characterized by simplicity and the possibility to solve nonlinear problems with a correction of the operator in accordance with the solution obtained in the previous iterative process.
-
Solving Differential Equations in R: Package deSolve
In this paper we present the R package deSolve to solve initial value problems (IVP) written as ordinary differential equations (ODE), differential algebraic equations (DAE) of index 0 or 1 and partial differential equations (PDE), the latter solved using the method of lines appr...
-
Introduction to partial differential equations with applications
Zachmanoglou, E C
1988-01-01
This text explores the essentials of partial differential equations as applied to engineering and the physical sciences. Discusses ordinary differential equations, integral curves and surfaces of vector fields, the Cauchy-Kovalevsky theory, more. Problems and answers.
-
Numerical Integration of Stiff System of Ordinary Differential ...
African Journals Online (AJOL)
The goal of this work is to develop, analyse and implement a K-step Implicit Rational Runge-Kutta schemes for Integration of Stiff system of Ordinary differential Equations. Its development adopted Taylor and Binomial series expansion Techniques to generate its parameters. The analysis of its basic properties adopted ...
-
DEFF Research Database (Denmark)
Mikkelsen, Frederik Vissing
eective computational tools for estimating unknown structures in dynamical systems, such as gene regulatory networks, which may be used to predict downstream eects of interventions in the system. A recommended algorithm based on the computational tools is presented and thoroughly tested in various......Broadly speaking, this thesis is devoted to model selection applied to ordinary dierential equations and risk estimation under model selection. A model selection framework was developed for modelling time course data by ordinary dierential equations. The framework is accompanied by the R software...... package, episode. This package incorporates a collection of sparsity inducing penalties into two types of loss functions: a squared loss function relying on numerically solving the equations and an approximate loss function based on inverse collocation methods. The goal of this framework is to provide...
-
Partial differential equations for scientists and engineers
Farlow, Stanley J
1993-01-01
Most physical phenomena, whether in the domain of fluid dynamics, electricity, magnetism, mechanics, optics, or heat flow, can be described in general by partial differential equations. Indeed, such equations are crucial to mathematical physics. Although simplifications can be made that reduce these equations to ordinary differential equations, nevertheless the complete description of physical systems resides in the general area of partial differential equations.This highly useful text shows the reader how to formulate a partial differential equation from the physical problem (constructing th
-
Solving Differential Equations in R: Package deSolve
Soetaert, K.E.R.; Petzoldt, T.; Setzer, R.W.
2010-01-01
In this paper we present the R package deSolve to solve initial value problems (IVP) written as ordinary differential equations (ODE), differential algebraic equations (DAE) of index 0 or 1 and partial differential equations (PDE), the latter solved using the method of lines approach. The
-
Tang, Kwong-Tin
2007-01-01
Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student-oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.
-
Lie symmetries in differential equations
International Nuclear Information System (INIS)
Pleitez, V.
1979-01-01
A study of ordinary and Partial Differential equations using the symmetries of Lie groups is made. Following such a study, an application to the Helmholtz, Line-Gordon, Korleweg-de Vries, Burguer, Benjamin-Bona-Mahony and wave equations is carried out [pt
-
LSODE, 1. Order Stiff or Non-Stiff Ordinary Differential Equations System Initial Value Problems
International Nuclear Information System (INIS)
Hindmarsh, A.C.; Petzold, L.R.
2005-01-01
1 - Description of program or function: LSODE (Livermore Solver for Ordinary Differential Equations) solves stiff and non-stiff systems of the form dy/dt = f. In the stiff case, it treats the Jacobian matrix df/dy as either a dense (full) or a banded matrix, and as either user-supplied or internally approximated by difference quotients. It uses Adams methods (predictor-corrector) in the non-stiff case, and Backward Differentiation Formula (BDF) methods (the Gear methods) in the stiff case. The linear systems that arise are solved by direct methods (LU factor/solve). The LSODE source is commented extensively to facilitate modification. Both a single-precision version and a double-precision version are available. 2 - Methods: It is assumed that the ODEs are given explicitly, so that the system can be written in the form dy/dt = f(t,y), where y is the vector of dependent variables, and t is the independent variable. LSODE contains two variable-order, variable- step (with interpolatory step-changing) integration methods. The first is the implicit Adams or non-stiff method, of orders one through twelve. The second is the backward differentiation or stiff method (or BDF method, or Gear's method), of orders one through five. 3 - Restrictions on the complexity of the problem: The differential equations must be given in explicit form, i.e., dy/dt = f(y,t). Problems with intermittent high-speed transients may cause inefficient or unstable performance
-
Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation
Directory of Open Access Journals (Sweden)
Hamidreza Rezazadeh
2014-05-01
Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.
-
Fuchs indices and the first integrals of nonlinear differential equations
International Nuclear Information System (INIS)
Kudryashov, Nikolai A.
2005-01-01
New method of finding the first integrals of nonlinear differential equations in polynomial form is presented. Basic idea of our approach is to use the scaling of solution of nonlinear differential equation and to find the dimensions of arbitrary constants in the Laurent expansion of the general solution. These dimensions allows us to obtain the scalings of members for the first integrals of nonlinear differential equations. Taking the polynomials with unknown coefficients into account we present the algorithm of finding the first integrals of nonlinear differential equations in the polynomial form. Our method is applied to look for the first integrals of eight nonlinear ordinary differential equations of the fourth order. The general solution of one of the fourth order ordinary differential equations is given
-
Generalized ordinary differential equations not absolutely continuous solutions
Kurzweil, Jaroslav
2012-01-01
This book provides a systematic treatment of the Volterra integral equation by means of a modern integration theory which extends considerably the field of differential equations. It contains many new concepts and results in the framework of a unifying theory. In particular, this new approach is suitable in situations where fast oscillations occur.
-
Solution of differential equations by application of transformation groups
Driskell, C. N., Jr.; Gallaher, L. J.; Martin, R. H., Jr.
1968-01-01
Report applies transformation groups to the solution of systems of ordinary differential equations and partial differential equations. Lies theorem finds an integrating factor for appropriate invariance group or groups can be found and can be extended to partial differential equations.
-
On a complex differential Riccati equation
International Nuclear Information System (INIS)
Khmelnytskaya, Kira V; Kravchenko, Vladislav V
2008-01-01
We consider a nonlinear partial differential equation for complex-valued functions which is related to the two-dimensional stationary Schroedinger equation and enjoys many properties similar to those of the ordinary differential Riccati equation such as the famous Euler theorems, the Picard theorem and others. Besides these generalizations of the classical 'one-dimensional' results, we discuss new features of the considered equation including an analogue of the Cauchy integral theorem
-
Differential equations from the algebraic standpoint
Ritt, Joseph Fels
1932-01-01
This book can be viewed as a first attempt to systematically develop an algebraic theory of nonlinear differential equations, both ordinary and partial. The main goal of the author was to construct a theory of elimination, which "will reduce the existence problem for a finite or infinite system of algebraic differential equations to the application of the implicit function theorem taken with Cauchy's theorem in the ordinary case and Riquier's in the partial." In his 1934 review of the book, J. M. Thomas called it "concise, readable, original, precise, and stimulating", and his words still rema
-
Solving Nonlinear Coupled Differential Equations
Mitchell, L.; David, J.
1986-01-01
Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.
-
The mixed BVP for second order nonlinear ordinary differential equation at resonance
Czech Academy of Sciences Publication Activity Database
Mukhigulashvili, Sulkhan
2017-01-01
Roč. 290, 2-3 (2017), s. 393-400 ISSN 0025-584X Institutional support: RVO:67985840 Keywords : mixed problem at resonance * nonlinear ordinary differencial equation Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.742, year: 2016
-
Linear measure functional differential equations with infinite delay
Monteiro, G. (Giselle Antunes); Slavík, A.
2014-01-01
We use the theory of generalized linear ordinary differential equations in Banach spaces to study linear measure functional differential equations with infinite delay. We obtain new results concerning the existence, uniqueness, and continuous dependence of solutions. Even for equations with a finite delay, our results are stronger than the existing ones. Finally, we present an application to functional differential equations with impulses.
-
ODEion--a software module for structural identification of ordinary differential equations.
Gennemark, Peter; Wedelin, Dag
2014-02-01
In the systems biology field, algorithms for structural identification of ordinary differential equations (ODEs) have mainly focused on fixed model spaces like S-systems and/or on methods that require sufficiently good data so that derivatives can be accurately estimated. There is therefore a lack of methods and software that can handle more general models and realistic data. We present ODEion, a software module for structural identification of ODEs. Main characteristic features of the software are: • The model space is defined by arbitrary user-defined functions that can be nonlinear in both variables and parameters, such as for example chemical rate reactions. • ODEion implements computationally efficient algorithms that have been shown to efficiently handle sparse and noisy data. It can run a range of realistic problems that previously required a supercomputer. • ODEion is easy to use and provides SBML output. We describe the mathematical problem, the ODEion system itself, and provide several examples of how the system can be used. Available at: http://www.odeidentification.org.
-
Hasegawa, Chihiro; Duffull, Stephen B
2018-02-01
Pharmacokinetic-pharmacodynamic systems are often expressed with nonlinear ordinary differential equations (ODEs). While there are numerous methods to solve such ODEs these methods generally rely on time-stepping solutions (e.g. Runge-Kutta) which need to be matched to the characteristics of the problem at hand. The primary aim of this study was to explore the performance of an inductive approximation which iteratively converts nonlinear ODEs to linear time-varying systems which can then be solved algebraically or numerically. The inductive approximation is applied to three examples, a simple nonlinear pharmacokinetic model with Michaelis-Menten elimination (E1), an integrated glucose-insulin model and an HIV viral load model with recursive feedback systems (E2 and E3, respectively). The secondary aim of this study was to explore the potential advantages of analytically solving linearized ODEs with two examples, again E3 with stiff differential equations and a turnover model of luteinizing hormone with a surge function (E4). The inductive linearization coupled with a matrix exponential solution provided accurate predictions for all examples with comparable solution time to the matched time-stepping solutions for nonlinear ODEs. The time-stepping solutions however did not perform well for E4, particularly when the surge was approximated by a square wave. In circumstances when either a linear ODE is particularly desirable or the uncertainty in matching the integrator to the ODE system is of potential risk, then the inductive approximation method coupled with an analytical integration method would be an appropriate alternative.
-
Numerical solution of ordinary differential equations. For classical, relativistic and nano systems
International Nuclear Information System (INIS)
Greenspan, D.
2006-01-01
An up-to-date survey on numerical solutions with theory, intuition and applications. Ordinary differential equations (ODE) play a significant role in mathematics, physics and engineering sciences, and thus are part of relevant college and university courses. Many problems, however, both traditional and modern, do not possess exact solutions, and must be treated numerically. Usually this is done with software packages, but for this to be efficient requires a sound understanding of the mathematics involved. This work meets the need for an affordable textbook that helps in understanding numerical solutions of ODE. Carefully structured by an experienced textbook author, it provides a survey of ODE for various applications, both classical and modern, including such special applications as relativistic and nano systems. The examples are carefully explained and compiled into an algorithm, each of which is presented generically, independent of a specific programming language, while each chapter is rounded off with exercises. The text meets the demands of MA200 courses and of the newly created Numerical Solution of Differential Equations courses, making it ideal for both students and lecturers in physics, mathematics, mechanical engineering, electrical engineering, as well as for physicists, mathematicians, engineers, and electrical engineers. From the Contents - Euler's Method - Runge-Kutta Methods - The Method of Taylor Expansions - Large Second Order Systems with Application to Nano Systems - Completely Conservative, Covariant Numerical Methodology - Instability - Numerical Solution of Tridiagonal Linear Algebraic Systems and Related Nonlinear Systems - Approximate Solution of Boundary Value Problems - Special Relativistic Motion - Special Topics - Appendix: Basic Matrix Operations - Bibliography. (orig.) (orig.)
-
On a quaternionic generalisation of the Riccati differential equation
Kravchenko, Viktor; Kravchenko, Vladislav; Williams, Benjamin
2001-01-01
A quaternionic partial differential equation is shown to be a generalisation of the Riccati ordinary differential equation and its relationship with the Schrodinger equation is established. Various approaches to the problem of finding particular solutions are explored, and the generalisations of two theorems of Euler on the Riccati differential equation, which correspond to the quaternionic equation, are given.
-
Legendre-tau approximations for functional differential equations
Ito, K.; Teglas, R.
1986-01-01
The numerical approximation of solutions to linear retarded functional differential equations are considered using the so-called Legendre-tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time-differentiation. The approximate solution is then represented as a truncated Legendre series with time-varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximation is made.
-
Ding, Xiao-Li; Nieto, Juan J.
2017-11-01
In this paper, we consider the analytical solutions of coupling fractional partial differential equations (FPDEs) with Dirichlet boundary conditions on a finite domain. Firstly, the method of successive approximations is used to obtain the analytical solutions of coupling multi-term time fractional ordinary differential equations. Then, the technique of spectral representation of the fractional Laplacian operator is used to convert the coupling FPDEs to the coupling multi-term time fractional ordinary differential equations. By applying the obtained analytical solutions to the resulting multi-term time fractional ordinary differential equations, the desired analytical solutions of the coupling FPDEs are given. Our results are applied to derive the analytical solutions of some special cases to demonstrate their applicability.
-
Differential equations inverse and direct problems
Favini, Angelo
2006-01-01
DEGENERATE FIRST ORDER IDENTIFICATION PROBLEMS IN BANACH SPACES A NONISOTHERMAL DYNAMICAL GINZBURG-LANDAU MODEL OF SUPERCONDUCTIVITY. EXISTENCE AND UNIQUENESS THEOREMSSOME GLOBAL IN TIME RESULTS FOR INTEGRODIFFERENTIAL PARABOLIC INVERSE PROBLEMSFOURTH ORDER ORDINARY DIFFERENTIAL OPERATORS WITH GENERAL WENTZELL BOUNDARY CONDITIONSTUDY OF ELLIPTIC DIFFERENTIAL EQUATIONS IN UMD SPACESDEGENERATE INTEGRODIFFERENTIAL EQUATIONS OF PARABOLIC TYPE EXPONENTIAL ATTRACTORS FOR SEMICONDUCTOR EQUATIONSCONVERGENCE TO STATIONARY STATES OF SOLUTIONS TO THE SEMILINEAR EQUATION OF VISCOELASTICITY ASYMPTOTIC BEHA
-
Complex centers of polynomial differential equations
Directory of Open Access Journals (Sweden)
Mohamad Ali M. Alwash
2007-07-01
Full Text Available We present some results on the existence and nonexistence of centers for polynomial first order ordinary differential equations with complex coefficients. In particular, we show that binomial differential equations without linear terms do not have complex centers. Classes of polynomial differential equations, with more than two terms, are presented that do not have complex centers. We also study the relation between complex centers and the Pugh problem. An algorithm is described to solve the Pugh problem for equations without complex centers. The method of proof involves phase plane analysis of the polar equations and a local study of periodic solutions.
-
Phase integral approximation for coupled ordinary differential equations of the Schroedinger type
International Nuclear Information System (INIS)
Skorupski, Andrzej A.
2008-01-01
Four generalizations of the phase integral approximation (PIA) to sets of ordinary differential equations of Schroedinger type [u j '' (x)+Σ k=1 N R jk (x)u k (x)=0, j=1,2,...,N] are described. The recurrence relations for higher order corrections are given in a form valid to arbitrary order and for the matrix R(x)[≡(R jk (x))] either Hermitian or non-Hermitian. For Hermitian and negative definite R(x) matrices, a Wronskian conserving PIA theory is formulated, which generalizes Fulling's current conserving theory pertinent to positive definite R(x) matrices. The idea of a modification of the PIA, which is well known for one equation [u '' (x)+R(x)u(x)=0], is generalized to sets. A simplification of Wronskian or current conserving theories is proposed which in each order eliminates one integration from the formulas for higher order corrections. If the PIA is generated by a nondegenerate eigenvalue of the R(x) matrix, the eliminated integration is the only one present. In that case, the simplified theory becomes fully algorithmic and is generalized to non-Hermitian R(x) matrices. The general theory is illustrated by a few examples automatically generated by using the author's program in MATHEMATICA published in e-print arXiv:0710.5406 [math-ph
-
Variable-mesh method of solving differential equations
Van Wyk, R.
1969-01-01
Multistep predictor-corrector method for numerical solution of ordinary differential equations retains high local accuracy and convergence properties. In addition, the method was developed in a form conducive to the generation of effective criteria for the selection of subsequent step sizes in step-by-step solution of differential equations.
-
Differential equations, mechanics, and computation
Palais, Richard S
2009-01-01
This book provides a conceptual introduction to the theory of ordinary differential equations, concentrating on the initial value problem for equations of evolution and with applications to the calculus of variations and classical mechanics, along with a discussion of chaos theory and ecological models. It has a unified and visual introduction to the theory of numerical methods and a novel approach to the analysis of errors and stability of various numerical solution algorithms based on carefully chosen model problems. While the book would be suitable as a textbook for an undergraduate or elementary graduate course in ordinary differential equations, the authors have designed the text also to be useful for motivated students wishing to learn the material on their own or desiring to supplement an ODE textbook being used in a course they are taking with a text offering a more conceptual approach to the subject.
-
Bartels, Robert E.
2002-01-01
A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.
-
A Line-Tau Collocation Method for Partial Differential Equations ...
African Journals Online (AJOL)
This paper deals with the numerical solution of second order linear partial differential equations with the use of the method of lines coupled with the tau collocation method. The method of lines is used to convert the partial differential equation (PDE) to a sequence of ordinary differential equations (ODEs) which is then ...
-
Workshop on Numerical Methods for Ordinary Differential Equations
Gear, Charles; Russo, Elvira
1989-01-01
Developments in numerical initial value ode methods were the focal topic of the meeting at L'Aquila which explord the connections between the classical background and new research areas such as differental-algebraic equations, delay integral and integro-differential equations, stability properties, continuous extensions (interpolants for Runge-Kutta methods and their applications, effective stepsize control, parallel algorithms for small- and large-scale parallel architectures). The resulting proceedings address many of these topics in both research and survey papers.
-
Exponentially Convergent Algorithms for Abstract Differential Equations
Gavrilyuk, Ivan; Vasylyk, Vitalii
2011-01-01
This book presents new accurate and efficient exponentially convergent methods for abstract differential equations with unbounded operator coefficients in Banach space. These methods are highly relevant for the practical scientific computing since the equations under consideration can be seen as the meta-models of systems of ordinary differential equations (ODE) as well as the partial differential equations (PDEs) describing various applied problems. The framework of functional analysis allows one to obtain very general but at the same time transparent algorithms and mathematical results which
-
Introduction to numerical methods for time dependent differential equations
Kreiss, Heinz-Otto
2014-01-01
Introduces both the fundamentals of time dependent differential equations and their numerical solutions Introduction to Numerical Methods for Time Dependent Differential Equations delves into the underlying mathematical theory needed to solve time dependent differential equations numerically. Written as a self-contained introduction, the book is divided into two parts to emphasize both ordinary differential equations (ODEs) and partial differential equations (PDEs). Beginning with ODEs and their approximations, the authors provide a crucial presentation of fundamental notions, such as the t
-
Savoye, Philippe
2009-01-01
In recent years, I started covering difference equations and z transform methods in my introductory differential equations course. This allowed my students to extend the "classical" methods for (ordinary differential equation) ODE's to discrete time problems arising in many applications.
-
The Numerical Solution of an Abelian Ordinary Differential Equation ...
African Journals Online (AJOL)
In this paper we present a relatively new technique call theNew Hybrid of Adomian decomposition method (ADM) for solution of an Abelian Differential equation. The numerical results of the equation have been obtained in terms of convergent series with easily computable component. These methods are applied to solve ...
-
Directory of Open Access Journals (Sweden)
Linares OA
2015-07-01
Full Text Available Oscar A Linares,1 William E Schiesser,2 Jeffrey Fudin,3–6 Thien C Pham,6 Jeffrey J Bettinger,6 Roy O Mathew,6 Annemarie L Daly7 1Translational Genomic Medicine Lab, Plymouth Pharmacokinetic Modeling Study Group, Plymouth, MI, 2Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, PA, 3University of Connecticut School of Pharmacy, Storrs, CT, 4Western New England College of Pharmacy, Springfield, MA, 5Albany College of Pharmacy and Health Sciences, Albany, NY, 6Stratton VA Medical Center, Albany, NY, 7Grace Hospice of Ann Arbor, Ann Arbor, MI, USA Background: There is a need to have a model to study methadone’s losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. Aim: To build a one-dimensional (1-D, hollow-fiber geometry, ordinary differential equation (ODE and partial differential equation (PDE countercurrent hemodialyzer model (ODE/PDE model. Methodology: We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone’s overall intradialytic mass transfer rate coefficient, km; and, methadone’s removal rate, jME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. Results: The ODE/PDE model revealed a significant increase in the change of methadone’s mass transfer with increased dialysate flow rate, %Δ km=18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11. This was accompanied by a small significant increase in methadone’s mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11. The ODE/PDE model accurately predicted methadone’s removal during dialysis. The absolute value
-
Goldston, J. W.
This unit introduces analytic solutions of ordinary differential equations. The objective is to enable the student to decide whether a given function solves a given differential equation. Examples of problems from biology and chemistry are covered. Problem sets, quizzes, and a model exam are included, and answers to all items are provided. The…
-
Cao, Jiguo; Huang, Jianhua Z.; Wu, Hulin
2012-01-01
Ordinary differential equations (ODEs) are widely used in biomedical research and other scientific areas to model complex dynamic systems. It is an important statistical problem to estimate parameters in ODEs from noisy observations. In this article we propose a method for estimating the time-varying coefficients in an ODE. Our method is a variation of the nonlinear least squares where penalized splines are used to model the functional parameters and the ODE solutions are approximated also using splines. We resort to the implicit function theorem to deal with the nonlinear least squares objective function that is only defined implicitly. The proposed penalized nonlinear least squares method is applied to estimate a HIV dynamic model from a real dataset. Monte Carlo simulations show that the new method can provide much more accurate estimates of functional parameters than the existing two-step local polynomial method which relies on estimation of the derivatives of the state function. Supplemental materials for the article are available online.
-
Stability analysis of impulsive functional differential equations
Stamova, Ivanka
2009-01-01
This book is devoted to impulsive functional differential equations which are a natural generalization of impulsive ordinary differential equations (without delay) and of functional differential equations (without impulses). At the present time the qualitative theory of such equationsis under rapid development. After a presentation of the fundamental theory of existence, uniqueness and continuability of solutions, a systematic development of stability theory for that class of problems is given which makes the book unique. It addresses to a wide audience such as mathematicians, applied research
-
Linear measure functional differential equations with infinite delay
Czech Academy of Sciences Publication Activity Database
Monteiro, Giselle Antunes; Slavík, A.
2014-01-01
Roč. 287, 11-12 (2014), s. 1363-1382 ISSN 0025-584X Institutional support: RVO:67985840 Keywords : measure functional differential equations * generalized ordinary differential equations * Kurzweil-Stieltjes integral Subject RIV: BA - General Mathematics Impact factor: 0.683, year: 2014 http://onlinelibrary.wiley.com/doi/10.1002/mana.201300048/abstract
-
Selected papers on analysis and differential equations
Society, American Mathematical
2010-01-01
This volume contains translations of papers that originally appeared in the Japanese journal Sūgaku. These papers range over a variety of topics in ordinary and partial differential equations, and in analysis. Many of them are survey papers presenting new results obtained in the last few years. This volume is suitable for graduate students and research mathematicians interested in analysis and differential equations.
-
Rosenblatt, Marcus; Timmer, Jens; Kaschek, Daniel
2016-01-01
Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g., by maximum-likelihood estimation. In particular, models of biological systems contain a large number of parameters. To reduce the dimensionality of the parameter space, steady-state information is incorporated in the parameter estimation process. For non-linear models, analytical steady-state calculation typically leads to higher-order polynomial equations for which no closed-form solutions can be obtained. This can be circumvented by solving the steady-state equations for kinetic parameters, which results in a linear equation system with comparatively simple solutions. At the same time multiplicity of steady-state solutions is avoided, which otherwise is problematic for optimization. When solved for kinetic parameters, however, steady-state constraints tend to become negative for particular model specifications, thus, generating new types of optimization problems. Here, we present an algorithm based on graph theory that derives non-negative, analytical steady-state expressions by stepwise removal of cyclic dependencies between dynamical variables. The algorithm avoids multiple steady-state solutions by construction. We show that our method is applicable to most common classes of biochemical reaction networks containing inhibition terms, mass-action and Hill-type kinetic equations. Comparing the performance of parameter estimation for different analytical and numerical methods of incorporating steady-state information, we show that our approach is especially well-tailored to guarantee a high success rate of optimization.
-
Dorren, H.J.S.
1998-01-01
It is shown that the Korteweg–de Vries (KdV) equation can be transformed into an ordinary linear partial differential equation in the wave number domain. Explicit solutions of the KdV equation can be obtained by subsequently solving this linear differential equation and by applying a cascade of
-
Reduced minimax filtering by means of differential-algebraic equations
V. Mallet; S. Zhuk (Sergiy)
2011-01-01
htmlabstractA reduced minimax state estimation approach is proposed for high-dimensional models. It is based on the reduction of the ordinary differential equation with high state space dimension to the low-dimensional Differential-Algebraic Equation (DAE) and on the subsequent application of the
-
Modified Chebyshev Collocation Method for Solving Differential Equations
Directory of Open Access Journals (Sweden)
M Ziaul Arif
2015-05-01
Full Text Available This paper presents derivation of alternative numerical scheme for solving differential equations, which is modified Chebyshev (Vieta-Lucas Polynomial collocation differentiation matrices. The Scheme of modified Chebyshev (Vieta-Lucas Polynomial collocation method is applied to both Ordinary Differential Equations (ODEs and Partial Differential Equations (PDEs cases. Finally, the performance of the proposed method is compared with finite difference method and the exact solution of the example. It is shown that modified Chebyshev collocation method more effective and accurate than FDM for some example given.
-
International Conference on Differential and Difference Equations with Applications
Došlá, Zuzana; Došlý, Ondrej; Kloeden, Peter
2016-01-01
Aimed at the community of mathematicians working on ordinary and partial differential equations, difference equations, and functional equations, this book contains selected papers based on the presentations at the International Conference on Differential and Difference Equations and Applications (ICDDEA) 2015, dedicated to the memory of Professor Georg Sell. Contributions include new trends in the field of differential and difference equations, applications of differential and difference equations, as well as high-level survey results. The main aim of this recurring conference series is to promote, encourage, cooperate, and bring together researchers in the fields of differential and difference equations. All areas of differential and difference equations are represented, with special emphasis on applications.
-
Mathematical modeling based on ordinary differential equations: A promising approach to vaccinology.
Bonin, Carla Rezende Barbosa; Fernandes, Guilherme Cortes; Dos Santos, Rodrigo Weber; Lobosco, Marcelo
2017-02-01
New contributions that aim to accelerate the development or to improve the efficacy and safety of vaccines arise from many different areas of research and technology. One of these areas is computational science, which traditionally participates in the initial steps, such as the pre-screening of active substances that have the potential to become a vaccine antigen. In this work, we present another promising way to use computational science in vaccinology: mathematical and computational models of important cell and protein dynamics of the immune system. A system of Ordinary Differential Equations represents different immune system populations, such as B cells and T cells, antigen presenting cells and antibodies. In this way, it is possible to simulate, in silico, the immune response to vaccines under development or under study. Distinct scenarios can be simulated by varying parameters of the mathematical model. As a proof of concept, we developed a model of the immune response to vaccination against the yellow fever. Our simulations have shown consistent results when compared with experimental data available in the literature. The model is generic enough to represent the action of other diseases or vaccines in the human immune system, such as dengue and Zika virus.
-
Theoretical Analysis of Fas Ligand-Induced Apoptosis with an Ordinary Differential Equation Model.
Shi, Zhimin; Li, Yan; Liu, Zhihai; Mi, Jun; Wang, Renxiao
2012-12-01
Upon the treatment of Fas ligand, different types of cells exhibit different apoptotic mechanisms, which are determined by a complex network of biological pathways. In order to derive a quantitative interpretation of the cell sensitivity and apoptosis pathways, we have developed an ordinary differential equation model. Our model is intended to include all of the known major components in apoptosis pathways mediated by Fas receptor. It is composed of 29 equations using a total of 49 rate constants and 13 protein concentrations. All parameters used in our model were derived through nonlinear fitting to experimentally measured concentrations of four selected proteins in Jurkat T-cells, including caspase-3, caspase-8, caspase-9, and Bid. Our model is able to correctly interpret the role of kinetic parameters and protein concentrations in cell sensitivity to FasL. It reveals the possible reasons for the transition between type-I and type-II pathways and also provides some interesting predictions, such as the more decisive role of Fas over Bax in apoptosis pathway and a possible feedback mechanism between type-I and type-II pathways. But our model failed in predicting FasL-induced apoptotic mechanism of NCI-60 cells from their gene-expression levels. Limitations in our model are also discussed. Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
-
Algorithms For Integrating Nonlinear Differential Equations
Freed, A. D.; Walker, K. P.
1994-01-01
Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.
-
Trends in differential equations and applications
Neble, María; Galván, José
2016-01-01
This work collects the most important results presented at the Congress on Differential Equations and Applications/Congress on Applied Mathematics (CEDYA/CMA) in Cádiz (Spain) in 2015. It supports further research in differential equations, numerical analysis, mechanics, control and optimization. In particular, it helps readers gain an overview of specific problems of interest in the current mathematical research related to different branches of applied mathematics. This includes the analysis of nonlinear partial differential equations, exact solutions techniques for ordinary differential equations, numerical analysis and numerical simulation of some models arising in experimental sciences and engineering, control and optimization, and also trending topics on numerical linear Algebra, dynamical systems, and applied mathematics for Industry. This volume is mainly addressed to any researcher interested in the applications of mathematics, especially in any subject mentioned above. It may be also useful to PhD s...
-
A practical course in differential equations and mathematical modeling
Ibragimov , Nail H
2009-01-01
A Practical Course in Differential Equations and Mathematical Modelling is a unique blend of the traditional methods of ordinary and partial differential equations with Lie group analysis enriched by the author's own theoretical developments. The book which aims to present new mathematical curricula based on symmetry and invariance principles is tailored to develop analytic skills and working knowledge in both classical and Lie's methods for solving linear and nonlinear equations. This approach helps to make courses in differential equations, mathematical modelling, distributions and fundame
-
Introductory course on differential equations
Gorain, Ganesh C
2014-01-01
Introductory Course on DIFFERENTIAL EQUATIONS provides an excellent exposition of the fundamentals of ordinary and partial differential equations and is ideally suited for a first course of undergraduate students of mathematics, physics and engineering. The aim of this book is to present the elementary theories of differential equations in the forms suitable for use of those students whose main interest in the subject are based on simple mathematical ideas. KEY FEATURES: Discusses the subject in a systematic manner without sacrificing mathematical rigour. A variety of exercises drill the students in problem solving in view of the mathematical theories explained in the book. Worked out examples illustrated according to the theories developed in the book with possible alternatives. Exhaustive collection of problems and the simplicity of presentation differentiate this book from several others. Material contained will help teachers as well as aspiring students of different competitive examinations.
-
The numerical solution of linear multi-term fractional differential equations: systems of equations
Edwards, John T.; Ford, Neville J.; Simpson, A. Charles
2002-11-01
In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.
-
Numerical methods for differential equations and applications
International Nuclear Information System (INIS)
Ixaru, L.G.
1984-01-01
This book is addressed to persons who, without being professionals in applied mathematics, are often faced with the problem of numerically solving differential equations. In each of the first three chapters a definite class of methods is discussed for the solution of the initial value problem for ordinary differential equations: multistep methods; one-step methods; and piecewise perturbation methods. The fourth chapter is mainly focussed on the boundary value problems for linear second-order equations, with a section devoted to the Schroedinger equation. In the fifth chapter the eigenvalue problem for the radial Schroedinger equation is solved in several ways, with computer programs included. (Auth.)
-
A procedure to construct exact solutions of nonlinear fractional differential equations.
Güner, Özkan; Cevikel, Adem C
2014-01-01
We use the fractional transformation to convert the nonlinear partial fractional differential equations with the nonlinear ordinary differential equations. The Exp-function method is extended to solve fractional partial differential equations in the sense of the modified Riemann-Liouville derivative. We apply the Exp-function method to the time fractional Sharma-Tasso-Olver equation, the space fractional Burgers equation, and the time fractional fmKdV equation. As a result, we obtain some new exact solutions.
-
Engelhardt, Benjamin; Kschischo, Maik; Fröhlich, Holger
2017-06-01
Ordinary differential equations (ODEs) are a popular approach to quantitatively model molecular networks based on biological knowledge. However, such knowledge is typically restricted. Wrongly modelled biological mechanisms as well as relevant external influence factors that are not included into the model are likely to manifest in major discrepancies between model predictions and experimental data. Finding the exact reasons for such observed discrepancies can be quite challenging in practice. In order to address this issue, we suggest a Bayesian approach to estimate hidden influences in ODE-based models. The method can distinguish between exogenous and endogenous hidden influences. Thus, we can detect wrongly specified as well as missed molecular interactions in the model. We demonstrate the performance of our Bayesian dynamic elastic-net with several ordinary differential equation models from the literature, such as human JAK-STAT signalling, information processing at the erythropoietin receptor, isomerization of liquid α -Pinene, G protein cycling in yeast and UV-B triggered signalling in plants. Moreover, we investigate a set of commonly known network motifs and a gene-regulatory network. Altogether our method supports the modeller in an algorithmic manner to identify possible sources of errors in ODE-based models on the basis of experimental data. © 2017 The Author(s).
-
Integration of differential equations by the pseudo-linear (PL) approximation
International Nuclear Information System (INIS)
Bonalumi, Riccardo A.
1998-01-01
A new method of integrating differential equations was originated with the technique of approximately calculating the integrals called the pseudo-linear (PL) procedure: this method is A-stable. This article contains the following examples: 1st order ordinary differential equations (ODEs), 2nd order linear ODEs, stiff system of ODEs (neutron kinetics), one-dimensional parabolic (diffusion) partial differential equations. In this latter case, this PL method coincides with the Crank-Nicholson method
-
Sobolev gradients and differential equations
Neuberger, J W
2010-01-01
A Sobolev gradient of a real-valued functional on a Hilbert space is a gradient of that functional taken relative to an underlying Sobolev norm. This book shows how descent methods using such gradients allow a unified treatment of a wide variety of problems in differential equations. For discrete versions of partial differential equations, corresponding Sobolev gradients are seen to be vastly more efficient than ordinary gradients. In fact, descent methods with these gradients generally scale linearly with the number of grid points, in sharp contrast with the use of ordinary gradients. Aside from the first edition of this work, this is the only known account of Sobolev gradients in book form. Most of the applications in this book have emerged since the first edition was published some twelve years ago. What remains of the first edition has been extensively revised. There are a number of plots of results from calculations and a sample MatLab code is included for a simple problem. Those working through a fair p...
-
Wu, Hulin; Xue, Hongqi; Kumar, Arun
2012-06-01
Differential equations are extensively used for modeling dynamics of physical processes in many scientific fields such as engineering, physics, and biomedical sciences. Parameter estimation of differential equation models is a challenging problem because of high computational cost and high-dimensional parameter space. In this article, we propose a novel class of methods for estimating parameters in ordinary differential equation (ODE) models, which is motivated by HIV dynamics modeling. The new methods exploit the form of numerical discretization algorithms for an ODE solver to formulate estimating equations. First, a penalized-spline approach is employed to estimate the state variables and the estimated state variables are then plugged in a discretization formula of an ODE solver to obtain the ODE parameter estimates via a regression approach. We consider three different order of discretization methods, Euler's method, trapezoidal rule, and Runge-Kutta method. A higher-order numerical algorithm reduces numerical error in the approximation of the derivative, which produces a more accurate estimate, but its computational cost is higher. To balance the computational cost and estimation accuracy, we demonstrate, via simulation studies, that the trapezoidal discretization-based estimate is the best and is recommended for practical use. The asymptotic properties for the proposed numerical discretization-based estimators are established. Comparisons between the proposed methods and existing methods show a clear benefit of the proposed methods in regards to the trade-off between computational cost and estimation accuracy. We apply the proposed methods t an HIV study to further illustrate the usefulness of the proposed approaches. © 2012, The International Biometric Society.
-
On nonlinear differential equation with exact solutions having various pole orders
International Nuclear Information System (INIS)
Kudryashov, N.A.
2015-01-01
We consider a nonlinear ordinary differential equation having solutions with various movable pole order on the complex plane. We show that the pole order of exact solution is determined by values of parameters of the equation. Exact solutions in the form of the solitary waves for the second order nonlinear differential equation are found taking into account the method of the logistic function. Exact solutions of differential equations are discussed and analyzed
-
Fractional Complex Transform and exp-Function Methods for Fractional Differential Equations
Directory of Open Access Journals (Sweden)
Ahmet Bekir
2013-01-01
Full Text Available The exp-function method is presented for finding the exact solutions of nonlinear fractional equations. New solutions are constructed in fractional complex transform to convert fractional differential equations into ordinary differential equations. The fractional derivatives are described in Jumarie's modified Riemann-Liouville sense. We apply the exp-function method to both the nonlinear time and space fractional differential equations. As a result, some new exact solutions for them are successfully established.
-
Asymptotic behavior and stability of second order neutral delay differential equations
Chen, G.L.; van Gaans, O.W.; Verduyn Lunel, Sjoerd
2014-01-01
We study the asymptotic behavior of a class of second order neutral delay differential equations by both a spectral projection method and an ordinary differential equation method approach. We discuss the relation of these two methods and illustrate some features using examples. Furthermore, a fixed
-
Nonlinear differential equations with exact solutions expressed via the Weierstrass function
Kudryashov, NA
2004-01-01
A new problem is studied, that is to find nonlinear differential equations with special solutions expressed via the Weierstrass function. A method is discussed to construct nonlinear ordinary differential equations with exact solutions. The main step of our method is the assumption that nonlinear
-
International Nuclear Information System (INIS)
Boisseau, Bruno; Forgacs, Peter; Giacomini, Hector
2007-01-01
A new (algebraic) approximation scheme to find global solutions of two-point boundary value problems of ordinary differential equations (ODEs) is presented. The method is applicable for both linear and nonlinear (coupled) ODEs whose solutions are analytic near one of the boundary points. It is based on replacing the original ODEs by a sequence of auxiliary first-order polynomial ODEs with constant coefficients. The coefficients in the auxiliary ODEs are uniquely determined from the local behaviour of the solution in the neighbourhood of one of the boundary points. The problem of obtaining the parameters of the global (connecting) solutions, analytic at one of the boundary points, reduces to find the appropriate zeros of algebraic equations. The power of the method is illustrated by computing the approximate values of the 'connecting parameters' for a number of nonlinear ODEs arising in various problems in field theory. We treat in particular the static and rotationally symmetric global vortex, the skyrmion, the Abrikosov-Nielsen-Olesen vortex, as well as the 't Hooft-Polyakov magnetic monopole. The total energy of the skyrmion and of the monopole is also computed by the new method. We also consider some ODEs coming from the exact renormalization group. The ground-state energy level of the anharmonic oscillator is also computed for arbitrary coupling strengths with good precision. (fast track communication)
-
DeVille, R. E. Lee; Harkin, Anthony; Holzer, Matt; Josić, Krešimir; Kaper, Tasso J.
2008-06-01
For singular perturbation problems, the renormalization group (RG) method of Chen, Goldenfeld, and Oono [Phys. Rev. E. 49 (1994) 4502-4511] has been shown to be an effective general approach for deriving reduced or amplitude equations that govern the long time dynamics of the system. It has been applied to a variety of problems traditionally analyzed using disparate methods, including the method of multiple scales, boundary layer theory, the WKBJ method, the Poincaré-Lindstedt method, the method of averaging, and others. In this article, we show how the RG method may be used to generate normal forms for large classes of ordinary differential equations. First, we apply the RG method to systems with autonomous perturbations, and we show that the reduced or amplitude equations generated by the RG method are equivalent to the classical Poincaré-Birkhoff normal forms for these systems up to and including terms of O(ɛ2), where ɛ is the perturbation parameter. This analysis establishes our approach and generalizes to higher order. Second, we apply the RG method to systems with nonautonomous perturbations, and we show that the reduced or amplitude equations so generated constitute time-asymptotic normal forms, which are based on KBM averages. Moreover, for both classes of problems, we show that the main coordinate changes are equivalent, up to translations between the spaces in which they are defined. In this manner, our results show that the RG method offers a new approach for deriving normal forms for nonautonomous systems, and it offers advantages since one can typically more readily identify resonant terms from naive perturbation expansions than from the nonautonomous vector fields themselves. Finally, we establish how well the solution to the RG equations approximates the solution of the original equations on time scales of O(1/ɛ).
-
Lectures on the practical solution of differential equations
International Nuclear Information System (INIS)
Dresner, L.
1979-11-01
This report comprises lectures on the practical solution of ordinary and partial differential equations given in the In-Hours Continuing Education Program for Scientific and Technical Personnel at Oak Ridge National Laboratory
-
POSITIVE SOLUTIONS TO SEMI-LINEAR SECOND-ORDER ORDINARY DIFFERENTIAL EQUATIONS IN BANACH SPACE
Institute of Scientific and Technical Information of China (English)
无
2008-01-01
In this paper,we study the existence of positive periodic solution to some second- order semi-linear differential equation in Banach space.By the fixed point index theory, we prove that the semi-linear differential equation has two positive periodic solutions.
-
Prague, Mélanie; Commenges, Daniel; Guedj, Jérémie; Drylewicz, Julia; Thiébaut, Rodolphe
2013-08-01
Models based on ordinary differential equations (ODE) are widespread tools for describing dynamical systems. In biomedical sciences, data from each subject can be sparse making difficult to precisely estimate individual parameters by standard non-linear regression but information can often be gained from between-subjects variability. This makes natural the use of mixed-effects models to estimate population parameters. Although the maximum likelihood approach is a valuable option, identifiability issues favour Bayesian approaches which can incorporate prior knowledge in a flexible way. However, the combination of difficulties coming from the ODE system and from the presence of random effects raises a major numerical challenge. Computations can be simplified by making a normal approximation of the posterior to find the maximum of the posterior distribution (MAP). Here we present the NIMROD program (normal approximation inference in models with random effects based on ordinary differential equations) devoted to the MAP estimation in ODE models. We describe the specific implemented features such as convergence criteria and an approximation of the leave-one-out cross-validation to assess the model quality of fit. In pharmacokinetics models, first, we evaluate the properties of this algorithm and compare it with FOCE and MCMC algorithms in simulations. Then, we illustrate NIMROD use on Amprenavir pharmacokinetics data from the PUZZLE clinical trial in HIV infected patients. Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.
-
Institute of Scientific and Technical Information of China (English)
LI Shoufu
2005-01-01
A series of stability, contractivity and asymptotic stability results of the solutions to nonlinear stiff Volterra functional differential equations (VFDEs) in Banach spaces is obtained, which provides the unified theoretical foundation for the stability analysis of solutions to nonlinear stiff problems in ordinary differential equations(ODEs), delay differential equations(DDEs), integro-differential equations(IDEs) and VFDEs of other type which appear in practice.
-
International Nuclear Information System (INIS)
Burde, G.I.
2002-01-01
A new approach to the use of the Lie group technique for partial and ordinary differential equations dependent on a small parameter is developed. In addition to determining approximate solutions to the perturbed equation, the approach allows constructing integrable equations that have solutions with (partially) prescribed features. Examples of application of the approach to partial differential equations are given
-
Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE
Ansmann, Gerrit
2018-04-01
We present a family of Python modules for the numerical integration of ordinary, delay, or stochastic differential equations. The key features are that the user enters the derivative symbolically and it is just-in-time-compiled, allowing the user to efficiently integrate differential equations from a higher-level interpreted language. The presented modules are particularly suited for large systems of differential equations such as those used to describe dynamics on complex networks. Through the selected method of input, the presented modules also allow almost complete automatization of the process of estimating regular as well as transversal Lyapunov exponents for ordinary and delay differential equations. We conceptually discuss the modules' design, analyze their performance, and demonstrate their capabilities by application to timely problems.
-
Functional analysis in the study of differential and integral equations
International Nuclear Information System (INIS)
Sell, G.R.
1976-01-01
This paper illustrates the use of functional analysis in the study of differential equations. Our particular starting point, the theory of flows or dynamical systems, originated with the work of H. Poincare, who is the founder of the qualitative theory of ordinary differential equations. In the qualitative theory one tries to describe the behaviour of a solution, or a collection of solutions, without ''solving'' the differential equation. As a starting point one assumes the existence, and sometimes the uniqueness, of solutions and then one tries to describe the asymptotic behaviour, as time t→+infinity, of these solutions. We compare the notion of a flow with that of a C 0 -group of bounded linear operators on a Banach space. We shall show how the concept C 0 -group, or more generally a C 0 -semigroup, can be used to study the behaviour of solutions of certain differential and integral equations. Our main objective is to show how the concept of a C 0 -group and especially the notion of weak-compactness can be used to prove the existence of an invariant measure for a flow on a compact Hausdorff space. Applications to the theory of ordinary differential equations are included. (author)
-
Symmetries of stochastic differential equations: A geometric approach
Energy Technology Data Exchange (ETDEWEB)
De Vecchi, Francesco C., E-mail: francesco.devecchi@unimi.it; Ugolini, Stefania, E-mail: stefania.ugolini@unimi.it [Dipartimento di Matematica, Università degli Studi di Milano, via Saldini 50, Milano (Italy); Morando, Paola, E-mail: paola.morando@unimi.it [DISAA, Università degli Studi di Milano, via Celoria 2, Milano (Italy)
2016-06-15
A new notion of stochastic transformation is proposed and applied to the study of both weak and strong symmetries of stochastic differential equations (SDEs). The correspondence between an algebra of weak symmetries for a given SDE and an algebra of strong symmetries for a modified SDE is proved under suitable regularity assumptions. This general approach is applied to a stochastic version of a two dimensional symmetric ordinary differential equation and to the case of two dimensional Brownian motion.
-
Modeling and Prediction Using Stochastic Differential Equations
DEFF Research Database (Denmark)
Juhl, Rune; Møller, Jan Kloppenborg; Jørgensen, John Bagterp
2016-01-01
Pharmacokinetic/pharmakodynamic (PK/PD) modeling for a single subject is most often performed using nonlinear models based on deterministic ordinary differential equations (ODEs), and the variation between subjects in a population of subjects is described using a population (mixed effects) setup...... deterministic and can predict the future perfectly. A more realistic approach would be to allow for randomness in the model due to e.g., the model be too simple or errors in input. We describe a modeling and prediction setup which better reflects reality and suggests stochastic differential equations (SDEs...
-
Pythagoras, Binomial, and de Moivre Revisited Through Differential Equations
Singh, Jitender; Bajaj, Renu
2018-01-01
The classical Pythagoras theorem, binomial theorem, de Moivre's formula, and numerous other deductions are made using the uniqueness theorem for the initial value problems in linear ordinary differential equations.
-
Variations in the Solution of Linear First-Order Differential Equations. Classroom Notes
Seaman, Brian; Osler, Thomas J.
2004-01-01
A special project which can be given to students of ordinary differential equations is described in detail. Students create new differential equations by changing the dependent variable in the familiar linear first-order equation (dv/dx)+p(x)v=q(x) by means of a substitution v=f(y). The student then creates a table of the new equations and…
-
ODEPACK, Initial Value Problems of Ordinary Differential Equation System
International Nuclear Information System (INIS)
Hindmarsh, A.C.; Petzold, L.R.
2005-01-01
I - Description of program or function: ODEPACK is a collection of Fortran solvers for the initial value problem for ordinary differential equation systems. It consists of nine solvers, namely a basic solver called LSODE and eight variants of it -- LSODES, LSODA, LSODAR, LSODPK, LSODKR, LSODI, LSOIBT, and LSODIS. The collection is suitable for both stiff and non-stiff systems. It includes solvers for systems given in explicit form, dy/dt = f(t,y), and also solvers for systems given in linearly implicit form, A(t,y) dy/dt = g(t,y). Two of the solvers use general sparse matrix solvers for the linear systems that arise. Two others use iterative (preconditioned Krylov) methods instead of direct methods for these linear systems. The most recent addition is LSODIS, which solves implicit problems with general sparse treatment of all matrices involved. The ODEPACK solvers are written in standard Fortran 77, with a few exceptions, and with minimal machine dependencies. There are separate double and single precision versions of ODEPACK. The actual solver names are those given above with a prefix of D- or S- for the double or single precision version, respectively, i.e. DLSODE/SLSODE, etc. Each solver consists of a main driver subroutine having the same name as the solver and some number of subordinate routines. For each solver, there is also a demonstration program, which solves one or two simple problems in a somewhat self-checking manner. A. Solvers for explicitly given systems. For each of the following solvers, it is assumed that the ODEs are given explicitly, so that the system can be written in the form dy/dt = f(t,y), where y is the vector of dependent variables, and t is the independent variable. 1. LSODE (Livermore Solver for Ordinary Differential Equations) is the basic solver of the collection. It solves stiff and non-stiff systems of the form dy/dt = f. In the stiff case, it treats the Jacobian matrix df/dy as either a dense (full) or a banded matrix, and as
-
Regarding on the exact solutions for the nonlinear fractional differential equations
Directory of Open Access Journals (Sweden)
Kaplan Melike
2016-01-01
Full Text Available In this work, we have considered the modified simple equation (MSE method for obtaining exact solutions of nonlinear fractional-order differential equations. The space-time fractional equal width (EW and the modified equal width (mEW equation are considered for illustrating the effectiveness of the algorithm. It has been observed that all exact solutions obtained in this paper verify the nonlinear ordinary differential equations which was obtained from nonlinear fractional-order differential equations under the terms of wave transformation relationship. The obtained results are shown graphically.
-
Subroutine for series solutions of linear differential equations
International Nuclear Information System (INIS)
Tasso, H.; Steuerwald, J.
1976-02-01
A subroutine for Taylor series solutions of systems of ordinary linear differential equations is descriebed. It uses the old idea of Lie series but allows simple implementation and is time-saving for symbolic manipulations. (orig.) [de
-
The Local Brewery: A Project for Use in Differential Equations Courses
Starling, James K.; Povich, Timothy J.; Findlay, Michael
2016-01-01
We describe a modeling project designed for an ordinary differential equations (ODEs) course using first-order and systems of first-order differential equations to model the fermentation process in beer. The project aims to expose the students to the modeling process by creating and solving a mathematical model and effectively communicating their…
-
Lie symmetries and differential galois groups of linear equations
Oudshoorn, W.R.; Put, M. van der
2002-01-01
For a linear ordinary differential equation the Lie algebra of its infinitesimal Lie symmetries is compared with its differential Galois group. For this purpose an algebraic formulation of Lie symmetries is developed. It turns out that there is no direct relation between the two above objects. In
-
Partial differential equations & boundary value problems with Maple
Articolo, George A
2009-01-01
Partial Differential Equations and Boundary Value Problems with Maple presents all of the material normally covered in a standard course on partial differential equations, while focusing on the natural union between this material and the powerful computational software, Maple. The Maple commands are so intuitive and easy to learn, students can learn what they need to know about the software in a matter of hours- an investment that provides substantial returns. Maple''s animation capabilities allow students and practitioners to see real-time displays of the solutions of partial differential equations. Maple files can be found on the books website. Ancillary list: Maple files- http://www.elsevierdirect.com/companion.jsp?ISBN=9780123747327 Provides a quick overview of the software w/simple commands needed to get startedIncludes review material on linear algebra and Ordinary Differential equations, and their contribution in solving partial differential equationsIncorporates an early introduction to Sturm-L...
-
Samples of noncommutative products in certain differential equations
International Nuclear Information System (INIS)
Legare, M
2010-01-01
A set of associative noncommutative products is considered in different differential equations of the ordinary and partial types. A method of separation of variables is considered for a large set of those systems. The products involved include for example some * products and some products based on Nijenhuis tensors, which are embedded in the differential equations of the Laplace/Poisson, Lax and Schroedinger styles. A comment on the *-products of Reshetikhin-Jambor-Sykora type is also given in relation to *-products of Vey type.
-
Surveys in differential-algebraic equations IV
Reis, Timo
2017-01-01
The present volume comprises survey articles on various fields of Differential-Algebraic Equations (DAEs) which have widespread applications in controlled dynamical systems, especially in mechanical and electrical engineering and a strong relation to (ordinary) differential equations. The individual chapters provide reviews, presentations of the current state of research and new concepts in - History of DAEs - DAE aspects of mechanical multibody systems - Model reduction of DAEs - Observability for DAEs - Numerical Analysis for DAEs The results are presented in an accessible style, making this book suitable not only for active researchers but also for graduate students (with a good knowledge of the basic principles of DAEs) for self-study.
-
Surveys in differential-algebraic equations III
Reis, Timo
2015-01-01
The present volume comprises survey articles on various fields of Differential-Algebraic Equations (DAEs), which have widespread applications in controlled dynamical systems, especially in mechanical and electrical engineering and a strong relation to (ordinary) differential equations. The individual chapters provide reviews, presentations of the current state of research and new concepts in - Flexibility of DAE formulations - Reachability analysis and deterministic global optimization - Numerical linear algebra methods - Boundary value problems The results are presented in an accessible style, making this book suitable not only for active researchers but also for graduate students (with a good knowledge of the basic principles of DAEs) for self-study.
-
Exact solutions of nonlinear fractional differential equations by (G′/G)-expansion method
International Nuclear Information System (INIS)
Bekir Ahmet; Güner Özkan
2013-01-01
In this paper, we use the fractional complex transform and the (G′/G)-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie's modified Riemann—Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations
-
Zhang, Xinyu; Cao, Jiguo; Carroll, Raymond J
2015-03-01
We consider model selection and estimation in a context where there are competing ordinary differential equation (ODE) models, and all the models are special cases of a "full" model. We propose a computationally inexpensive approach that employs statistical estimation of the full model, followed by a combination of a least squares approximation (LSA) and the adaptive Lasso. We show the resulting method, here called the LSA method, to be an (asymptotically) oracle model selection method. The finite sample performance of the proposed LSA method is investigated with Monte Carlo simulations, in which we examine the percentage of selecting true ODE models, the efficiency of the parameter estimation compared to simply using the full and true models, and coverage probabilities of the estimated confidence intervals for ODE parameters, all of which have satisfactory performances. Our method is also demonstrated by selecting the best predator-prey ODE to model a lynx and hare population dynamical system among some well-known and biologically interpretable ODE models. © 2014, The International Biometric Society.
-
DEFF Research Database (Denmark)
Tornøe, Christoffer Wenzel; Overgaard, Rune Viig; Agerso, H.
2005-01-01
of noise: a measurement and a system noise term. The measurement noise represents uncorrelated error due to, for example, assay error while the system noise accounts for structural misspecifications, approximations of the dynamical model, and true random physiological fluctuations. Since the system noise...... degarelix. Conclusions. The EKF-based algorithm was successfully implemented in NONMEM for parameter estimation in population PK/PD models described by systems of SDEs. The example indicated that it was possible to pinpoint structural model deficiencies, and that valuable information may be obtained......Purpose. The objective of the present analysis was to explore the use of stochastic differential equations (SDEs) in population pharmacokinetic/pharmacodynamic (PK/PD) modeling. Methods. The intra-individual variability in nonlinear mixed-effects models based on SDEs is decomposed into two types...
-
A New Factorisation of a General Second Order Differential Equation
Clegg, Janet
2006-01-01
A factorisation of a general second order ordinary differential equation is introduced from which the full solution to the equation can be obtained by performing two integrations. The method is compared with traditional methods for solving these type of equations. It is shown how the Green's function can be derived directly from the factorisation…
-
First order linear ordinary differential equations in associative algebras
Directory of Open Access Journals (Sweden)
Gordon Erlebacher
2004-01-01
Full Text Available In this paper, we study the linear differential equation $$ frac{dx}{dt}=sum_{i=1}^n a_i(t x b_i(t + f(t $$ in an associative but non-commutative algebra $mathcal{A}$, where the $b_i(t$ form a set of commuting $mathcal{A}$-valued functions expressed in a time-independent spectral basis consisting of mutually annihilating idempotents and nilpotents. Explicit new closed solutions are derived, and examples are presented to illustrate the theory.
-
An accurate scheme by block method for third order ordinary ...
African Journals Online (AJOL)
problems of ordinary differential equations is presented in this paper. The approach of collocation approximation is adopted in the derivation of the scheme and then the scheme is applied as simultaneous integrator to special third order initial value problem of ordinary differential equations. This implementation strategy is ...
-
Universal formats for nonlinear ordinary differential systems
International Nuclear Information System (INIS)
Kerner, E.H.
1981-01-01
It is shown that very general nonlinear ordinary differential systems (embracing all that arise in practice) may, first, be brought down to polynomial systems (where the nonlinearities occur only as polynomials in the dependent variables) by introducing suitable new variables into the original system; second, that polynomial systems are reducible to ''Riccati systems,'' where the nonlinearities are quadratic at most; third, that Riccati systems may be brought to elemental universal formats containing purely quadratic terms with simple arrays of coefficients that are all zero or unity. The elemental systems have representations as novel types of matrix Riccati equations. Different starting systems and their associated Riccati systems differ from one another, at the final elemental level, in order and in initial data, but not in format
-
International Nuclear Information System (INIS)
Zhang Huiqun
2009-01-01
By using some exact solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct the exact complex solutions for nonlinear partial differential equations. The method is implemented for the NLS equation, a new Hamiltonian amplitude equation, the coupled Schrodinger-KdV equations and the Hirota-Maccari equations. New exact complex solutions are obtained.
-
Modelling Evolutionary Algorithms with Stochastic Differential Equations.
Heredia, Jorge Pérez
2017-11-20
There has been renewed interest in modelling the behaviour of evolutionary algorithms (EAs) by more traditional mathematical objects, such as ordinary differential equations or Markov chains. The advantage is that the analysis becomes greatly facilitated due to the existence of well established methods. However, this typically comes at the cost of disregarding information about the process. Here, we introduce the use of stochastic differential equations (SDEs) for the study of EAs. SDEs can produce simple analytical results for the dynamics of stochastic processes, unlike Markov chains which can produce rigorous but unwieldy expressions about the dynamics. On the other hand, unlike ordinary differential equations (ODEs), they do not discard information about the stochasticity of the process. We show that these are especially suitable for the analysis of fixed budget scenarios and present analogues of the additive and multiplicative drift theorems from runtime analysis. In addition, we derive a new more general multiplicative drift theorem that also covers non-elitist EAs. This theorem simultaneously allows for positive and negative results, providing information on the algorithm's progress even when the problem cannot be optimised efficiently. Finally, we provide results for some well-known heuristics namely Random Walk (RW), Random Local Search (RLS), the (1+1) EA, the Metropolis Algorithm (MA), and the Strong Selection Weak Mutation (SSWM) algorithm.
-
The ATOMFT integrator - Using Taylor series to solve ordinary differential equations
Berryman, Kenneth W.; Stanford, Richard H.; Breckheimer, Peter J.
1988-01-01
This paper discusses the application of ATOMFT, an integration package based on Taylor series solution with a sophisticated user interface. ATOMFT has the capabilities to allow the implementation of user defined functions and the solution of stiff and algebraic equations. Detailed examples, including the solutions to several astrodynamics problems, are presented. Comparisons with its predecessor ATOMCC and other modern integrators indicate that ATOMFT is a fast, accurate, and easy method to use to solve many differential equation problems.
-
ON ASYMTOTIC APPROXIMATIONS OF FIRST INTEGRALS FOR DIFFERENTIAL AND DIFFERENCE EQUATIONS
Directory of Open Access Journals (Sweden)
W.T. van Horssen
2007-04-01
Full Text Available In this paper the concept of integrating factors for differential equations and the concept of invariance factors for difference equations to obtain first integrals or invariants will be presented. It will be shown that all integrating factors have to satisfya system of partial differential equations, and that all invariance factors have to satisfy a functional equation. In the period 1997-2001 a perturbation method based on integrating vectors was developed to approximate first integrals for systems of ordinary differential equations. This perturbation method will be reviewed shortly. Also in the paper the first results in the development of a perturbation method for difference equations based on invariance factors will be presented.
-
Dynamic data analysis modeling data with differential equations
Ramsay, James
2017-01-01
This text focuses on the use of smoothing methods for developing and estimating differential equations following recent developments in functional data analysis and building on techniques described in Ramsay and Silverman (2005) Functional Data Analysis. The central concept of a dynamical system as a buffer that translates sudden changes in input into smooth controlled output responses has led to applications of previously analyzed data, opening up entirely new opportunities for dynamical systems. The technical level has been kept low so that those with little or no exposure to differential equations as modeling objects can be brought into this data analysis landscape. There are already many texts on the mathematical properties of ordinary differential equations, or dynamic models, and there is a large literature distributed over many fields on models for real world processes consisting of differential equations. However, a researcher interested in fitting such a model to data, or a statistician interested in...
-
Differential Equations Models to Study Quorum Sensing.
Pérez-Velázquez, Judith; Hense, Burkhard A
2018-01-01
Mathematical models to study quorum sensing (QS) have become an important tool to explore all aspects of this type of bacterial communication. A wide spectrum of mathematical tools and methods such as dynamical systems, stochastics, and spatial models can be employed. In this chapter, we focus on giving an overview of models consisting of differential equations (DE), which can be used to describe changing quantities, for example, the dynamics of one or more signaling molecule in time and space, often in conjunction with bacterial growth dynamics. The chapter is divided into two sections: ordinary differential equations (ODE) and partial differential equations (PDE) models of QS. Rates of change are represented mathematically by derivatives, i.e., in terms of DE. ODE models allow describing changes in one independent variable, for example, time. PDE models can be used to follow changes in more than one independent variable, for example, time and space. Both types of models often consist of systems (i.e., more than one equation) of equations, such as equations for bacterial growth and autoinducer concentration dynamics. Almost from the onset, mathematical modeling of QS using differential equations has been an interdisciplinary endeavor and many of the works we revised here will be placed into their biological context.
-
Application of the Generalized Differential Quadrature Method in Solving Burgers' Equations
International Nuclear Information System (INIS)
Mokhtari, R.; Toodar, A. Samadi; Chegini, N.G.
2011-01-01
The aim of this paper is to obtain numerical solutions of the one-dimensional, two-dimensional and coupled Burgers' equations through the generalized differential quadrature method (GDQM). The polynomial-based differential quadrature (PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta (TVD-RK) method. The numerical solutions are satisfactorily coincident with the exact solutions. The method can compete against the methods applied in the literature. (general)
-
Methods of mathematical modelling continuous systems and differential equations
Witelski, Thomas
2015-01-01
This book presents mathematical modelling and the integrated process of formulating sets of equations to describe real-world problems. It describes methods for obtaining solutions of challenging differential equations stemming from problems in areas such as chemical reactions, population dynamics, mechanical systems, and fluid mechanics. Chapters 1 to 4 cover essential topics in ordinary differential equations, transport equations and the calculus of variations that are important for formulating models. Chapters 5 to 11 then develop more advanced techniques including similarity solutions, matched asymptotic expansions, multiple scale analysis, long-wave models, and fast/slow dynamical systems. Methods of Mathematical Modelling will be useful for advanced undergraduate or beginning graduate students in applied mathematics, engineering and other applied sciences.
-
Avellar, J.; Duarte, L. G. S.; da Mota, L. A. C. P.
2012-10-01
We present a set of software routines in Maple 14 for solving first order ordinary differential equations (FOODEs). The package implements the Prelle-Singer method in its original form together with its extension to include integrating factors in terms of elementary functions. The package also presents a theoretical extension to deal with all FOODEs presenting Liouvillian solutions. Applications to ODEs taken from standard references show that it solves ODEs which remain unsolved using Maple's standard ODE solution routines. New version program summary Program title: PSsolver Catalogue identifier: ADPR_v2_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADPR_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 2302 No. of bytes in distributed program, including test data, etc.: 31962 Distribution format: tar.gz Programming language: Maple 14 (also tested using Maple 15 and 16). Computer: Intel Pentium Processor P6000, 1.86 GHz. Operating system: Windows 7. RAM: 4 GB DDR3 Memory Classification: 4.3. Catalogue identifier of previous version: ADPR_v1_0 Journal reference of previous version: Comput. Phys. Comm. 144 (2002) 46 Does the new version supersede the previous version?: Yes Nature of problem: Symbolic solution of first order differential equations via the Prelle-Singer method. Solution method: The method of solution is based on the standard Prelle-Singer method, with extensions for the cases when the FOODE contains elementary functions. Additionally, an extension of our own which solves FOODEs with Liouvillian solutions is included. Reasons for new version: The program was not running anymore due to changes in the latest versions of Maple. Additionally, we corrected/changed some bugs/details that were hampering the smoother functioning of the routines. Summary
-
DEFF Research Database (Denmark)
Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode
2009-01-01
are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE1 approximation......The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model...... development, J. Pharmacokinet. Pharmacodyn. 32 (February(l)) (2005) 109-141; C.W. Tornoe, R.V Overgaard, H. Agerso, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8...
-
Analytical approach for the Floquet theory of delay differential equations.
Simmendinger, C; Wunderlin, A; Pelster, A
1999-05-01
We present an analytical approach to deal with nonlinear delay differential equations close to instabilities of time periodic reference states. To this end we start with approximately determining such reference states by extending the Poincaré-Lindstedt and the Shohat expansions, which were originally developed for ordinary differential equations. Then we systematically elaborate a linear stability analysis around a time periodic reference state. This allows us to approximately calculate the Floquet eigenvalues and their corresponding eigensolutions by using matrix valued continued fractions.
-
A Second-Year Undergraduate Course in Applied Differential Equations.
Fahidy, Thomas Z.
1991-01-01
Presents the framework for a chemical engineering course using ordinary differential equations to solve problems with the underlying strategy of concisely discussing the theory behind each solution technique without extensions to formal proofs. Includes typical class illustrations, student responses to this strategy, and reaction of the…
-
On Robust Stability of Systems of Differential-Algebraic Equations
Directory of Open Access Journals (Sweden)
A. Shcheglova
2016-06-01
The sufficient conditions of robust stability for index-one and index-two systems are obtained. We use the values of real and complex stability radii obtained for system of ordinary differential equations solved with respect to the derivatives. We consider the example illustrating the obtained results.
-
An ordinary differential equation model for full thickness wounds and the effects of diabetes.
Bowden, L G; Maini, P K; Moulton, D E; Tang, J B; Wang, X T; Liu, P Y; Byrne, H M
2014-11-21
Wound healing is a complex process in which a sequence of interrelated phases contributes to a reduction in wound size. For diabetic patients, many of these processes are compromised, so that wound healing slows down. In this paper we present a simple ordinary differential equation model for wound healing in which attention focusses on the dominant processes that contribute to closure of a full thickness wound. Asymptotic analysis of the resulting model reveals that normal healing occurs in stages: the initial and rapid elastic recoil of the wound is followed by a longer proliferative phase during which growth in the dermis dominates healing. At longer times, fibroblasts exert contractile forces on the dermal tissue, the resulting tension stimulating further dermal tissue growth and enhancing wound closure. By fitting the model to experimental data we find that the major difference between normal and diabetic healing is a marked reduction in the rate of dermal tissue growth for diabetic patients. The model is used to estimate the breakdown of dermal healing into two processes: tissue growth and contraction, the proportions of which provide information about the quality of the healed wound. We show further that increasing dermal tissue growth in the diabetic wound produces closure times similar to those associated with normal healing and we discuss the clinical implications of this hypothesised treatment. Copyright © 2014 Elsevier Ltd. All rights reserved.
-
Differential equations and integrable models: the SU(3) case
International Nuclear Information System (INIS)
Dorey, Patrick; Tateo, Roberto
2000-01-01
We exhibit a relationship between the massless a 2 (2) integrable quantum field theory and a certain third-order ordinary differential equation, thereby extending a recent result connecting the massless sine-Gordon model to the Schroedinger equation. This forms part of a more general correspondence involving A 2 -related Bethe ansatz systems and third-order differential equations. A non-linear integral equation for the generalised spectral problem is derived, and some numerical checks are performed. Duality properties are discussed, and a simple variant of the non-linear equation is suggested as a candidate to describe the finite volume ground state energies of minimal conformal field theories perturbed by the operators phi 12 , phi 21 and phi 15 . This is checked against previous results obtained using the thermodynamic Bethe ansatz
-
International Nuclear Information System (INIS)
Yao Ruo-Xia; Wang Wei; Chen Ting-Hua
2014-01-01
Motivated by the widely used ansätz method and starting from the modified Riemann—Liouville derivative together with a fractional complex transformation that can be utilized to transform nonlinear fractional partial differential equations to nonlinear ordinary differential equations, new types of exact traveling wave solutions to three important nonlinear space- and time-fractional partial differential equations are obtained simultaneously in terms of solutions of a Riccati equation. The results are new and first reported in this paper. (general)
-
Solving differential equations with unknown constitutive relations as recurrent neural networks
Energy Technology Data Exchange (ETDEWEB)
Hagge, Tobias J.; Stinis, Panagiotis; Yeung, Enoch H.; Tartakovsky, Alexandre M.
2017-12-08
We solve a system of ordinary differential equations with an unknown functional form of a sink (reaction rate) term. We assume that the measurements (time series) of state variables are partially available, and use a recurrent neural network to “learn” the reaction rate from this data. This is achieved by including discretized ordinary differential equations as part of a recurrent neural network training problem. We extend TensorFlow’s recurrent neural network architecture to create a simple but scalable and effective solver for the unknown functions, and apply it to a fedbatch bioreactor simulation problem. Use of techniques from recent deep learning literature enables training of functions with behavior manifesting over thousands of time steps. Our networks are structurally similar to recurrent neural networks, but differ in purpose, and require modified training strategies.
-
Differential equation for genus-two characters in arbitrary rational conformal field theories
International Nuclear Information System (INIS)
Mathur, S.D.; Sen, A.
1989-01-01
We develop a general method for deriving ordinary differential equations for the genus-two ''characters'' of an arbitrary rational conformal field theory using the hyperelliptic representation of the genus-two moduli space. We illustrate our method by explicitly deriving the character differential equations for k=1 SU(2), G 2 , and F 4 WZW models. Our method provides an intrinsic definition of conformal field theories on higher genus Riemann surfaces. (orig.)
-
Domoshnitsky, Alexander; Maghakyan, Abraham; Berezansky, Leonid
2017-01-01
In this paper a method for studying stability of the equation [Formula: see text] not including explicitly the first derivative is proposed. We demonstrate that although the corresponding ordinary differential equation [Formula: see text] is not exponentially stable, the delay equation can be exponentially stable.
-
Solving differential-algebraic equation systems by means of index reduction methodology
DEFF Research Database (Denmark)
Sørensen, Kim; Houbak, Niels; Condra, Thomas Joseph
2006-01-01
of a number of differential equations and algebraic equations - a so called DAE system. Two of the DAE systems are of index 1 and they can be solved by means of standard DAE-solvers. For the actual application, the equation systems are integrated by means of MATLAB’s solver: ode23t, that solves moderately...... stiff ODE’s and index 1 DAE’s by means of the trapezoidal rule. The last sub-model that models the boilers steam drum consist of two differential and three algebraic equations. The index of this model is greater than 1, which means that ode23t cannot integrate this equation system. In this paper......, it is shown how the equation system, by means of an index reduction methodology, can be reduced to a system of Ordinary- Differential-Equations - ODE’s....
-
Pratt, D. T.
1984-01-01
Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.
-
Numerical solution of second-order stochastic differential equations with Gaussian random parameters
Directory of Open Access Journals (Sweden)
Rahman Farnoosh
2014-07-01
Full Text Available In this paper, we present the numerical solution of ordinary differential equations (or SDEs, from each orderespecially second-order with time-varying and Gaussian random coefficients. We indicate a complete analysisfor second-order equations in specially case of scalar linear second-order equations (damped harmonicoscillators with additive or multiplicative noises. Making stochastic differential equations system from thisequation, it could be approximated or solved numerically by different numerical methods. In the case oflinear stochastic differential equations system by Computing fundamental matrix of this system, it could becalculated based on the exact solution of this system. Finally, this stochastic equation is solved by numericallymethod like E.M. and Milstein. Also its Asymptotic stability and statistical concepts like expectationand variance of solutions are discussed.
-
Directory of Open Access Journals (Sweden)
Meng-Rong Li
2018-01-01
Full Text Available Considering the phenomenon of the mean reversion and the different speeds of stock prices in the bull market and in the bear market, we propose four dynamic models each of which is represented by a parameterized ordinary differential equation in this study. Based on existing studies, the models are in the form of either the logistic growth or the law of Newton’s cooling. We solve the models by dynamic integration and apply them to the daily closing prices of the Taiwan stock index, Taiwan Stock Exchange Capitalization Weighted Stock Index. The empirical study shows that some of the models fit the prices well and the forecasting ability of the best model is acceptable even though the martingale forecasts the prices slightly better. To increase the forecasting ability and to broaden the scope of applications of the dynamic models, we will model the coefficients of the dynamic models in the future. Applying the models to the market without the price limit is also our future work.
-
Collage-based approaches for elliptic partial differential equations inverse problems
Yodzis, Michael; Kunze, Herb
2017-01-01
The collage method for inverse problems has become well-established in the literature in recent years. Initial work developed a collage theorem, based upon Banach's fixed point theorem, for treating inverse problems for ordinary differential equations (ODEs). Amongst the subsequent work was a generalized collage theorem, based upon the Lax-Milgram representation theorem, useful for treating inverse problems for elliptic partial differential equations (PDEs). Each of these two different approaches can be applied to elliptic PDEs in one space dimension. In this paper, we explore and compare how the two different approaches perform for the estimation of the diffusivity for a steady-state heat equation.
-
Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode; Overgaard, Rune Viig; Madsen, Henrik
2009-06-01
The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model development, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 109-141; C.W. Tornøe, R.V. Overgaard, H. Agersø, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8)) (2005) 1247-1258; R.V. Overgaard, N. Jonsson, C.W. Tornøe, H. Madsen, Non-linear mixed-effects models with stochastic differential equations: implementation of an estimation algorithm, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 85-107; U. Picchini, S. Ditlevsen, A. De Gaetano, Maximum likelihood estimation of a time-inhomogeneous stochastic differential model of glucose dynamics, Math. Med. Biol. 25 (June(2)) (2008) 141-155]. PK/PD models are traditionally based ordinary differential equations (ODEs) with an observation link that incorporates noise. This state-space formulation only allows for observation noise and not for system noise. Extending to SDEs allows for a Wiener noise component in the system equations. This additional noise component enables handling of autocorrelated residuals originating from natural variation or systematic model error. Autocorrelated residuals are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE(1) approximation to the population likelihood which is generated from the individual likelihoods that are approximated using the Extended Kalman Filter's one-step predictions.
-
Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review
Kennedy, Christopher A.; Carpenter, Mark H.
2016-01-01
A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.
-
Zero singularities of codimension two and three in delay differential equations
International Nuclear Information System (INIS)
Campbell, Sue Ann; Yuan Yuan
2008-01-01
We give conditions under which a general class of delay differential equations has a point of Bogdanov–Takens or a triple zero bifurcation. We show how a centre manifold projection of the delay equations reduces the dynamics to two- or three-dimensional systems of ordinary differential equations. We put these equations in normal form and determine how the coefficients of the normal forms depend on the original parameters in the model. Finally we apply our results to two neural models and compare the predictions of the theory with numerical bifurcation analysis of the full equations. One model involves a transcritical bifurcation, hence we derive and analyse the appropriate unfoldings for this case
-
Neural network error correction for solving coupled ordinary differential equations
Shelton, R. O.; Darsey, J. A.; Sumpter, B. G.; Noid, D. W.
1992-01-01
A neural network is presented to learn errors generated by a numerical algorithm for solving coupled nonlinear differential equations. The method is based on using a neural network to correctly learn the error generated by, for example, Runge-Kutta on a model molecular dynamics (MD) problem. The neural network programs used in this study were developed by NASA. Comparisons are made for training the neural network using backpropagation and a new method which was found to converge with fewer iterations. The neural net programs, the MD model and the calculations are discussed.
-
International Nuclear Information System (INIS)
Elmhirst, Toby; Stewart, Ian; Doebeli, Michael
2008-01-01
We present a class of systems of ordinary differential equations (ODEs), which we call 'pod systems', that offers a new perspective on dynamical systems defined on a spatial domain. Such systems are typically studied as partial differential equations, but pod systems bring the analytic techniques of ODE theory to bear on the problems, and are thus able to study behaviours and bifurcations that are not easily accessible to the standard methods. In particular, pod systems are specifically designed to study spatial dynamical systems that exhibit multi-modal solutions. A pod system is essentially a linear combination of parametrized functions in which the coefficients and parameters are variables whose dynamics are specified by a system of ODEs. That is, pod systems are concerned with the dynamics of functions of the form Ψ(s, t) = y 1 (t) φ(s; x 1 (t)) + ··· + y N (t) φ(s; x N (t)), where s in R n is the spatial variable and φ: R n × R d → R. The parameters x i in R d and coefficients y i in R are dynamic variables which evolve according to some system of ODEs, x-dot i = G i (x, y) and y-dot i = H i (x, y), for i = 1, ..., N. The dynamics of Ψ in function space can then be studied through the dynamics of the x and y in finite dimensions. A vital feature of pod systems is that the ODEs that specify the dynamics of the x and y variables are not arbitrary; restrictions on G i and H i are required to guarantee that the dynamics of Ψ in function space are well defined (that is, that trajectories are unique). One important restriction is symmetry in the ODEs which arises because Ψ is invariant under permutations of the indices of the (x i , y i ) pairs. However, this is not the whole story, and the primary goal of this paper is to determine the necessary structure of the ODEs explicitly to guarantee that the dynamics of Ψ are well defined
-
Directory of Open Access Journals (Sweden)
Eman Ali Hussain
2015-01-01
Full Text Available Absract In this project A new method for solving Stochastic Differential Equations SDEs deriving by Wiener process numerically will be construct and implement using Accelerated Genetic Algorithm AGA. An SDE is a differential equation in which one or more of the terms and hence the solutions itself is a stochastic process. Solving stochastic differential equations requires going away from the recognizable deterministic setting of ordinary and partial differential equations into a world where the evolution of a quantity has an inherent random component and where the expected behavior of this quantity can be described in terms of probability distributions. We applied our method on the Ito formula which is equivalent to the SDE to find approximation solution of the SDEs. Numerical experiments illustrate the behavior of the proposed method.
-
Dorodnitsyn, Vladimir A.; Kozlov, Roman; Meleshko, Sergey V.; Winternitz, Pavel
2018-05-01
A recent article was devoted to an analysis of the symmetry properties of a class of first-order delay ordinary differential systems (DODSs). Here we concentrate on linear DODSs, which have infinite-dimensional Lie point symmetry groups due to the linear superposition principle. Their symmetry algebra always contains a two-dimensional subalgebra realized by linearly connected vector fields. We identify all classes of linear first-order DODSs that have additional symmetries, not due to linearity alone, and we present representatives of each class. These additional symmetries are then used to construct exact analytical particular solutions using symmetry reduction.
-
Improved stochastic approximation methods for discretized parabolic partial differential equations
Guiaş, Flavius
2016-12-01
We present improvements of the stochastic direct simulation method, a known numerical scheme based on Markov jump processes which is used for approximating solutions of ordinary differential equations. This scheme is suited especially for spatial discretizations of evolution partial differential equations (PDEs). By exploiting the full path simulation of the stochastic method, we use this first approximation as a predictor and construct improved approximations by Picard iterations, Runge-Kutta steps, or a combination. This has as consequence an increased order of convergence. We illustrate the features of the improved method at a standard benchmark problem, a reaction-diffusion equation modeling a combustion process in one space dimension (1D) and two space dimensions (2D).
-
Furihata, Daisuke
2010-01-01
Nonlinear Partial Differential Equations (PDEs) have become increasingly important in the description of physical phenomena. Unlike Ordinary Differential Equations, PDEs can be used to effectively model multidimensional systems. The methods put forward in Discrete Variational Derivative Method concentrate on a new class of ""structure-preserving numerical equations"" which improves the qualitative behaviour of the PDE solutions and allows for stable computing. The authors have also taken care to present their methods in an accessible manner, which means that the book will be useful to engineer
-
(Ln-bar, g)-spaces. Ordinary and tensor differentials
International Nuclear Information System (INIS)
Manoff, S.; Dimitrov, B.
1998-01-01
Different types of differentials as special cases of differential operators acting on tensor fields over (L n bar, g)-spaces are considered. The ordinary differential, the covariant differential as a special case of the covariant differential operator, and the Lie differential as a special case of the Lie differential operator are investigated. The tensor differential and its special types (Covariant tensor differential, and Lie tensor differential) are determined and their properties are discussed. Covariant symmetric and antisymmetric (external) tensor differentials, Lie symmetric, and Lie antisymmetric (external) tensor differentials are determined and considered over (L n bar, g)-spaces
-
Tornøe, Christoffer W; Overgaard, Rune V; Agersø, Henrik; Nielsen, Henrik A; Madsen, Henrik; Jonsson, E Niclas
2005-08-01
The objective of the present analysis was to explore the use of stochastic differential equations (SDEs) in population pharmacokinetic/pharmacodynamic (PK/PD) modeling. The intra-individual variability in nonlinear mixed-effects models based on SDEs is decomposed into two types of noise: a measurement and a system noise term. The measurement noise represents uncorrelated error due to, for example, assay error while the system noise accounts for structural misspecifications, approximations of the dynamical model, and true random physiological fluctuations. Since the system noise accounts for model misspecifications, the SDEs provide a diagnostic tool for model appropriateness. The focus of the article is on the implementation of the Extended Kalman Filter (EKF) in NONMEM for parameter estimation in SDE models. Various applications of SDEs in population PK/PD modeling are illustrated through a systematic model development example using clinical PK data of the gonadotropin releasing hormone (GnRH) antagonist degarelix. The dynamic noise estimates were used to track variations in model parameters and systematically build an absorption model for subcutaneously administered degarelix. The EKF-based algorithm was successfully implemented in NONMEM for parameter estimation in population PK/PD models described by systems of SDEs. The example indicated that it was possible to pinpoint structural model deficiencies, and that valuable information may be obtained by tracking unexplained variations in parameters.
-
On Robust Stability of Differential-Algebraic Equations with Structured Uncertainty
Directory of Open Access Journals (Sweden)
A. Kononov
2018-03-01
Full Text Available We consider a linear time-invariant system of differential-algebraic equations (DAE, which can be written as a system of ordinary differential equations with non-invertible coefficients matrices. An important characteristic of DAE is the unsolvability index, which reflects the complexity of the internal structure of the system. The question of the asymptotic stability of DAE containing the uncertainty given by the matrix norm is investigated. We consider a perturbation in the structured uncertainty case. It is assumed that the initial nominal system is asymptotically stable. For the analysis, the original equation is reduced to the structural form, in which the differential and algebraic subsystems are separated. This structural form is equivalent to the input system in the sense of coincidence of sets of solutions, and the operator transforming the DAE into the structural form possesses the inverse operator. The conversion to structural form does not use a change of variables. Regularity of matrix pencil of the source equation is the necessary and sufficient condition of structural form existence. Sufficient conditions have been obtained that perturbations do not break the internal structure of the nominal system. Under these conditions robust stability of the DAE with structured uncertainty is investigated. Estimates for the stability radius of the perturbed DAE system are obtained. The text of the article is from the simpler case, in which the perturbation is present only for an unknown function, to a more complex one, under which the perturbation is also present in the derivative of the unknown function. We used values of the real and the complex stability radii of explicit ordinary differential equations for obtaining the results. We consider the example illustrating the obtained results.
-
An implementation of Kovacic's algorithm for solving ordinary differential equations in FORMAC
International Nuclear Information System (INIS)
Zharkov, A.Yu.
1987-01-01
An implementation of Kovacic's algorithm for finding Liouvillian solutions of the differential equations y'' + a(x)y' + b(x)y = 0 with rational coefficients a(x) and b(x) in the Computer Algebra System FORMAC is described. The algorithm description is presented in such a way that one can easily implement it in a suitable Computer Algebra System
-
WKB: an interactive code for solving differential equations using phase integral methods
International Nuclear Information System (INIS)
White, R.B.
1978-01-01
A small code for the analysis of ordinary differential equations interactively through the use of Phase Integral Methods (WKB) has been written for use on the DEC 10. This note is a descriptive manual for those interested in using the code
-
Directory of Open Access Journals (Sweden)
Resat Yilmazer
2016-02-01
Full Text Available In this work; we present a method for solving the second-order linear ordinary differential equation of hypergeometric type. The solutions of this equation are given by the confluent hypergeometric functions (CHFs. Unlike previous studies, we obtain some different new solutions of the equation without using the CHFs. Therefore, we obtain new discrete fractional solutions of the homogeneous and non-homogeneous confluent hypergeometric differential equation (CHE by using a discrete fractional Nabla calculus operator. Thus, we obtain four different new discrete complex fractional solutions for these equations.
-
Stability with respect to initial time difference for generalized delay differential equations
Directory of Open Access Journals (Sweden)
Ravi Agarwal
2015-02-01
Full Text Available Stability with initial data difference for nonlinear delay differential equations is introduced. This type of stability generalizes the known concept of stability in the literature. It gives us the opportunity to compare the behavior of two nonzero solutions when both initial values and initial intervals are different. Several sufficient conditions for stability and for asymptotic stability with initial time difference are obtained. Lyapunov functions as well as comparison results for scalar ordinary differential equations are employed. Several examples are given to illustrate the theory.
-
International Conference on Differential Equations and Mathematical Physics
Saitō, Yoshimi
1987-01-01
The meeting in Birmingham, Alabama, provided a forum for the discussion of recent developments in the theory of ordinary and partial differential equations, both linear and non-linear, with particular reference to work relating to the equations of mathematical physics. The meeting was attended by about 250 mathematicians from 22 countries. The papers in this volume all involve new research material, with at least outline proofs; some papers also contain survey material. Topics covered include: Schrödinger theory, scattering and inverse scattering, fluid mechanics (including conservative systems and inertial manifold theory attractors), elasticity, non-linear waves, and feedback control theory.
-
Painleve test and discrete Boltzmann equations
International Nuclear Information System (INIS)
Euler, N.; Steeb, W.H.
1989-01-01
The Painleve test for various discrete Boltzmann equations is performed. The connection with integrability is discussed. Furthermore the Lie symmetry vector fields are derived and group-theoretical reduction of the discrete Boltzmann equations to ordinary differentiable equations is performed. Lie Backlund transformations are gained by performing the Painleve analysis for the ordinary differential equations. 16 refs
-
International Nuclear Information System (INIS)
Fakhar, K.; Kara, A. H.
2012-01-01
We study the symmetries, conservation laws and reduction of third-order equations that evolve from a prior reduction of models that arise in fluid phenomena. These could be the ordinary differential equations (ODEs) that are reductions of partial differential equations (PDEs) or, alternatively, PDEs related to given ODEs. In this class, the analysis includes the well-known Blasius, Chazy, and other associated third-order ODEs. (general)
-
Operational method of solution of linear non-integer ordinary and partial differential equations.
Zhukovsky, K V
2016-01-01
We propose operational method with recourse to generalized forms of orthogonal polynomials for solution of a variety of differential equations of mathematical physics. Operational definitions of generalized families of orthogonal polynomials are used in this context. Integral transforms and the operational exponent together with some special functions are also employed in the solutions. The examples of solution of physical problems, related to such problems as the heat propagation in various models, evolutional processes, Black-Scholes-like equations etc. are demonstrated by the operational technique.
-
DEFF Research Database (Denmark)
Mejlbro, Leif
1997-01-01
An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians.......An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians....
-
Institute of Scientific and Technical Information of China (English)
2008-01-01
Using functional derivative technique in quantum field theory, the algebraic dy-namics approach for solution of ordinary differential evolution equations was gen-eralized to treat partial differential evolution equations. The partial differential evo-lution equations were lifted to the corresponding functional partial differential equations in functional space by introducing the time translation operator. The functional partial differential evolution equations were solved by algebraic dynam-ics. The algebraic dynamics solutions are analytical in Taylor series in terms of both initial functions and time. Based on the exact analytical solutions, a new nu-merical algorithm—algebraic dynamics algorithm was proposed for partial differ-ential evolution equations. The difficulty of and the way out for the algorithm were discussed. The application of the approach to and computer numerical experi-ments on the nonlinear Burgers equation and meteorological advection equation indicate that the algebraic dynamics approach and algebraic dynamics algorithm are effective to the solution of nonlinear partial differential evolution equations both analytically and numerically.
-
Delay differential equations for mode-locked semiconductor lasers.
Vladimirov, Andrei G; Turaev, Dmitry; Kozyreff, Gregory
2004-06-01
We propose a new model for passive mode locking that is a set of ordinary delay differential equations. We assume a ring-cavity geometry and Lorentzian spectral filtering of the pulses but do not use small gain and loss and weak saturation approximations. By means of a continuation method, we study mode-locking solutions and their stability. We find that stable mode locking can exist even when the nonlasing state between pulses becomes unstable.
-
Analytic, Algebraic and Geometric Aspects of Differential Equations
Haraoka, Yoshishige; Michalik, Sławomir
2017-01-01
This volume consists of invited lecture notes, survey papers and original research papers from the AAGADE school and conference held in Będlewo, Poland in September 2015. The contributions provide an overview of the current level of interaction between algebra, geometry and analysis and demonstrate the manifold aspects of the theory of ordinary and partial differential equations, while also pointing out the highly fruitful interrelations between those aspects. These interactions continue to yield new developments, not only in the theory of differential equations but also in several related areas of mathematics and physics such as differential geometry, representation theory, number theory and mathematical physics. The main goal of the volume is to introduce basic concepts, techniques, detailed and illustrative examples and theorems (in a manner suitable for non-specialists), and to present recent developments in the field, together with open problems for more advanced and experienced readers. It will be of i...
-
International Nuclear Information System (INIS)
Supriyono; Miyoshi, T.
1997-01-01
NECD Method and runge-Kutta method for large system of second order ordinary differential equations in comparing algorithm. The paper introduce a extrapolation method used for solving the large system of second order ordinary differential equation. We call this method the modified extrapolated central difference (MECD) method. for the accuracy and efficiency MECD method. we compare the method with 4-th order runge-Kutta method. The comparison results show that, this method has almost the same accuracy as the 4-th order runge-Kutta method, but the computation time is about half of runge-Kutta. The MECD was declare by the author and Tetsuhiko Miyoshi of the Dept. Applied Science Yamaguchi University Japan
-
Numerical methods for stochastic partial differential equations with white noise
Zhang, Zhongqiang
2017-01-01
This book covers numerical methods for stochastic partial differential equations with white noise using the framework of Wong-Zakai approximation. The book begins with some motivational and background material in the introductory chapters and is divided into three parts. Part I covers numerical stochastic ordinary differential equations. Here the authors start with numerical methods for SDEs with delay using the Wong-Zakai approximation and finite difference in time. Part II covers temporal white noise. Here the authors consider SPDEs as PDEs driven by white noise, where discretization of white noise (Brownian motion) leads to PDEs with smooth noise, which can then be treated by numerical methods for PDEs. In this part, recursive algorithms based on Wiener chaos expansion and stochastic collocation methods are presented for linear stochastic advection-diffusion-reaction equations. In addition, stochastic Euler equations are exploited as an application of stochastic collocation methods, where a numerical compa...
-
International Nuclear Information System (INIS)
Mokhtari, R.; Toodar, A. Samadi; Chegini, N. G.
2011-01-01
We the extend application of the generalized differential quadrature method (GDQM) to solve some coupled nonlinear Schrödinger equations. The cosine-based GDQM is employed and the obtained system of ordinary differential equations is solved via the fourth order Runge—Kutta method. The numerical solutions coincide with the exact solutions in desired machine precision and invariant quantities are conserved sensibly. Some comparisons with the methods applied in the literature are carried out. (general)
-
Periodic Solutions of a System of Delay Differential Equations for a Small Delay
Directory of Open Access Journals (Sweden)
Adu A.M. Wasike
2002-06-01
Full Text Available We prove the existence of an asymptotically stable periodic solution of a system of delay differential equations with a small time delay t > 0. To achieve this, we transform the system of equations into a system of perturbed ordinary differential equations and then use perturbation results to show the existence of an asymptotically stable periodic solution. This approach is contingent on the fact that the system of equations with t = 0 has a stable limit cycle. We also provide a comparative study of the solutions of the original system and the perturbed system. This comparison lays the ground for proving the existence of periodic solutions of the original system by Schauder's fixed point theorem.
-
Rusyaman, E.; Parmikanti, K.; Chaerani, D.; Asefan; Irianingsih, I.
2018-03-01
One of the application of fractional ordinary differential equation is related to the viscoelasticity, i.e., a correlation between the viscosity of fluids and the elasticity of solids. If the solution function develops into function with two or more variables, then its differential equation must be changed into fractional partial differential equation. As the preliminary study for two variables viscoelasticity problem, this paper discusses about convergence analysis of function sequence which is the solution of the homogenous fractional partial differential equation. The method used to solve the problem is Homotopy Analysis Method. The results show that if given two real number sequences (αn) and (βn) which converge to α and β respectively, then the solution function sequences of fractional partial differential equation with order (αn, βn) will also converge to the solution function of fractional partial differential equation with order (α, β).
-
Fractal-Based Methods and Inverse Problems for Differential Equations: Current State of the Art
Directory of Open Access Journals (Sweden)
Herb E. Kunze
2014-01-01
Full Text Available We illustrate, in this short survey, the current state of the art of fractal-based techniques and their application to the solution of inverse problems for ordinary and partial differential equations. We review several methods based on the Collage Theorem and its extensions. We also discuss two innovative applications: the first one is related to a vibrating string model while the second one considers a collage-based approach for solving inverse problems for partial differential equations on a perforated domain.
-
Directory of Open Access Journals (Sweden)
Tarikul Islam
2018-03-01
Full Text Available In this article, the analytical solutions to the space-time fractional foam drainage equation and the space-time fractional symmetric regularized long wave (SRLW equation are successfully examined by the recently established rational (G′/G-expansion method. The suggested equations are reduced into the nonlinear ordinary differential equations with the aid of the fractional complex transform. Consequently, the theories of the ordinary differential equations are implemented effectively. Three types closed form traveling wave solutions, such as hyperbolic function, trigonometric function and rational, are constructed by using the suggested method in the sense of conformable fractional derivative. The obtained solutions might be significant to analyze the depth and spacing of parallel subsurface drain and small-amplitude long wave on the surface of the water in a channel. It is observed that the performance of the rational (G′/G-expansion method is reliable and will be used to establish new general closed form solutions for any other NPDEs of fractional order.
-
Directory of Open Access Journals (Sweden)
Ahmet Bekir
2014-09-01
Full Text Available In this paper, the fractional partial differential equations are defined by modified Riemann–Liouville fractional derivative. With the help of fractional derivative and traveling wave transformation, these equations can be converted into the nonlinear nonfractional ordinary differential equations. Then G′G-expansion method is applied to obtain exact solutions of the space-time fractional Burgers equation, the space-time fractional KdV-Burgers equation and the space-time fractional coupled Burgers’ equations. As a result, many exact solutions are obtained including hyperbolic function solutions, trigonometric function solutions and rational solutions. These results reveal that the proposed method is very effective and simple in performing a solution to the fractional partial differential equation.
-
Differential Equations Compatible with KZ Equations
International Nuclear Information System (INIS)
Felder, G.; Markov, Y.; Tarasov, V.; Varchenko, A.
2000-01-01
We define a system of 'dynamical' differential equations compatible with the KZ differential equations. The KZ differential equations are associated to a complex simple Lie algebra g. These are equations on a function of n complex variables z i taking values in the tensor product of n finite dimensional g-modules. The KZ equations depend on the 'dual' variable in the Cartan subalgebra of g. The dynamical differential equations are differential equations with respect to the dual variable. We prove that the standard hypergeometric solutions of the KZ equations also satisfy the dynamical equations. As an application we give a new determinant formula for the coordinates of a basis of hypergeometric solutions
-
Hartman-Wintner growth results for sublinear functional differential equations
Directory of Open Access Journals (Sweden)
John A. D. Appleby
2017-01-01
Full Text Available This article determines the rate of growth to infinity of scalar autonomous nonlinear functional and Volterra differential equations. In these equations, the right-hand side is a positive continuous linear functional of f(x. We assume f grows sublinearly, leading to subexponential growth in the solutions. The main results show that the solution of the functional differential equations are asymptotic to that of an auxiliary autonomous ordinary differential equation with right-hand side proportional to f. This happens provided f grows more slowly than l(x=x/log(x. The linear-logarithmic growth rate is also shown to be critical: if f grows more rapidly than l, the ODE dominates the FDE; if f is asymptotic to a constant multiple of l, the FDE and ODE grow at the same rate, modulo a constant non-unit factor; if f grows more slowly than l, the ODE and FDE grow at exactly the same rate. A partial converse of the last result is also proven. In the case when the growth rate is slower than that of the ODE, sharp bounds on the growth rate are determined. The Volterra and finite memory equations can have differing asymptotic behaviour and we explore the source of these differences.
-
Dielectric metasurfaces solve differential and integro-differential equations.
Abdollahramezani, Sajjad; Chizari, Ata; Dorche, Ali Eshaghian; Jamali, Mohammad Vahid; Salehi, Jawad A
2017-04-01
Leveraging subwavelength resonant nanostructures, plasmonic metasurfaces have recently attracted much attention as a breakthrough concept for engineering optical waves both spatially and spectrally. However, inherent ohmic losses concomitant with low coupling efficiencies pose fundamental impediments over their practical applications. Not only can all-dielectric metasurfaces tackle such substantial drawbacks, but also their CMOS-compatible configurations support both Mie resonances that are invariant to the incident angle. Here, we report on a transmittive metasurface comprising arrayed silicon nanodisks embedded in a homogeneous dielectric medium to manipulate phase and amplitude of incident light locally and almost independently. By taking advantage of the interplay between the electric/magnetic resonances and employing general concepts of spatial Fourier transformation, a highly efficient metadevice is proposed to perform mathematical operations including solution of ordinary differential and integro-differential equations with constant coefficients. Our findings further substantiate dielectric metasurfaces as promising candidates for miniaturized, two-dimensional, and planar optical analog computing systems that are much thinner than their conventional lens-based counterparts.
-
Tisdell, Christopher C.
2017-01-01
Knowing an equation has a unique solution is important from both a modelling and theoretical point of view. For over 70 years, the approach to learning and teaching "well posedness" of initial value problems (IVPs) for second- and higher-order ordinary differential equations has involved transforming the problem and its analysis to a…
-
International Nuclear Information System (INIS)
Jimenez, J.C.
2009-06-01
Local Linearization (LL) methods conform a class of one-step explicit integrators for ODEs derived from the following primary and common strategy: the vector field of the differential equation is locally (piecewise) approximated through a first-order Taylor expansion at each time step, thus obtaining successive linear equations that are explicitly integrated. Hereafter, the LL approach may include some additional strategies to improve that basic affine approximation. Theoretical and practical results have shown that the LL integrators have a number of convenient properties. These include arbitrary order of convergence, A-stability, linearization preserving, regularity under quite general conditions, preservation of the dynamics of the exact solution around hyperbolic equilibrium points and periodic orbits, integration of stiff and high-dimensional equations, low computational cost, and others. In this paper, a review of the LL methods and their properties is presented. (author)
-
Soliton-like solutions to the ordinary Schroedinger equation
International Nuclear Information System (INIS)
Zamboni-Rached, Michel; Recami, Erasmo
2011-01-01
In recent times it has been paid attention to the fact that (linear) wave equations admit of soliton-like solutions, known as Localized Waves or Non-diffracting Waves, which propagate without distortion in one direction. Such Localized Solutions (existing also for K-G or Dirac equations) are a priori suitable, more than Gaussian's, for describing elementary particle motion. In this paper we show that, mutatis mutandis, Localized Solutions exist even for the ordinary Schroedinger equation within standard Quantum Mechanics; and we obtain both approximate and exact solutions, also setting forth for them particular examples. In the ideal case such solutions bear infinite energy, as well as plane or spherical waves: we show therefore how to obtain nite-energy solutions. At last, we briefly consider solutions for a particle moving in the presence of a potential. (author)
-
Soliton-like solutions to the ordinary Schroedinger equation
Energy Technology Data Exchange (ETDEWEB)
Zamboni-Rached, Michel [Universidade Estadual de Campinas (DMO/FEEC/UNICAMP), Campinas, SP (Brazil). Fac. de Engenharia Eletrica e de Computacao. Dept. de Microondas e Optica; Recami, Erasmo, E-mail: recami@mi.infn.i [Universita Statale di Bergamo, Bergamo (Italy). Facolta di Ingegneria
2011-07-01
In recent times it has been paid attention to the fact that (linear) wave equations admit of soliton-like solutions, known as Localized Waves or Non-diffracting Waves, which propagate without distortion in one direction. Such Localized Solutions (existing also for K-G or Dirac equations) are a priori suitable, more than Gaussian's, for describing elementary particle motion. In this paper we show that, mutatis mutandis, Localized Solutions exist even for the ordinary Schroedinger equation within standard Quantum Mechanics; and we obtain both approximate and exact solutions, also setting forth for them particular examples. In the ideal case such solutions bear infinite energy, as well as plane or spherical waves: we show therefore how to obtain nite-energy solutions. At last, we briefly consider solutions for a particle moving in the presence of a potential. (author)
-
Sustainability in a Differential Equations Course: A Case Study of Easter Island
Koss, Lorelei
2011-01-01
Easter Island is a fascinating example of resource depletion and population collapse, and its relatively short period of human habitation combined with its isolation lends itself well to investigation by students in a first-semester ordinary differential equations course. This article describes curricular materials for a semester-long case study…
-
International Nuclear Information System (INIS)
Lima, M.L.; Mignaco, J.A.
1983-01-01
The power law potentials in the Schroedinger equation solved recently are shown to come from the classical treatment of the singularities of a linear, second order differential equation. This allows to enlarge the class of solvable power law potentials. (Author) [pt
-
Singular Initial Value Problem for Certain Classes of Systems of Ordinary Differential Equations
Directory of Open Access Journals (Sweden)
Josef Diblík
2013-01-01
dimension of the set of initial data generating such solutions is estimated. An asymptotic behavior of solutions is determined as well and relevant asymptotic formulas are derived. The method of functions defined implicitly and the topological method (Ważewski's method are used in the proofs. The results generalize some previous ones on singular initial value problems for differential equations.
-
Analytic solutions of hydrodynamics equations
International Nuclear Information System (INIS)
Coggeshall, S.V.
1991-01-01
Many similarity solutions have been found for the equations of one-dimensional (1-D) hydrodynamics. These special combinations of variables allow the partial differential equations to be reduced to ordinary differential equations, which must then be solved to determine the physical solutions. Usually, these reduced ordinary differential equations are solved numerically. In some cases it is possible to solve these reduced equations analytically to obtain explicit solutions. In this work a collection of analytic solutions of the 1-D hydrodynamics equations is presented. These can be used for a variety of purposes, including (i) numerical benchmark problems, (ii) as a basis for analytic models, and (iii) to provide insight into more complicated solutions
-
Gaik Tay, Kim; Cheong, Tau Han; Foong Lee, Ming; Kek, Sie Long; Abdul-Kahar, Rosmila
2017-08-01
In the previous work on Euler’s spreadsheet calculator for solving an ordinary differential equation, the Visual Basic for Application (VBA) programming was used, however, a graphical user interface was not developed to capture users input. This weakness may make users confuse on the input and output since those input and output are displayed in the same worksheet. Besides, the existing Euler’s spreadsheet calculator is not interactive as there is no prompt message if there is a mistake in inputting the parameters. On top of that, there are no users’ instructions to guide users to input the derivative function. Hence, in this paper, we improved previous limitations by developing a user-friendly and interactive graphical user interface. This improvement is aimed to capture users’ input with users’ instructions and interactive prompt error messages by using VBA programming. This Euler’s graphical user interface spreadsheet calculator is not acted as a black box as users can click on any cells in the worksheet to see the formula used to implement the numerical scheme. In this way, it could enhance self-learning and life-long learning in implementing the numerical scheme in a spreadsheet and later in any programming language.
-
Numerical Solution of Fuzzy Differential Equations with Z-numbers Using Bernstein Neural Networks
Directory of Open Access Journals (Sweden)
Raheleh Jafari
2017-01-01
Full Text Available The uncertain nonlinear systems can be modeled with fuzzy equations or fuzzy differential equations (FDEs by incorporating the fuzzy set theory. The solutions of them are applied to analyze many engineering problems. However, it is very difficult to obtain solutions of FDEs. In this paper, the solutions of FDEs are approximated by two types of Bernstein neural networks. Here, the uncertainties are in the sense of Z-numbers. Initially the FDE is transformed into four ordinary differential equations (ODEs with Hukuhara differentiability. Then neural models are constructed with the structure of ODEs. With modified back propagation method for Z- number variables, the neural networks are trained. The theory analysis and simulation results show that these new models, Bernstein neural networks, are effective to estimate the solutions of FDEs based on Z-numbers.
-
On symmetries and exact solutions of the Einstein–Maxwell field equations via the symmetry approach
International Nuclear Information System (INIS)
Kaur, Lakhveer; Gupta, R K
2013-01-01
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. (paper)
-
Filimonov, M. Yu.
2017-12-01
The method of special series with recursively calculated coefficients is used to solve nonlinear partial differential equations. The recurrence of finding the coefficients of the series is achieved due to a special choice of functions, in powers of which the solution is expanded in a series. We obtain a sequence of linear partial differential equations to find the coefficients of the series constructed. In many cases, one can deal with a sequence of linear ordinary differential equations. We construct classes of solutions in the form of convergent series for a certain class of nonlinear evolution equations. A new class of solutions of generalized Boussinesque equation with an arbitrary function in the form of a convergent series is constructed.
-
Energy Technology Data Exchange (ETDEWEB)
Lorenzini, R.; Passoni, L. [ENEA, Centro Ricerche Ezio Clementel, Bologna (Italy). Dipt. Ambiente
1999-07-01
The integration of ordinary differential equations systems (ODEs) is of significant concern to tropospheric and stratospheric chemistry modelers. The solution of the ODEs requires a large computational effort because of their stiff nature; in a three-dimensional photochemical model the solution of the ODEs required at least 70% of the total CPU time. Several numerical integration techniques exist which attempt to provide accurate and computationally efficient solutions. In this work it is presented a comparison of some of the techniques in terms of solution accuracy and required computational time. It has been compared the Hybrid Solver (Young and Boris, 1977), the Quasi Steady-State Approximation method (Hesstvedt et al., 1978) and the Chemical Solver for Ordinary Differential Equations (Aro, 1996), by using the CALGRID photochemical model. The accuracy is evaluated by comparing the results of every method with the solutions obtained by the Livermore Solver for Ordinary Differential Equations (Hindmarsh, 1980). The comparison has been made varing the parameters of the error tolerances, and taking into account the trade-off between solution accuracy and computational efficiency. [Italian] L'integrazione di sistemi di equazioni differenziali ordinarie (ODEs), e' un problema significativo per i modellisti della chimica troposferica e stratosferica. A causa della loro natura stiff la soluzione degli ODEs richiese un notevole sforzo computazionale; in un modello fotochimico tridimensionale la soluzione degli ODEs richiede almeno il 70% del tempo totale di CPU. Esistono diverse tecniche di integrazione numerica che possono fornire soluzioni accurate e computazionalmente efficienti: in questo lavoro presentiamo un confronto fra alcune tecniche in termini di accuratezza della soluzione e tempo computazionale richiesto. Si sono confrontati il Solver Ibrido (Young and Boris, 1977), il metodo Quasi Steady-State Approximation (Hesstvedt et al., 1978) ed il Chemical
-
Partial Differential Equations
1988-01-01
The volume contains a selection of papers presented at the 7th Symposium on differential geometry and differential equations (DD7) held at the Nankai Institute of Mathematics, Tianjin, China, in 1986. Most of the contributions are original research papers on topics including elliptic equations, hyperbolic equations, evolution equations, non-linear equations from differential geometry and mechanics, micro-local analysis.
-
Asymptotic behavior of monodromy singularly perturbed differential equations on a Riemann surface
Simpson, Carlos
1991-01-01
This book concerns the question of how the solution of a system of ODE's varies when the differential equation varies. The goal is to give nonzero asymptotic expansions for the solution in terms of a parameter expressing how some coefficients go to infinity. A particular classof families of equations is considered, where the answer exhibits a new kind of behavior not seen in most work known until now. The techniques include Laplace transform and the method of stationary phase, and a combinatorial technique for estimating the contributions of terms in an infinite series expansion for the solution. Addressed primarily to researchers inalgebraic geometry, ordinary differential equations and complex analysis, the book will also be of interest to applied mathematicians working on asymptotics of singular perturbations and numerical solution of ODE's.
-
Datta, D P
2003-01-01
A new class of finitely differentiable scale free solutions to the simplest class of ordinary differential equations is presented. Consequently, the real number set gets replaced by an extended physical set, each element of which is endowed with an equivalence class of infinitesimally separated neighbours in the form of random fluctuations. We show how a sense of time and evolution is intrinsically defined by the infinite continued fraction of the golden mean irrational number (Radical radicand 5 -1)/2, which plays a key role in this extended SL(2,R) formalism of calculus analogous to El Naschie's theory of E sup ( supinfinity sup ) spacetime manifold. Time may thereby undergo random inversions generating well defined random scales, thus allowing a dynamical system to evolve self similarly over the set of multiple scales. The late time stochastic fluctuations of a dynamical system enjoys the generic 1/f spectrum. A universal form of the related probability density is also derived. We prove that the golden mea...
-
Differential equations, dynamical systems, and an introduction to chaos
Smale, Stephen; Devaney, Robert L
2003-01-01
Thirty years in the making, this revised text by three of the world''s leading mathematicians covers the dynamical aspects of ordinary differential equations. it explores the relations between dynamical systems and certain fields outside pure mathematics, and has become the standard textbook for graduate courses in this area. The Second Edition now brings students to the brink of contemporary research, starting from a background that includes only calculus and elementary linear algebra.The authors are tops in the field of advanced mathematics, including Steve Smale who is a recipient of the Field''s Medal for his work in dynamical systems.* Developed by award-winning researchers and authors* Provides a rigorous yet accessible introduction to differential equations and dynamical systems* Includes bifurcation theory throughout* Contains numerous explorations for students to embark uponNEW IN THIS EDITION* New contemporary material and updated applications* Revisions throughout the text, including simplification...
-
Directory of Open Access Journals (Sweden)
Musa Danjuma SHEHU
2008-06-01
Full Text Available This paper lays emphasis on formulation of two dimensional differential games via optimal control theory and consideration of control systems whose dynamics is described by a system of Ordinary Differential equation in the form of linear equation under the influence of two controls U(. and V(.. Base on this, strategies were constructed. Hence we determine the optimal strategy for a control say U(. under a perturbation generated by the second control V(. within a given manifold M.
-
Fiedler, Anna; Raeth, Sebastian; Theis, Fabian J; Hausser, Angelika; Hasenauer, Jan
2016-08-22
Ordinary differential equation (ODE) models are widely used to describe (bio-)chemical and biological processes. To enhance the predictive power of these models, their unknown parameters are estimated from experimental data. These experimental data are mostly collected in perturbation experiments, in which the processes are pushed out of steady state by applying a stimulus. The information that the initial condition is a steady state of the unperturbed process provides valuable information, as it restricts the dynamics of the process and thereby the parameters. However, implementing steady-state constraints in the optimization often results in convergence problems. In this manuscript, we propose two new methods for solving optimization problems with steady-state constraints. The first method exploits ideas from optimization algorithms on manifolds and introduces a retraction operator, essentially reducing the dimension of the optimization problem. The second method is based on the continuous analogue of the optimization problem. This continuous analogue is an ODE whose equilibrium points are the optima of the constrained optimization problem. This equivalence enables the use of adaptive numerical methods for solving optimization problems with steady-state constraints. Both methods are tailored to the problem structure and exploit the local geometry of the steady-state manifold and its stability properties. A parameterization of the steady-state manifold is not required. The efficiency and reliability of the proposed methods is evaluated using one toy example and two applications. The first application example uses published data while the second uses a novel dataset for Raf/MEK/ERK signaling. The proposed methods demonstrated better convergence properties than state-of-the-art methods employed in systems and computational biology. Furthermore, the average computation time per converged start is significantly lower. In addition to the theoretical results, the
-
Solving Differential Equations in R: Package deSolve
Directory of Open Access Journals (Sweden)
Karline Soetaert
2010-02-01
Full Text Available In this paper we present the R package deSolve to solve initial value problems (IVP written as ordinary differential equations (ODE, differential algebraic equations (DAE of index 0 or 1 and partial differential equations (PDE, the latter solved using the method of lines approach. The differential equations can be represented in R code or as compiled code. In the latter case, R is used as a tool to trigger the integration and post-process the results, which facilitates model development and application, whilst the compiled code significantly increases simulation speed. The methods implemented are efficient, robust, and well documented public-domain Fortran routines. They include four integrators from the ODEPACK package (LSODE, LSODES, LSODA, LSODAR, DVODE and DASPK2.0. In addition, a suite of Runge-Kutta integrators and special-purpose solvers to efficiently integrate 1-, 2- and 3-dimensional partial differential equations are available. The routines solve both stiff and non-stiff systems, and include many options, e.g., to deal in an efficient way with the sparsity of the Jacobian matrix, or finding the root of equations. In this article, our objectives are threefold: (1 to demonstrate the potential of using R for dynamic modeling, (2 to highlight typical uses of the different methods implemented and (3 to compare the performance of models specified in R code and in compiled code for a number of test cases. These comparisons demonstrate that, if the use of loops is avoided, R code can efficiently integrate problems comprising several thousands of state variables. Nevertheless, the same problem may be solved from 2 to more than 50 times faster by using compiled code compared to an implementation using only R code. Still, amongst the benefits of R are a more flexible and interactive implementation, better readability of the code, and access to R’s high-level procedures. deSolve is the successor of package odesolve which will be deprecated in
-
Biala, T A; Jator, S N
2015-01-01
In this article, the boundary value method is applied to solve three dimensional elliptic and hyperbolic partial differential equations. The partial derivatives with respect to two of the spatial variables (y, z) are discretized using finite difference approximations to obtain a large system of ordinary differential equations (ODEs) in the third spatial variable (x). Using interpolation and collocation techniques, a continuous scheme is developed and used to obtain discrete methods which are applied via the Block unification approach to obtain approximations to the resulting large system of ODEs. Several test problems are investigated to elucidate the solution process.
-
International Nuclear Information System (INIS)
Silva O, G.; Garcia G, P.
2004-01-01
In this work we describe the procedure to obtain all the family of third order ordinary differential equations connected by a contact transformation such that in their spaces of solutions is defined a conformal three dimensional Minkowski metric. (Author)
-
Benhammouda, Brahim; Vazquez-Leal, Hector
2016-01-01
This work presents an analytical solution of some nonlinear delay differential equations (DDEs) with variable delays. Such DDEs are difficult to treat numerically and cannot be solved by existing general purpose codes. A new method of steps combined with the differential transform method (DTM) is proposed as a powerful tool to solve these DDEs. This method reduces the DDEs to ordinary differential equations that are then solved by the DTM. Furthermore, we show that the solutions can be improved by Laplace-Padé resummation method. Two examples are presented to show the efficiency of the proposed technique. The main advantage of this technique is that it possesses a simple procedure based on a few straight forward steps and can be combined with any analytical method, other than the DTM, like the homotopy perturbation method.
-
Stability by fixed point theory for functional differential equations
Burton, T A
2006-01-01
This book is the first general introduction to stability of ordinary and functional differential equations by means of fixed point techniques. It contains an extensive collection of new and classical examples worked in detail and presented in an elementary manner. Most of this text relies on three principles: a complete metric space, the contraction mapping principle, and an elementary variation of parameters formula. The material is highly accessible to upper-level undergraduate students in the mathematical sciences, as well as working biologists, chemists, economists, engineers, mathematicia
-
A method for solving the KDV equation and some numerical experiments
International Nuclear Information System (INIS)
Chang Jinjiang.
1993-01-01
In this paper, by means of difference method for discretization of space partial derivatives of KDV equation, an initial value problem in ordinary differential equations of large dimensions is produced. By using this ordinary differential equations the existence and the uniqueness of the solution of the KDV equation and the conservation of scheme are proved. This ordinary differential equation can be solved by using implicit Runge-Kutta methods, so a new method for finding the numerical solution of the KDV equation is presented. Numerical experiments not only describe in detail the procedure of two solitons collision, soliton reflex and soliton produce, but also show that this method is very effective. (author). 7 refs, 3 figs
-
Differential equations with Matlab exploration, applications, and theory
McKibben, Mark
2014-01-01
ORDINARY DIFFERENTIAL EQUATIONS Welcome! Introduction This Book Is a Field Guide. What Does That Mean for YOU? Mired in Jargon - A Quick Language Lesson! Introducing MATLAB A First Look at Some Elementary Mathematical Models A Basic Analysis Toolbox Some Basic Mathematical Shorthand Set Algebra Functions The Space (R; j_j) A Closer Look at Sequences in (R; j_j) The Spaces (RN; k_kRN ) and (MN(R); k_kMN(R)Calculus of RN-valued and MN(R)-valued FunctionsSome Elementary ODEs Looking Ahead A First Wave of Mathematical Models Newton's Law of Heating and Cooling-Revisited Pharmocokinetics Uniform Mi
-
International Nuclear Information System (INIS)
Winkler, E.
1991-01-01
The general theory of inhomogeneous compartments with age-dependent elimination rates is illustrated by examples. Mathematically, it turns out that models consisting of partial differential equations include ordinary, delayed and integro-differential equations, a general fact which is treated here in the context of linear tracer kinetics. The examples include standard compartments as a degenerate case, systems of standard compartments (compartment blocks), models resulting in special residence time distributions, models with pipes, and systems with heterogeneous particles. (orig./BBR) [de
-
Limiting precision in differential equation solvers. II Sources of trouble and starting a code
International Nuclear Information System (INIS)
Shampine, L.F.
1978-01-01
The reasons a class of codes for solving ordinary differential equations might want to use an extremely small step size are investigated. For this class the likelihood of precision difficulties is evaluated and remedies examined. The investigations suggests a way of selecting automatically an initial step size which should be reliably on scale
-
International Nuclear Information System (INIS)
Donolato, C
2005-01-01
A relation between the gradient of the electric field and mean curvature of equipotential surfaces (Green's differential equation) is applied to a two-dimensional free-boundary problem arising in the study of ion sheaths around wedge-shaped cathodes. With the assumption that the equipotential lines are hyperbolae, this relation leads to a nonlinear ordinary differential equation for the potential along the bisector line of the wedge. An approximate solution is found, which yields, in particular, the sheath width along this line as a function of the wedge angle. The resulting values are in good agreement with published results obtained by numerically solving Poisson's equation
-
Similarity Solutions for Multiterm Time-Fractional Diffusion Equation
Elsaid, A.; Abdel Latif, M. S.; Maneea, M.
2016-01-01
Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on ...
-
Tisdell, Christopher C.
2017-07-01
Knowing an equation has a unique solution is important from both a modelling and theoretical point of view. For over 70 years, the approach to learning and teaching 'well posedness' of initial value problems (IVPs) for second- and higher-order ordinary differential equations has involved transforming the problem and its analysis to a first-order system of equations. We show that this excursion is unnecessary and present a direct approach regarding second- and higher-order problems that does not require an understanding of systems.
-
Differential equations for dummies
Holzner, Steven
2008-01-01
The fun and easy way to understand and solve complex equations Many of the fundamental laws of physics, chemistry, biology, and economics can be formulated as differential equations. This plain-English guide explores the many applications of this mathematical tool and shows how differential equations can help us understand the world around us. Differential Equations For Dummies is the perfect companion for a college differential equations course and is an ideal supplemental resource for other calculus classes as well as science and engineering courses. It offers step-by-step techniques, practical tips, numerous exercises, and clear, concise examples to help readers improve their differential equation-solving skills and boost their test scores.
-
Functional equations with causal operators
Corduneanu, C
2003-01-01
Functional equations encompass most of the equations used in applied science and engineering: ordinary differential equations, integral equations of the Volterra type, equations with delayed argument, and integro-differential equations of the Volterra type. The basic theory of functional equations includes functional differential equations with causal operators. Functional Equations with Causal Operators explains the connection between equations with causal operators and the classical types of functional equations encountered by mathematicians and engineers. It details the fundamentals of linear equations and stability theory and provides several applications and examples.
-
Murphy, K. A.
1990-01-01
A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.
-
Differential equations with applications in cancer diseases.
Ilea, M; Turnea, M; Rotariu, M
2013-01-01
Mathematical modeling is a process by which a real world problem is described by a mathematical formulation. The cancer modeling is a highly challenging problem at the frontier of applied mathematics. A variety of modeling strategies have been developed, each focusing on one or more aspects of cancer. The vast majority of mathematical models in cancer diseases biology are formulated in terms of differential equations. We propose an original mathematical model with small parameter for the interactions between these two cancer cell sub-populations and the mathematical model of a vascular tumor. We work on the assumption that, the quiescent cells' nutrient consumption is long. One the equations system includes small parameter epsilon. The smallness of epsilon is relative to the size of the solution domain. MATLAB simulations obtained for transition rate from the quiescent cells' nutrient consumption is long, we show a similar asymptotic behavior for two solutions of the perturbed problem. In this system, the small parameter is an asymptotic variable, different from the independent variable. The graphical output for a mathematical model of a vascular tumor shows the differences in the evolution of the tumor populations of proliferating, quiescent and necrotic cells. The nutrient concentration decreases sharply through the viable rim and tends to a constant level in the core due to the nearly complete necrosis in this region. Many mathematical models can be quantitatively characterized by ordinary differential equations or partial differential equations. The use of MATLAB in this article illustrates the important role of informatics in research in mathematical modeling. The study of avascular tumor growth cells is an exciting and important topic in cancer research and will profit considerably from theoretical input. Interpret these results to be a permanent collaboration between math's and medical oncologists.
-
Directory of Open Access Journals (Sweden)
Sufia Zulfa Ahmad
2016-01-01
Full Text Available We derived a two-step, four-stage, and fifth-order semi-implicit hybrid method which can be used for solving special second-order ordinary differential equations. The method is then trigonometrically fitted so that it is suitable for solving problems which are oscillatory in nature. The methods are then used for solving oscillatory delay differential equations. Numerical results clearly show the efficiency of the new method when compared to the existing explicit and implicit methods in the scientific literature.
-
DEFF Research Database (Denmark)
Overgaard, Rune Viig; Holford, Nick; Rytved, K. A.
2007-01-01
Purpose To describe the pharmacodynamic effects of recombinant human interleukin-21 (IL-21) on core body temperature in cynomolgus monkeys using basic mechanisms of heat regulation. A major effort was devoted to compare the use of ordinary differential equations (ODEs) with stochastic differentia...
-
National Research Council Canada - National Science Library
Mitchell, Jason
2002-01-01
A method is presented for the generation of exact numerical coefficients found in two families of implicit Chebyshev methods for the numerical integration of first- and second-order ordinary differential equations...
-
Partial differential equations
Evans, Lawrence C
2010-01-01
This text gives a comprehensive survey of modern techniques in the theoretical study of partial differential equations (PDEs) with particular emphasis on nonlinear equations. The exposition is divided into three parts: representation formulas for solutions; theory for linear partial differential equations; and theory for nonlinear partial differential equations. Included are complete treatments of the method of characteristics; energy methods within Sobolev spaces; regularity for second-order elliptic, parabolic, and hyperbolic equations; maximum principles; the multidimensional calculus of variations; viscosity solutions of Hamilton-Jacobi equations; shock waves and entropy criteria for conservation laws; and, much more.The author summarizes the relevant mathematics required to understand current research in PDEs, especially nonlinear PDEs. While he has reworked and simplified much of the classical theory (particularly the method of characteristics), he primarily emphasizes the modern interplay between funct...
-
A new auxiliary equation and exact travelling wave solutions of nonlinear equations
International Nuclear Information System (INIS)
Sirendaoreji
2006-01-01
A new auxiliary ordinary differential equation and its solutions are used for constructing exact travelling wave solutions of nonlinear partial differential equations in a unified way. The main idea of this method is to take full advantage of the auxiliary equation which has more new exact solutions. More new exact travelling wave solutions are obtained for the quadratic nonlinear Klein-Gordon equation, the combined KdV and mKdV equation, the sine-Gordon equation and the Whitham-Broer-Kaup equations
-
Alternans promotion in cardiac electrophysiology models by delay differential equations.
Gomes, Johnny M; Dos Santos, Rodrigo Weber; Cherry, Elizabeth M
2017-09-01
Cardiac electrical alternans is a state of alternation between long and short action potentials and is frequently associated with harmful cardiac conditions. Different dynamic mechanisms can give rise to alternans; however, many cardiac models based on ordinary differential equations are not able to reproduce this phenomenon. A previous study showed that alternans can be induced by the introduction of delay differential equations (DDEs) in the formulations of the ion channel gating variables of a canine myocyte model. The present work demonstrates that this technique is not model-specific by successfully promoting alternans using DDEs for five cardiac electrophysiology models that describe different types of myocytes, with varying degrees of complexity. By analyzing results across the different models, we observe two potential requirements for alternans promotion via DDEs for ionic gates: (i) the gate must have a significant influence on the action potential duration and (ii) a delay must significantly impair the gate's recovery between consecutive action potentials.
-
Alternans promotion in cardiac electrophysiology models by delay differential equations
Gomes, Johnny M.; dos Santos, Rodrigo Weber; Cherry, Elizabeth M.
2017-09-01
Cardiac electrical alternans is a state of alternation between long and short action potentials and is frequently associated with harmful cardiac conditions. Different dynamic mechanisms can give rise to alternans; however, many cardiac models based on ordinary differential equations are not able to reproduce this phenomenon. A previous study showed that alternans can be induced by the introduction of delay differential equations (DDEs) in the formulations of the ion channel gating variables of a canine myocyte model. The present work demonstrates that this technique is not model-specific by successfully promoting alternans using DDEs for five cardiac electrophysiology models that describe different types of myocytes, with varying degrees of complexity. By analyzing results across the different models, we observe two potential requirements for alternans promotion via DDEs for ionic gates: (i) the gate must have a significant influence on the action potential duration and (ii) a delay must significantly impair the gate's recovery between consecutive action potentials.
-
Chow, Sy-Miin; Ou, Lu; Ciptadi, Arridhana; Prince, Emily B; You, Dongjun; Hunter, Michael D; Rehg, James M; Rozga, Agata; Messinger, Daniel S
2018-06-01
A growing number of social scientists have turned to differential equations as a tool for capturing the dynamic interdependence among a system of variables. Current tools for fitting differential equation models do not provide a straightforward mechanism for diagnosing evidence for qualitative shifts in dynamics, nor do they provide ways of identifying the timing and possible determinants of such shifts. In this paper, we discuss regime-switching differential equation models, a novel modeling framework for representing abrupt changes in a system of differential equation models. Estimation was performed by combining the Kim filter (Kim and Nelson State-space models with regime switching: classical and Gibbs-sampling approaches with applications, MIT Press, Cambridge, 1999) and a numerical differential equation solver that can handle both ordinary and stochastic differential equations. The proposed approach was motivated by the need to represent discrete shifts in the movement dynamics of [Formula: see text] mother-infant dyads during the Strange Situation Procedure (SSP), a behavioral assessment where the infant is separated from and reunited with the mother twice. We illustrate the utility of a novel regime-switching differential equation model in representing children's tendency to exhibit shifts between the goal of staying close to their mothers and intermittent interest in moving away from their mothers to explore the room during the SSP. Results from empirical model fitting were supplemented with a Monte Carlo simulation study to evaluate the use of information criterion measures to diagnose sudden shifts in dynamics.
-
Interpreting experimental data on egg production--applications of dynamic differential equations.
France, J; Lopez, S; Kebreab, E; Dijkstra, J
2013-09-01
This contribution focuses on applying mathematical models based on systems of ordinary first-order differential equations to synthesize and interpret data from egg production experiments. Models based on linear systems of differential equations are contrasted with those based on nonlinear systems. Regression equations arising from analytical solutions to linear compartmental schemes are considered as candidate functions for describing egg production curves, together with aspects of parameter estimation. Extant candidate functions are reviewed, a role for growth functions such as the Gompertz equation suggested, and a function based on a simple new model outlined. Structurally, the new model comprises a single pool with an inflow and an outflow. Compartmental simulation models based on nonlinear systems of differential equations, and thus requiring numerical solution, are next discussed, and aspects of parameter estimation considered. This type of model is illustrated in relation to development and evaluation of a dynamic model of calcium and phosphorus flows in layers. The model consists of 8 state variables representing calcium and phosphorus pools in the crop, stomachs, plasma, and bone. The flow equations are described by Michaelis-Menten or mass action forms. Experiments that measure Ca and P uptake in layers fed different calcium concentrations during shell-forming days are used to evaluate the model. In addition to providing a useful management tool, such a simulation model also provides a means to evaluate feeding strategies aimed at reducing excretion of potential pollutants in poultry manure to the environment.
-
Singular stochastic differential equations
Cherny, Alexander S
2005-01-01
The authors introduce, in this research monograph on stochastic differential equations, a class of points termed isolated singular points. Stochastic differential equations possessing such points (called singular stochastic differential equations here) arise often in theory and in applications. However, known conditions for the existence and uniqueness of a solution typically fail for such equations. The book concentrates on the study of the existence, the uniqueness, and, what is most important, on the qualitative behaviour of solutions of singular stochastic differential equations. This is done by providing a qualitative classification of isolated singular points, into 48 possible types.
-
Mei, Shu-Li; Lv, Hong-Liang; Ma, Qin
2008-01-01
Based on restricted variational principle, a novel method for interval wavelet construction is proposed. For the excellent local property of quasi-Shannon wavelet, its interval wavelet is constructed, and then applied to solve ordinary differential equations. Parameter choices for the interval wavelet method are discussed and its numerical performance is demonstrated.
-
Block Hybrid Collocation Method with Application to Fourth Order Differential Equations
Directory of Open Access Journals (Sweden)
Lee Ken Yap
2015-01-01
Full Text Available The block hybrid collocation method with three off-step points is proposed for the direct solution of fourth order ordinary differential equations. The interpolation and collocation techniques are applied on basic polynomial to generate the main and additional methods. These methods are implemented in block form to obtain the approximation at seven points simultaneously. Numerical experiments are conducted to illustrate the efficiency of the method. The method is also applied to solve the fourth order problem from ship dynamics.
-
Differential Equation over Banach Algebra
Kleyn, Aleks
2018-01-01
In the book, I considered differential equations of order $1$ over Banach $D$-algebra: differential equation solved with respect to the derivative; exact differential equation; linear homogeneous equation. In noncommutative Banach algebra, initial value problem for linear homogeneous equation has infinitely many solutions.
-
An algorithm for solving initial value problems of third order ordinary ...
African Journals Online (AJOL)
Abstract. We propose an implicit multi-step method for the solution of initial value problems (IVPs) of third order ordinary differential equations (ODE) which does not require reducing the ODE to first order before solving. The development of the method is based on collocation of the differential system and interpolation of the ...
-
Uncertain differential equations
Yao, Kai
2016-01-01
This book introduces readers to the basic concepts of and latest findings in the area of differential equations with uncertain factors. It covers the analytic method and numerical method for solving uncertain differential equations, as well as their applications in the field of finance. Furthermore, the book provides a number of new potential research directions for uncertain differential equation. It will be of interest to researchers, engineers and students in the fields of mathematics, information science, operations research, industrial engineering, computer science, artificial intelligence, automation, economics, and management science.
-
Equations of motion for a rotor blade, including gravity, pitch action and rotor speed variations
DEFF Research Database (Denmark)
Kallesøe, Bjarne Skovmose
2007-01-01
This paper extends Hodges-Dowell's partial differential equations of blade motion, by including the effects from gravity, pitch action and varying rotor speed. New equations describing the pitch action and rotor speeds are also derived. The physical interpretation of the individual terms...... in the equations is discussed. The partial differential equations of motion are approximated by ordinary differential equations of motion using an assumed mode method. The ordinary differential equations are used to simulate a sudden pitch change of a rotating blade. This work is a part of a project on pitch blade...
-
International Nuclear Information System (INIS)
Ignat'yev, A O
2003-01-01
A system of ordinary differential equations with impulse action at fixed moments of time is considered. The system is assumed to have the zero solution. It is shown that the existence of a corresponding Lyapunov function is a necessary and sufficient condition for the uniform asymptotic stability of the zero solution. Restrictions on perturbations of the right-hand sides of differential equations and impulse actions are obtained under which the uniform asymptotic stability of the zero solution of the 'unperturbed' system implies the uniform asymptotic stability of the zero solution of the 'perturbed' system
-
Algorithmic Verification of Linearizability for Ordinary Differential Equations
Lyakhov, Dmitry A.; Gerdt, Vladimir P.; Michels, Dominik L.
2017-01-01
one by a point transformation of the dependent and independent variables. The first algorithm is based on a construction of the Lie point symmetry algebra and on the computation of its derived algebra. The second algorithm exploits the differential
-
A Tutorial on RxODE: Simulating Differential Equation Pharmacometric Models in R.
Wang, W; Hallow, K M; James, D A
2016-01-01
This tutorial presents the application of an R package, RxODE, that facilitates quick, efficient simulations of ordinary differential equation models completely within R. Its application is illustrated through simulation of design decision effects on an adaptive dosing regimen. The package provides an efficient, versatile way to specify dosing scenarios and to perform simulation with variability with minimal custom coding. Models can be directly translated to Rshiny applications to facilitate interactive, real-time evaluation/iteration on simulation scenarios.
-
Computational issues of solving the 1D steady gradually varied flow equation
Directory of Open Access Journals (Sweden)
Artichowicz Wojciech
2014-09-01
Full Text Available In this paper a problem of multiple solutions of steady gradually varied flow equation in the form of the ordinary differential energy equation is discussed from the viewpoint of its numerical solution. Using the Lipschitz theorem dealing with the uniqueness of solution of an initial value problem for the ordinary differential equation it was shown that the steady gradually varied flow equation can have more than one solution. This fact implies that the nonlinear algebraic equation approximating the ordinary differential energy equation, which additionally coincides with the wellknown standard step method usually applied for computing of the flow profile, can have variable number of roots. Consequently, more than one alternative solution corresponding to the same initial condition can be provided. Using this property it is possible to compute the water flow profile passing through the critical stage.
-
Two-scale approach to oscillatory singularly perturbed transport equations
Frénod, Emmanuel
2017-01-01
This book presents the classical results of the two-scale convergence theory and explains – using several figures – why it works. It then shows how to use this theory to homogenize ordinary differential equations with oscillating coefficients as well as oscillatory singularly perturbed ordinary differential equations. In addition, it explores the homogenization of hyperbolic partial differential equations with oscillating coefficients and linear oscillatory singularly perturbed hyperbolic partial differential equations. Further, it introduces readers to the two-scale numerical methods that can be built from the previous approaches to solve oscillatory singularly perturbed transport equations (ODE and hyperbolic PDE) and demonstrates how they can be used efficiently. This book appeals to master’s and PhD students interested in homogenization and numerics, as well as to the Iter community.
-
Generalized differential transform method to differential-difference equation
International Nuclear Information System (INIS)
Zou Li; Wang Zhen; Zong Zhi
2009-01-01
In this Letter, we generalize the differential transform method to solve differential-difference equation for the first time. Two simple but typical examples are applied to illustrate the validity and the great potential of the generalized differential transform method in solving differential-difference equation. A Pade technique is also introduced and combined with GDTM in aim of extending the convergence area of presented series solutions. Comparisons are made between the results of the proposed method and exact solutions. Then we apply the differential transform method to the discrete KdV equation and the discrete mKdV equation, and successfully obtain solitary wave solutions. The results reveal that the proposed method is very effective and simple. We should point out that generalized differential transform method is also easy to be applied to other nonlinear differential-difference equation.
-
Introduction to partial differential equations
Greenspan, Donald
2000-01-01
Designed for use in a one-semester course by seniors and beginning graduate students, this rigorous presentation explores practical methods of solving differential equations, plus the unifying theory underlying the mathematical superstructure. Topics include basic concepts, Fourier series, second-order partial differential equations, wave equation, potential equation, heat equation, approximate solution of partial differential equations, and more. Exercises appear at the ends of most chapters. 1961 edition.
-
Darboux transformation for the NLS equation
International Nuclear Information System (INIS)
Aktosun, Tuncay; Mee, Cornelis van der
2010-01-01
We analyze a certain class of integral equations associated with Marchenko equations and Gel'fand-Levitan equations. Such integral equations arise through a Fourier transformation on various ordinary differential equations involving a spectral parameter. When the integral operator is perturbed by a finite-rank perturbation, we explicitly evaluate the change in the solution in terms of the unperturbed quantities and the finite-rank perturbation. We show that this result provides a fundamental approach to derive Darboux transformations for various systems of ordinary differential operators. We illustrate our theory by providing the explicit Darboux transformation for the Zakharov-Shabat system and show how the potential and wave function change when a simple discrete eigenvalue is added to the spectrum, and thus we also provide a one-parameter family of Darboux transformations for the nonlinear Schroedinger equation.
-
Solution of the reactor point kinetics equations by MATLAB computing
Directory of Open Access Journals (Sweden)
Singh Sudhansu S.
2015-01-01
Full Text Available The numerical solution of the point kinetics equations in the presence of Newtonian temperature feedback has been a challenging issue for analyzing the reactor transients. Reactor point kinetics equations are a system of stiff ordinary differential equations which need special numerical treatments. Although a plethora of numerical intricacies have been introduced to solve the point kinetics equations over the years, some of the simple and straightforward methods still work very efficiently with extraordinary accuracy. As an example, it has been shown recently that the fundamental backward Euler finite difference algorithm with its simplicity has proven to be one of the most effective legacy methods. Complementing the back-ward Euler finite difference scheme, the present work demonstrates the application of ordinary differential equation suite available in the MATLAB software package to solve the stiff reactor point kinetics equations with Newtonian temperature feedback effects very effectively by analyzing various classic benchmark cases. Fair accuracy of the results implies the efficient application of MATLAB ordinary differential equation suite for solving the reactor point kinetics equations as an alternate method for future applications.
-
Similarity Solutions for Multiterm Time-Fractional Diffusion Equation
Directory of Open Access Journals (Sweden)
A. Elsaid
2016-01-01
Full Text Available Similarity method is employed to solve multiterm time-fractional diffusion equation. The orders of the fractional derivatives belong to the interval (0,1] and are defined in the Caputo sense. We illustrate how the problem is reduced from a multiterm two-variable fractional partial differential equation to a multiterm ordinary fractional differential equation. Power series solution is obtained for the resulting ordinary problem and the convergence of the series solution is discussed. Based on the obtained results, we propose a definition for a multiterm error function with generalized coefficients.
-
Symmetry reduction for nonlinear wave equations in Riemannian and pseudo-Riemannian spaces
International Nuclear Information System (INIS)
Grundland, A.M.; Harnad, J.; Winternitz, P.
1984-01-01
The authors show how group theory can be systematically employed to reduce nonlinear partial differential equations in n independent variables to partial differential equations in fewer variables and in particular, to ordinary differential equations. (Auth.)
-
Fractional Order Stochastic Differential Equation with Application in European Option Pricing
Directory of Open Access Journals (Sweden)
Qing Li
2014-01-01
Full Text Available Memory effect is an important phenomenon in financial systems, and a number of research works have been carried out to study the long memory in the financial markets. In recent years, fractional order ordinary differential equation is used as an effective instrument for describing the memory effect in complex systems. In this paper, we establish a fractional order stochastic differential equation (FSDE model to describe the effect of trend memory in financial pricing. We, then, derive a European option pricing formula based on the FSDE model and prove the existence of the trend memory (i.e., the mean value function in the option pricing formula when the Hurst index is between 0.5 and 1. In addition, we make a comparison analysis between our proposed model, the classic Black-Scholes model, and the stochastic model with fractional Brownian motion. Numerical results suggest that our model leads to more accurate and lower standard deviation in the empirical study.
-
Energy Technology Data Exchange (ETDEWEB)
Ho, C.-L. [Department of Physics, Tamkang University, Tamsui 25137, Taiwan (China); Lee, C.-C., E-mail: chieh.no27@gmail.com [Center of General Education, Aletheia University, Tamsui 25103, Taiwan (China)
2016-01-15
We consider solvability of the generalized reaction–diffusion equation with both space- and time-dependent diffusion and reaction terms by means of the similarity method. By introducing the similarity variable, the reaction–diffusion equation is reduced to an ordinary differential equation. Matching the resulting ordinary differential equation with known exactly solvable equations, one can obtain corresponding exactly solvable reaction–diffusion systems. Several representative examples of exactly solvable reaction–diffusion equations are presented.
-
Khan, Faiz M; Schmitz, Ulf; Nikolov, Svetoslav; Engelmann, David; Pützer, Brigitte M; Wolkenhauer, Olaf; Vera, Julio
2014-01-01
A decade of successful results indicates that systems biology is the appropriate approach to investigate the regulation of complex biochemical networks involving transcriptional and post-transcriptional regulations. It becomes mandatory when dealing with highly interconnected biochemical networks, composed of hundreds of compounds, or when networks are enriched in non-linear motifs like feedback and feedforward loops. An emerging dilemma is to conciliate models of massive networks and the adequate description of non-linear dynamics in a suitable modeling framework. Boolean networks are an ideal representation of massive networks that are humble in terms of computational complexity and data demand. However, they are inappropriate when dealing with nested feedback/feedforward loops, structural motifs common in biochemical networks. On the other hand, models of ordinary differential equations (ODEs) cope well with these loops, but they require enormous amounts of quantitative data for a full characterization of the model. Here we propose hybrid models, composed of ODE and logical sub-modules, as a strategy to handle large scale, non-linear biochemical networks that include transcriptional and post-transcriptional regulations. We illustrate the construction of this kind of models using as example a regulatory network centered on E2F1, a transcription factor involved in cancer. The hybrid modeling approach proposed is a good compromise between quantitative/qualitative accuracy and scalability when considering large biochemical networks with a small highly interconnected core, and module of transcriptionally regulated genes that are not part of critical regulatory loops. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications. Guest Editor: Yudong Cai. Copyright © 2013 Elsevier B.V. All rights reserved.
-
A hybrid Pade-Galerkin technique for differential equations
Geer, James F.; Andersen, Carl M.
1993-01-01
A three-step hybrid analysis technique, which successively uses the regular perturbation expansion method, the Pade expansion method, and then a Galerkin approximation, is presented and applied to some model boundary value problems. In the first step of the method, the regular perturbation method is used to construct an approximation to the solution in the form of a finite power series in a small parameter epsilon associated with the problem. In the second step of the method, the series approximation obtained in step one is used to construct a Pade approximation in the form of a rational function in the parameter epsilon. In the third step, the various powers of epsilon which appear in the Pade approximation are replaced by new (unknown) parameters (delta(sub j)). These new parameters are determined by requiring that the residual formed by substituting the new approximation into the governing differential equation is orthogonal to each of the perturbation coordinate functions used in step one. The technique is applied to model problems involving ordinary or partial differential equations. In general, the technique appears to provide good approximations to the solution even when the perturbation and Pade approximations fail to do so. The method is discussed and topics for future investigations are indicated.
-
Automatic differentiation bibliography
Energy Technology Data Exchange (ETDEWEB)
Corliss, G.F. [comp.
1992-07-01
This is a bibliography of work related to automatic differentiation. Automatic differentiation is a technique for the fast, accurate propagation of derivative values using the chain rule. It is neither symbolic nor numeric. Automatic differentiation is a fundamental tool for scientific computation, with applications in optimization, nonlinear equations, nonlinear least squares approximation, stiff ordinary differential equation, partial differential equations, continuation methods, and sensitivity analysis. This report is an updated version of the bibliography which originally appeared in Automatic Differentiation of Algorithms: Theory, Implementation, and Application.
-
Introduction to differential equations
Taylor, Michael E
2011-01-01
The mathematical formulations of problems in physics, economics, biology, and other sciences are usually embodied in differential equations. The analysis of the resulting equations then provides new insight into the original problems. This book describes the tools for performing that analysis. The first chapter treats single differential equations, emphasizing linear and nonlinear first order equations, linear second order equations, and a class of nonlinear second order equations arising from Newton's laws. The first order linear theory starts with a self-contained presentation of the exponen
-
Wu, Jiayang; Cao, Pan; Hu, Xiaofeng; Jiang, Xinhong; Pan, Ting; Yang, Yuxing; Qiu, Ciyuan; Tremblay, Christine; Su, Yikai
2014-10-20
We propose and experimentally demonstrate an all-optical temporal differential-equation solver that can be used to solve ordinary differential equations (ODEs) characterizing general linear time-invariant (LTI) systems. The photonic device implemented by an add-drop microring resonator (MRR) with two tunable interferometric couplers is monolithically integrated on a silicon-on-insulator (SOI) wafer with a compact footprint of ~60 μm × 120 μm. By thermally tuning the phase shifts along the bus arms of the two interferometric couplers, the proposed device is capable of solving first-order ODEs with two variable coefficients. The operation principle is theoretically analyzed, and system testing of solving ODE with tunable coefficients is carried out for 10-Gb/s optical Gaussian-like pulses. The experimental results verify the effectiveness of the fabricated device as a tunable photonic ODE solver.
-
Traveling solitary wave solutions to evolution equations with nonlinear terms of any order
International Nuclear Information System (INIS)
Feng Zhaosheng
2003-01-01
Many physical phenomena in one- or higher-dimensional space can be described by nonlinear evolution equations, which can be reduced to ordinary differential equations such as the Lienard equation. Thus, to study those ordinary differential equations is of significance not only in mathematics itself, but also in physics. In this paper, a kind of explicit exact solutions to the Lienard equation is obtained. The applications of the solutions to the nonlinear RR-equation and the compound KdV-type equation are presented, which extend the results obtained in the previous literature
-
International Nuclear Information System (INIS)
Foroutan, A.
1992-05-01
The essential mathematical challenge in transport theory is based on the nonlinearity of the integro-differential equations governing classical thermodynamic systems on molecular kinetic level. It is the aim of this thesis to gain exact analytical solutions to the model Boltzmann equation suggested by Tjon and Wu. Such solutions afford a deeper insight into the dynamics of rarefied gases. Tjon and Wu have provided a stochastic model of a Boltzmann equation. Its transition probability depends only on the relative speed of the colliding particles. This assumption leads in the case of two translational degrees of freedom to an integro-differential equation of convolution type. According to this convolution structure the integro-differential equation is Laplace transformed. The result is a nonlinear partial differential equation. The investigation of the symmetries of this differential equation by means of Lie groups of transformation enables us to transform the originally nonlinear partial differential equation into ordinary differential equation into ordinary differential equations of Bernoulli type. (author)
-
BOOK REVIEW: Partial Differential Equations in General Relativity
Halburd, Rodney G.
2008-11-01
Although many books on general relativity contain an overview of the relevant background material from differential geometry, very little attention is usually paid to background material from the theory of differential equations. This is understandable in a first course on relativity but it often limits the kinds of problems that can be studied rigorously. Einstein's field equations lie at the heart of general relativity. They are a system of partial differential equations (PDEs) relating the curvature of spacetime to properties of matter. A central part of most problems in general relativity is to extract information about solutions of these equations. Most standard texts achieve this by studying exact solutions or numerical and analytical approximations. In the book under review, Alan Rendall emphasises the role of rigorous qualitative methods in general relativity. There has long been a need for such a book, giving a broad overview of the relevant background from the theory of partial differential equations, and not just from differential geometry. It should be noted that the book also covers the basic theory of ordinary differential equations. Although there are many good books on the rigorous theory of PDEs, methods related to the Einstein equations deserve special attention, not only because of the complexity and importance of these equations, but because these equations do not fit into any of the standard classes of equations (elliptic, parabolic, hyperbolic) that one typically encounters in a course on PDEs. Even specifying exactly what ones means by a Cauchy problem in general relativity requires considerable care. The main problem here is that the manifold on which the solution is defined is determined by the solution itself. This means that one does not simply define data on a submanifold. Rendall's book gives a good overview of applications and results from the qualitative theory of PDEs to general relativity. It would be impossible to give detailed
-
Nonlinear evolution equations and Painlevé test
Steeb, Willi-Hans
1988-01-01
This book is an edited version of lectures given by the authors at a seminar at the Rand Afrikaans University. It gives a survey on the Painlevé test, Painlevé property and integrability. Both ordinary differential equations and partial differential equations are considered.
-
Arithmetic differential equations on $GL_n$, I: differential cocycles
Buium, Alexandru; Dupuy, Taylor
2013-01-01
The theory of differential equations has an arithmetic analogue in which derivatives are replaced by Fermat quotients. One can then ask what is the arithmetic analogue of a linear differential equation. The study of usual linear differential equations is the same as the study of the differential cocycle from $GL_n$ into its Lie algebra given by the logarithmic derivative. However we prove here that there are no such cocycles in the context of arithmetic differential equations. In sequels of t...
-
Modeling of delays in PKPD: classical approaches and a tutorial for delay differential equations.
Koch, Gilbert; Krzyzanski, Wojciech; Pérez-Ruixo, Juan Jose; Schropp, Johannes
2014-08-01
In pharmacokinetics/pharmacodynamics (PKPD) the measured response is often delayed relative to drug administration, individuals in a population have a certain lifespan until they maturate or the change of biomarkers does not immediately affects the primary endpoint. The classical approach in PKPD is to apply transit compartment models (TCM) based on ordinary differential equations to handle such delays. However, an alternative approach to deal with delays are delay differential equations (DDE). DDEs feature additional flexibility and properties, realize more complex dynamics and can complementary be used together with TCMs. We introduce several delay based PKPD models and investigate mathematical properties of general DDE based models, which serve as subunits in order to build larger PKPD models. Finally, we review current PKPD software with respect to the implementation of DDEs for PKPD analysis.
-
Lyapunov stability and its application to systems of ordinary differential equations
Kennedy, E. W.
1979-01-01
An outline and a brief introduction to some of the concepts and implications of Lyapunov stability theory are presented. Various aspects of the theory are illustrated by the inclusion of eight examples, including the Cartesian coordinate equations of the two-body problem, linear and nonlinear (Van der Pol's equation) oscillatory systems, and the linearized Kustaanheimo-Stiefel element equations for the unperturbed two-body problem.
-
Beginning partial differential equations
O'Neil, Peter V
2011-01-01
A rigorous, yet accessible, introduction to partial differential equations-updated in a valuable new edition Beginning Partial Differential Equations, Second Edition provides a comprehensive introduction to partial differential equations (PDEs) with a special focus on the significance of characteristics, solutions by Fourier series, integrals and transforms, properties and physical interpretations of solutions, and a transition to the modern function space approach to PDEs. With its breadth of coverage, this new edition continues to present a broad introduction to the field, while also addres
-
Application of the B-Determining Equations Method to One Problem of Free Turbulence
Directory of Open Access Journals (Sweden)
Oleg V. Kaptsov
2012-10-01
Full Text Available A three-dimensional model of the far turbulent wake behind a self-propelled body in a passively stratified medium is considered. The model is reduced to a system of ordinary differential equations by a similarity reduction and the B-determining equations method. The system of ordinary differential equations satisfying natural boundary conditions is solved numerically. The solutions obtained here are in close agreement with experimental data.
-
International Nuclear Information System (INIS)
Sidell, J.
1976-08-01
EXTRA is a program written for the Winfrith KDF9 enabling the user to solve first order initial value differential equations. In this report general numerical integration methods are discussed with emphasis on their application to the solution of stiff sets of equations. A method of particular applicability to stiff sets of equations is described. This method is incorporated in the program EXTRA and full instructions for its use are given. A comparison with other methods of computation is included. (author)
-
International Nuclear Information System (INIS)
Datta, Dhurjati Prasad
2003-01-01
A new class of finitely differentiable scale free solutions to the simplest class of ordinary differential equations is presented. Consequently, the real number set gets replaced by an extended physical set, each element of which is endowed with an equivalence class of infinitesimally separated neighbours in the form of random fluctuations. We show how a sense of time and evolution is intrinsically defined by the infinite continued fraction of the golden mean irrational number (Radical radicand 5 -1)/2, which plays a key role in this extended SL(2,R) formalism of calculus analogous to El Naschie's theory of E (∞) spacetime manifold. Time may thereby undergo random inversions generating well defined random scales, thus allowing a dynamical system to evolve self similarly over the set of multiple scales. The late time stochastic fluctuations of a dynamical system enjoys the generic 1/f spectrum. A universal form of the related probability density is also derived. We prove that the golden mean number is intrinsically random, letting all measurements in the physical universe fundamentally uncertain. The present analysis offers an explanation of the universal occurrence of the golden mean in diverse natural and biological processes as well as the mass spectrum of high energy particle physics
-
On differential operators generating iterative systems of linear ODEs of maximal symmetry algebra
Ndogmo, J. C.
2017-06-01
Although every iterative scalar linear ordinary differential equation is of maximal symmetry algebra, the situation is different and far more complex for systems of linear ordinary differential equations, and an iterative system of linear equations need not be of maximal symmetry algebra. We illustrate these facts by examples and derive families of vector differential operators whose iterations are all linear systems of equations of maximal symmetry algebra. Some consequences of these results are also discussed.
-
Hyperbolic partial differential equations
Witten, Matthew
1986-01-01
Hyperbolic Partial Differential Equations III is a refereed journal issue that explores the applications, theory, and/or applied methods related to hyperbolic partial differential equations, or problems arising out of hyperbolic partial differential equations, in any area of research. This journal issue is interested in all types of articles in terms of review, mini-monograph, standard study, or short communication. Some studies presented in this journal include discretization of ideal fluid dynamics in the Eulerian representation; a Riemann problem in gas dynamics with bifurcation; periodic M
-
Hemmerling, R.; Evers, J.B.; Smolenova, K.; Buck-Sorlin, G.H.; Kurth, W.
2013-01-01
In simulation models of plant development, physiological processes taking place in plants are typically described in terms of ODEs (Ordinary Differential Equations). On the one hand, those processes drive the development of the plant structure and on the other hand, the developed structure again
-
Raychaudhuri, D.
2007-01-01
The focus of this paper is on student interpretation and usage of the existence and uniqueness theorems for first-order ordinary differential equations. The inherent structure of the theorems is made explicit by the introduction of a framework of layers concepts-conditions-connectives-conclusions, and we discuss the manners in which students'…
-
On matrix fractional differential equations
Directory of Open Access Journals (Sweden)
Adem Kılıçman
2017-01-01
Full Text Available The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objective of this article is to discuss the Laplace transform method based on operational matrices of fractional derivatives for solving several kinds of linear fractional differential equations. Moreover, we present the operational matrices of fractional derivatives with Laplace transform in many applications of various engineering systems as control system. We present the analytical technique for solving fractional-order, multi-term fractional differential equation. In other words, we propose an efficient algorithm for solving fractional matrix equation.
-
Differential equations extended to superspace
International Nuclear Information System (INIS)
Torres, J.; Rosu, H.C.
2003-01-01
We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)
-
Differential equations extended to superspace
Energy Technology Data Exchange (ETDEWEB)
Torres, J. [Instituto de Fisica, Universidad de Guanajuato, A.P. E-143, Leon, Guanajuato (Mexico); Rosu, H.C. [Instituto Potosino de Investigacion Cientifica y Tecnologica, A.P. 3-74, Tangamanga, San Luis Potosi (Mexico)
2003-07-01
We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)
-
Beginning partial differential equations
O'Neil, Peter V
2014-01-01
A broad introduction to PDEs with an emphasis on specialized topics and applications occurring in a variety of fields Featuring a thoroughly revised presentation of topics, Beginning Partial Differential Equations, Third Edition provides a challenging, yet accessible,combination of techniques, applications, and introductory theory on the subjectof partial differential equations. The new edition offers nonstandard coverageon material including Burger's equation, the telegraph equation, damped wavemotion, and the use of characteristics to solve nonhomogeneous problems. The Third Edition is or
-
Solutions of the finite type of Sine-Gordon equation
International Nuclear Information System (INIS)
Zhao Guosong
1998-01-01
We use the technique of differential geometry to prove that the solutions of finite type of the sine-Gordon equation φ xx - φ yy = sin φ cosφ can be obtained from a system of ordinary differential equations
-
Differential equations problem solver
Arterburn, David R
2012-01-01
REA's Problem Solvers is a series of useful, practical, and informative study guides. Each title in the series is complete step-by-step solution guide. The Differential Equations Problem Solver enables students to solve difficult problems by showing them step-by-step solutions to Differential Equations problems. The Problem Solvers cover material ranging from the elementary to the advanced and make excellent review books and textbook companions. They're perfect for undergraduate and graduate studies.The Differential Equations Problem Solver is the perfect resource for any class, any exam, and
-
Direct 'delay' reductions of the Toda equation
International Nuclear Information System (INIS)
Joshi, Nalini
2009-01-01
A new direct method of obtaining reductions of the Toda equation is described. We find a canonical and complete class of all possible reductions under certain assumptions. The resulting equations are ordinary differential-difference equations, sometimes referred to as delay-differential equations. The representative equation of this class is hypothesized to be a new version of one of the classical Painleve equations. The Lax pair associated with this equation is obtained, also by reduction. (fast track communication)
-
International Nuclear Information System (INIS)
Salat, A.
1981-12-01
The existence of quasi-periodic eigensolutions of a linear second order ordinary differential equation with quasi-periodic coefficient f(ω 1 t,ω 2 t) is investigated numerically and graphically. For sufficiently incommensurate frequencies ω 1 , ω 2 a doubly indexed infinite sequence of eigenvalues and eigenmodes is obtained. The equation considered is a model for the magneto-hydrodynamic 'continuum' in general toroidal geometry. The result suggests that continuum modes exist at least on sufficiently irrational magnetic surfaces. (orig.)
-
DEFF Research Database (Denmark)
Tornøe, Christoffer Wenzel; Agersø, Henrik; Madsen, Henrik
2004-01-01
The standard software for non-linear mixed-effect analysis of pharmacokinetic/phar-macodynamic (PK/PD) data is NONMEM while the non-linear mixed-effects package NLME is an alternative as tong as the models are fairly simple. We present the nlmeODE package which combines the ordinary differential...... equation (ODE) solver package odesolve and the non-Linear mixed effects package NLME thereby enabling the analysis of complicated systems of ODEs by non-linear mixed-effects modelling. The pharmacokinetics of the anti-asthmatic drug theophylline is used to illustrate the applicability of the nlme...
-
Applications of Lie-group methods to the equations of magnetohydrodynamics
International Nuclear Information System (INIS)
Mandrekas, J.
1987-01-01
The invariance properties of various sets of magnetohydrodynamic (MHD) equations are studied using techniques from the theory of differential forms. Equations considered include the ideal MHD equations in different geometries and with different magnetic field configurations, the MHD equations in the presence of gravitational forces due to self-attraction or external fields, and the MHD equations including finite thermal conductivity and magnetic viscosity. The knowledge of the group structure of these equations is then used to introduce similarity variables to these equations. For each choice of similarity variables, the original set of partial differential equations is transformed into a set of ordinary differential equations and the most general form of the initial conditions is determined. Three cases are studied in detail and the corresponding sets of ordinary differential equations are solved numerically: the problem of a blast wave in an inhomogeneous atmosphere, the problem of a piston moving according to a power law in time, and the problem of a piston moving according to an exponential law in time
-
Linear determining equations for differential constraints
International Nuclear Information System (INIS)
Kaptsov, O V
1998-01-01
A construction of differential constraints compatible with partial differential equations is considered. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the classical determining equations used in the search for admissible Lie operators. As applications of this approach equations of an ideal incompressible fluid and non-linear heat equations are discussed
-
Functional differential equations—a reciprocity principle
Directory of Open Access Journals (Sweden)
Lloyd K. Williams
1986-01-01
Full Text Available The functional differential equations proposed for solution here are mainly ordinary differential equations with fairly general argument deviations. Included among them are equations with involutions and some with reflections of the argument. Solutions will be obtained by quadratures in terms of implicitly defined functions. They have a wide range of applicability from the stability theory of differential-difference equations to electrodynamics and biological models.
-
Directory of Open Access Journals (Sweden)
A. H. Bhrawy
2014-01-01
Full Text Available One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initial-boundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initial-boundary coupled hyperbolic PDEs with variable coefficients as well as initial-nonlocal conditions. Using triangular, soliton, and exponential-triangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate.
-
On double reductions from symmetries and conservation laws for a damped Boussinesq equation
International Nuclear Information System (INIS)
Gandarias, M.L.; Rosa, M.
2016-01-01
In this work, we study a Boussinesq equation with a strong damping term from the point of view of the Lie theory. We derive the classical Lie symmetries admitted by the equation as well as the reduced ordinary differential equations. Some nontrivial conservation laws are derived by using the multipliers method. Taking into account the relationship between symmetries and conservation laws and applying the double reduction method, we obtain a direct reduction of order of the ordinary differential equations and in particular a kink solution.
-
Numerical solutions of the Vlasov equation
International Nuclear Information System (INIS)
Satofuka, Nobuyuki; Morinishi, Koji; Nishida, Hidetoshi
1985-01-01
A numerical procedure is derived for the solutions of the one- and two-dimensional Vlasov-Poisson system equations. This numerical procedure consists of the phase space discretization and the integration of the resulting set of ordinary differential equations. In the phase space discretization, derivatives with respect to the phase space variable are approximated by a weighted sum of the values of the distribution function at properly chosen neighboring points. Then, the resulting set of ordinary differential equations is solved by using an appropriate time integration scheme. The results for linear Landau damping, nonlinear Landau damping and counter-streaming plasmas are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient. (author)
-
Algebraic entropy for differential-delay equations
Viallet, Claude M.
2014-01-01
We extend the definition of algebraic entropy to a class of differential-delay equations. The vanishing of the entropy, as a structural property of an equation, signals its integrability. We suggest a simple way to produce differential-delay equations with vanishing entropy from known integrable differential-difference equations.
-
Exact solutions of generalized Zakharov and Ginzburg-Landau equations
International Nuclear Information System (INIS)
Zhang Jinliang; Wang Mingliang; Gao Kequan
2007-01-01
By using the homogeneous balance principle, the exact solutions of the generalized Zakharov equations and generalized Ginzburg-Landau equation are obtained with the aid of a set of subsidiary higher-order ordinary differential equations (sub-equations for short)
-
Schiesser, William E
2014-01-01
Features a solid foundation of mathematical and computational tools to formulate and solve real-world PDE problems across various fields With a step-by-step approach to solving partial differential equations (PDEs), Differential Equation Analysis in Biomedical Science and Engineering: Partial Differential Equation Applications with R successfully applies computational techniques for solving real-world PDE problems that are found in a variety of fields, including chemistry, physics, biology, and physiology. The book provides readers with the necessary knowledge to reproduce and extend the com
-
Converting differential-equation models of biological systems to membrane computing.
Muniyandi, Ravie Chandren; Zin, Abdullah Mohd; Sanders, J W
2013-12-01
This paper presents a method to convert the deterministic, continuous representation of a biological system by ordinary differential equations into a non-deterministic, discrete membrane computation. The dynamics of the membrane computation is governed by rewrite rules operating at certain rates. That has the advantage of applying accurately to small systems, and to expressing rates of change that are determined locally, by region, but not necessary globally. Such spatial information augments the standard differentiable approach to provide a more realistic model. A biological case study of the ligand-receptor network of protein TGF-β is used to validate the effectiveness of the conversion method. It demonstrates the sense in which the behaviours and properties of the system are better preserved in the membrane computing model, suggesting that the proposed conversion method may prove useful for biological systems in particular. Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.
-
On matrix fractional differential equations
Adem Kılıçman; Wasan Ajeel Ahmood
2017-01-01
The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objec...
-
Local p-Adic Differential Equations
Put, Marius van der; Taelman, Lenny
2006-01-01
This paper studies divergence in solutions of p-adic linear local differential equations. Such divergence is related to the notion of p-adic Liouville numbers. Also, the influence of the divergence on the differential Galois groups of such differential equations is explored. A complete result is
-
Tricomi, FG
2013-01-01
Based on his extensive experience as an educator, F. G. Tricomi wrote this practical and concise teaching text to offer a clear idea of the problems and methods of the theory of differential equations. The treatment is geared toward advanced undergraduates and graduate students and addresses only questions that can be resolved with rigor and simplicity.Starting with a consideration of the existence and uniqueness theorem, the text advances to the behavior of the characteristics of a first-order equation, boundary problems for second-order linear equations, asymptotic methods, and diff
-
A Note on Unsteady Temperature Equation For Gravity Flow of A ...
African Journals Online (AJOL)
We present an analytical study of unsteady temperature energy equation for gravity of a fluid with non – Newtonian behaviour through a porous medium. For the case of radial axisymmetric flow, the governing partial differential equation is transformed into an ordinary differential equation through similarity variables.
-
Stability analysis for a delay differential equations model of a hydraulic turbine speed governor
Halanay, Andrei; Safta, Carmen A.; Dragoi, Constantin; Piraianu, Vlad F.
2017-01-01
The paper aims to study the dynamic behavior of a speed governor for a hydraulic turbine using a mathematical model. The nonlinear mathematical model proposed consists in a system of delay differential equations (DDE) to be compared with already established mathematical models of ordinary differential equations (ODE). A new kind of nonlinearity is introduced as a time delay. The delays can characterize different running conditions of the speed governor. For example, it is considered that spool displacement of hydraulic amplifier might be blocked due to oil impurities in the oil supply system and so the hydraulic amplifier has a time delay in comparison to the time control. Numerical simulations are presented in a comparative manner. A stability analysis of the hydraulic control system is performed, too. Conclusions of the dynamic behavior using the DDE model of a hydraulic turbine speed governor are useful in modeling and controlling hydropower plants.
-
International Nuclear Information System (INIS)
Tian Ye; Chen Jing; Zhang Zhifei
2012-01-01
In this paper, the separation transformation approach is extended to the (N + 1)-dimensional dispersive double sine-Gordon equation arising in many physical systems such as the spin dynamics in the B phase of 3 He superfluid. This equation is first reduced to a set of partial differential equations and a nonlinear ordinary differential equation. Then the general solutions of the set of partial differential equations are obtained and the nonlinear ordinary differential equation is solved by F-expansion method. Finally, many new exact solutions of the (N + 1)-dimensional dispersive double sine-Gordon equation are constructed explicitly via the separation transformation. For the case of N > 2, there is an arbitrary function in the exact solutions, which may reveal more novel nonlinear structures in the high-dimensional dispersive double sine-Gordon equation.
-
Sloss, J. M.; Kranzler, S. K.
1972-01-01
The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.
-
Exact solutions of the Drinfel'd–Sokolov–Wilson equation using ...
Indian Academy of Sciences (India)
1Department of Engineering Mathematics and Physics, Higher Institute of Engineering, .... ordinary differential equation (NLODE) of the form. P .... 3.1 The Bäcklund transformation of Riccati equation method applied to DSW equation.
-
A class of exact solutions to the Einstein field equations
International Nuclear Information System (INIS)
Goyal, Nisha; Gupta, R K
2012-01-01
The Einstein-Rosen metric is considered and a class of new exact solutions of the field equations for stationary axisymmetric Einstein-Maxwell fields is obtained. The Lie classical approach is applied to obtain exact solutions. By using the Lie classical method, we are able to derive symmetries that are used for reducing the coupled system of partial differential equations into ordinary differential equations. From reduced differential equations we have derived some new exact solutions of Einstein-Maxwell equations. (paper)
-
Double compactons in the Olver–Rosenau equation
Indian Academy of Sciences (India)
2013-03-01
Mar 1, 2013 ... are treated by the dynamical systems theory and a phase-space analysis ... The distinguishing feature of the systems ... This equation has attracted a lot of attention ... integration over ξ leads to an ordinary differential equation.
-
Leander, Jacob; Lundh, Torbjörn; Jirstrand, Mats
2014-05-01
In this paper we consider the problem of estimating parameters in ordinary differential equations given discrete time experimental data. The impact of going from an ordinary to a stochastic differential equation setting is investigated as a tool to overcome the problem of local minima in the objective function. Using two different models, it is demonstrated that by allowing noise in the underlying model itself, the objective functions to be minimized in the parameter estimation procedures are regularized in the sense that the number of local minima is reduced and better convergence is achieved. The advantage of using stochastic differential equations is that the actual states in the model are predicted from data and this will allow the prediction to stay close to data even when the parameters in the model is incorrect. The extended Kalman filter is used as a state estimator and sensitivity equations are provided to give an accurate calculation of the gradient of the objective function. The method is illustrated using in silico data from the FitzHugh-Nagumo model for excitable media and the Lotka-Volterra predator-prey system. The proposed method performs well on the models considered, and is able to regularize the objective function in both models. This leads to parameter estimation problems with fewer local minima which can be solved by efficient gradient-based methods. Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.
-
Abstract methods in partial differential equations
Carroll, Robert W
2012-01-01
Detailed, self-contained treatment examines modern abstract methods in partial differential equations, especially abstract evolution equations. Suitable for graduate students with some previous exposure to classical partial differential equations. 1969 edition.
-
ANALYTICAL SOLUTION OF THE K-TH ORDER AUTONOMOUS ORDINARY DIFFERENTIAL EQUATION
Directory of Open Access Journals (Sweden)
Ronald Orozco López
2017-04-01
Full Text Available The main objective of this paper is to find the analytical solution of the autonomous equation y(k = f (y and prove its convergence using autonomous polynomials of order k, define here in addition of the formula of Faá di Bruno for composition of functions and Bell polynomials. Autonomous polynomials of order k are defined in terms of the boundary values of the equation. Also special values of autonomous polynomials of order 1 are given.
-
Five-dimensional Monopole Equation with Hedge-Hog Ansatz and Abel's Differential Equation
Kihara, Hironobu
2008-01-01
We review the generalized monopole in the five-dimensional Euclidean space. A numerical solution with the Hedge-Hog ansatz is studied. The Bogomol'nyi equation becomes a second order autonomous non-linear differential equation. The equation can be translated into the Abel's differential equation of the second kind and is an algebraic differential equation.
-
Mathieu's Equation and its Generalizations: Overview of Stability Charts and their Features
DEFF Research Database (Denmark)
Kovacic, Ivana; Rand, Richard H.; Sah, Si Mohamed
2018-01-01
This work is concerned with Mathieu's equation - a classical differential equation, which has the form of a linear second-order ordinary differential equation with Cosine-type periodic forcing of the stiffness coefficient, and its different generalizations/extensions. These extensions include...... and features, and how it differs from that of the classical Mathieu's equation....
-
Computational partial differential equations using Matlab
Li, Jichun
2008-01-01
Brief Overview of Partial Differential Equations The parabolic equations The wave equations The elliptic equations Differential equations in broader areasA quick review of numerical methods for PDEsFinite Difference Methods for Parabolic Equations Introduction Theoretical issues: stability, consistence, and convergence 1-D parabolic equations2-D and 3-D parabolic equationsNumerical examples with MATLAB codesFinite Difference Methods for Hyperbolic Equations IntroductionSome basic difference schemes Dissipation and dispersion errors Extensions to conservation lawsThe second-order hyperbolic PDE
-
New exact solutions of the Einstein—Maxwell equations for magnetostatic fields
International Nuclear Information System (INIS)
Goyal, Nisha; Gupta, R.K.
2012-01-01
The symmetry reduction method based on the Fréchet derivative of differential operators is applied to investigate symmetries of the Einstein—Maxwell field equations for magnetostatic fields, which is a coupled system of nonlinear partial differential equations of the second order. The technique yields invariant transformations that reduce the given system of partial differential equations to a system of nonlinear ordinary differential equations. Some of the reduced systems are further studied to obtain the exact solutions
-
Fuchsia. A tool for reducing differential equations for Feynman master integral to epsilon form
Energy Technology Data Exchange (ETDEWEB)
Gituliar, Oleksandr [Hamburg Univ. (Germany). 2. Inst. fuer Theoretische Physik; Magerya, Vitaly
2017-01-15
We present Fuchsia - an implementation of the Lee algorithm, which for a given system of ordinary differential equations with rational coefficients ∂{sub x}f(x,ε)=A(x,ε)f(x,ε) finds a basis transformation T(x,ε), i.e., f(x,ε)=T(x,ε)g(x,ε), such that the system turns into the epsilon form: ∂{sub x}g(x,ε)=εS(x)g(x,ε), where S(x) is a Fuchsian matrix. A system of this form can be trivially solved in terms of polylogarithms as a Laurent series in the dimensional regulator ε. That makes the construction of the transformation T(x,ε) crucial for obtaining solutions of the initial system. In principle, Fuchsia can deal with any regular systems, however its primary task is to reduce differential equations for Feynman master integrals. It ensures that solutions contain only regular singularities due to the properties of Feynman integrals.
-
Fuchsia. A tool for reducing differential equations for Feynman master integral to epsilon form
International Nuclear Information System (INIS)
Gituliar, Oleksandr; Magerya, Vitaly
2017-01-01
We present Fuchsia - an implementation of the Lee algorithm, which for a given system of ordinary differential equations with rational coefficients ∂ x f(x,ε)=A(x,ε)f(x,ε) finds a basis transformation T(x,ε), i.e., f(x,ε)=T(x,ε)g(x,ε), such that the system turns into the epsilon form: ∂ x g(x,ε)=εS(x)g(x,ε), where S(x) is a Fuchsian matrix. A system of this form can be trivially solved in terms of polylogarithms as a Laurent series in the dimensional regulator ε. That makes the construction of the transformation T(x,ε) crucial for obtaining solutions of the initial system. In principle, Fuchsia can deal with any regular systems, however its primary task is to reduce differential equations for Feynman master integrals. It ensures that solutions contain only regular singularities due to the properties of Feynman integrals.
-
New traveling wave solutions to AKNS and SKdV equations
International Nuclear Information System (INIS)
Ozer, Teoman
2009-01-01
We analyze the traveling wave solutions of Ablowitz-Kaup-Newell-Segur (AKNS) and Schwarz-Korteweg-de Vries (SKdV) equations. As the solution method for differential equations we consider the improved tanh approach. This approach provides to transform the partial differential equation into the ordinary differential equation and then obtain the new families of exact solutions based on the solutions of the Riccati equation. The different values of the coefficients of the Riccati equation allow us to obtain new type of traveling wave solutions to AKNS and SKdV equations.
-
Exact traveling wave solutions of the Boussinesq equation
International Nuclear Information System (INIS)
Ding Shuangshuang; Zhao Xiqiang
2006-01-01
The repeated homogeneous balance method is used to construct exact traveling wave solutions of the Boussinesq equation, in which the homogeneous balance method is applied to solve the Riccati equation and the reduced nonlinear ordinary differential equation, respectively. Many new exact traveling wave solutions of the Boussinesq equation are successfully obtained
-
Model identification using stochastic differential equation grey-box models in diabetes.
Duun-Henriksen, Anne Katrine; Schmidt, Signe; Røge, Rikke Meldgaard; Møller, Jonas Bech; Nørgaard, Kirsten; Jørgensen, John Bagterp; Madsen, Henrik
2013-03-01
The acceptance of virtual preclinical testing of control algorithms is growing and thus also the need for robust and reliable models. Models based on ordinary differential equations (ODEs) can rarely be validated with standard statistical tools. Stochastic differential equations (SDEs) offer the possibility of building models that can be validated statistically and that are capable of predicting not only a realistic trajectory, but also the uncertainty of the prediction. In an SDE, the prediction error is split into two noise terms. This separation ensures that the errors are uncorrelated and provides the possibility to pinpoint model deficiencies. An identifiable model of the glucoregulatory system in a type 1 diabetes mellitus (T1DM) patient is used as the basis for development of a stochastic-differential-equation-based grey-box model (SDE-GB). The parameters are estimated on clinical data from four T1DM patients. The optimal SDE-GB is determined from likelihood-ratio tests. Finally, parameter tracking is used to track the variation in the "time to peak of meal response" parameter. We found that the transformation of the ODE model into an SDE-GB resulted in a significant improvement in the prediction and uncorrelated errors. Tracking of the "peak time of meal absorption" parameter showed that the absorption rate varied according to meal type. This study shows the potential of using SDE-GBs in diabetes modeling. Improved model predictions were obtained due to the separation of the prediction error. SDE-GBs offer a solid framework for using statistical tools for model validation and model development. © 2013 Diabetes Technology Society.
-
Skew differential fields, differential and difference equations
van der Put, M
2004-01-01
The central question is: Let a differential or difference equation over a field K be isomorphic to all its Galois twists w.r.t. the group Gal(K/k). Does the equation descend to k? For a number of categories of equations an answer is given.
-
Weak self-adjoint differential equations
International Nuclear Information System (INIS)
Gandarias, M L
2011-01-01
The concepts of self-adjoint and quasi self-adjoint equations were introduced by Ibragimov (2006 J. Math. Anal. Appl. 318 742-57; 2007 Arch. ALGA 4 55-60). In Ibragimov (2007 J. Math. Anal. Appl. 333 311-28), a general theorem on conservation laws was proved. In this paper, we generalize the concept of self-adjoint and quasi self-adjoint equations by introducing the definition of weak self-adjoint equations. We find a class of weak self-adjoint quasi-linear parabolic equations. The property of a differential equation to be weak self-adjoint is important for constructing conservation laws associated with symmetries of the differential equation. (fast track communication)
-
Yang, Ting; Dong, Jianji; Lu, Liangjun; Zhou, Linjie; Zheng, Aoling; Zhang, Xinliang; Chen, Jianping
2014-07-04
Photonic integrated circuits for photonic computing open up the possibility for the realization of ultrahigh-speed and ultra wide-band signal processing with compact size and low power consumption. Differential equations model and govern fundamental physical phenomena and engineering systems in virtually any field of science and engineering, such as temperature diffusion processes, physical problems of motion subject to acceleration inputs and frictional forces, and the response of different resistor-capacitor circuits, etc. In this study, we experimentally demonstrate a feasible integrated scheme to solve first-order linear ordinary differential equation with constant-coefficient tunable based on a single silicon microring resonator. Besides, we analyze the impact of the chirp and pulse-width of input signals on the computing deviation. This device can be compatible with the electronic technology (typically complementary metal-oxide semiconductor technology), which may motivate the development of integrated photonic circuits for optical computing.
-
International Nuclear Information System (INIS)
Winkler, E.
1991-01-01
Mathematical models in tracer kinetics are usually based on ordinary differential equations which correspond to a system of kinetically homogeneous compartments (standard compartments). A generalization is possible by the admission of inhomogeneities in the behaviour of the elements belonging to a compartment. The important special case of the age-dependence of elimination rates is treated in its deterministic version. It leads to partial different equations (i.e., systems with distributed coefficients) with the 'age' or the 'residence time' of an element of the compartment as a variable additional to 'time'. The basic equations for one generalized compartment and for systems of such compartments are given together with their general solutions. (orig.) [de
-
Partial differential equations of mathematical physics
Sobolev, S L
1964-01-01
Partial Differential Equations of Mathematical Physics emphasizes the study of second-order partial differential equations of mathematical physics, which is deemed as the foundation of investigations into waves, heat conduction, hydrodynamics, and other physical problems. The book discusses in detail a wide spectrum of topics related to partial differential equations, such as the theories of sets and of Lebesgue integration, integral equations, Green's function, and the proof of the Fourier method. Theoretical physicists, experimental physicists, mathematicians engaged in pure and applied math
-
Soliton surfaces via a zero-curvature representation of differential equations
International Nuclear Information System (INIS)
Grundland, A M; Post, S
2012-01-01
The main aim of this paper is to introduce a new version of the Fokas–Gel’fand formula for immersion of soliton surfaces in Lie algebras. The paper contains a detailed exposition of the technique for obtaining exact forms of 2D surfaces associated with any solution of a given nonlinear ordinary differential equation which can be written in the zero-curvature form. That is, for any generalized symmetry of the zero-curvature condition of the associated integrable model, it is possible to construct soliton surfaces whose Gauss–Mainardi–Codazzi equations are equivalent to infinitesimal deformations of the zero-curvature representation of the considered model. Conversely, it is shown (proposition 1) that for a given immersion function of a 2D soliton surface in a Lie algebra, it is possible to derive the associated generalized vector field in the evolutionary form which characterizes all symmetries of the zero-curvature condition. The theoretical considerations are illustrated via surfaces associated with the Painlevé equations P1, P2 and P3, including transcendental functions, the special cases of the rational and Airy solutions of P2 and the classical solutions of P3. (paper)
-
Similarity solutions of the Fokker–Planck equation with time-dependent coefficients
International Nuclear Information System (INIS)
Lin, W.-T.; Ho, C.-L.
2012-01-01
In this work, we consider the solvability of the Fokker–Planck equation with both time-dependent drift and diffusion coefficients by means of the similarity method. By the introduction of the similarity variable, the Fokker–Planck equation is reduced to an ordinary differential equation. Adopting the natural requirement that the probability current density vanishes at the boundary, the resulting ordinary differential equation turns out to be integrable, and the probability density function can be given in closed form. New examples of exactly solvable Fokker–Planck equations are presented, and their properties analyzed. - Highlights: ► Scaling form of the Fokker–Planck equation with time-dependent drift and diffusion coefficients is derived. ► Exact similarity solution of the Fokker–Planck equation is given in closed forms. ► New examples of Fokker–Planck equations exactly solvable by similarity methods are discussed.
-
On Degenerate Partial Differential Equations
Chen, Gui-Qiang G.
2010-01-01
Some of recent developments, including recent results, ideas, techniques, and approaches, in the study of degenerate partial differential equations are surveyed and analyzed. Several examples of nonlinear degenerate, even mixed, partial differential equations, are presented, which arise naturally in some longstanding, fundamental problems in fluid mechanics and differential geometry. The solution to these fundamental problems greatly requires a deep understanding of nonlinear degenerate parti...
-
Differential equations and finite groups
Put, Marius van der; Ulmer, Felix
2000-01-01
The classical solution of the Riemann-Hilbert problem attaches to a given representation of the fundamental group a regular singular linear differential equation. We present a method to compute this differential equation in the case of a representation with finite image. The approach uses Galois
-
Violon, D
2012-12-01
To develop a multicompartment model of only essential human body components that predicts the contrast medium concentration vs time curve in a chosen compartment after an intravenous injection. Also to show that the model can be used to time adequately contrast-enhanced CT series. A system of linked time delay instead of ordinary differential equations described the model and was solved with a Matlab program (Matlab v. 6.5; The Mathworks, Inc., Natick, MA). All the injection and physiological parameters were modified to cope with normal or pathological situations. In vivo time-concentration curves from the literature were recalculated to validate the model. The recalculated contrast medium time-concentration curves and parameters are given. The results of the statistical analysis of the study findings are expressed as the median prediction error and the median absolute prediction error values for both the time delay and ordinary differential equation systems; these are situated well below the generally accepted maximum 20% limit. The presented program correctly predicts the time-concentration curve of an intravenous contrast medium injection and, consequently, allows an individually tailored approach of CT examinations with optimised use of the injected contrast medium volume, as long as time delay instead of ordinary differential equations are used. The presented program offers good preliminary knowledge of the time-contrast medium concentration curve after any intravenous injection, allowing adequate timing of a CT examination, required by the short scan time of present-day scanners. The injected volume of contrast medium can be tailored to the individual patient with no more contrast medium than is strictly needed.
-
Discontinuous gradient differential equations and trajectories in the calculus of variations
International Nuclear Information System (INIS)
Bogaevskii, I A
2006-01-01
The concept of gradient of smooth functions is generalized for their sums with concave functions. An existence, uniqueness, and continuous dependence theorem for increasing time is formulated and proved for solutions of an ordinary differential equation the right-hand side of which is the gradient of the sum of a concave and a smooth function. With the use of this result a physically natural motion of particles, well defined even at discontinuities of the velocity field, is constructed in the variational problem of the minimal mechanical action in a space of arbitrary dimension. For such a motion of particles in the plane all typical cases of the birth and the interaction of point clusters of positive mass are described.
-
Discontinuous gradient differential equations and trajectories in the calculus of variations
Bogaevskii, I. A.
2006-12-01
The concept of gradient of smooth functions is generalized for their sums with concave functions. An existence, uniqueness, and continuous dependence theorem for increasing time is formulated and proved for solutions of an ordinary differential equation the right-hand side of which is the gradient of the sum of a concave and a smooth function. With the use of this result a physically natural motion of particles, well defined even at discontinuities of the velocity field, is constructed in the variational problem of the minimal mechanical action in a space of arbitrary dimension. For such a motion of particles in the plane all typical cases of the birth and the interaction of point clusters of positive mass are described.
-
Statistical Methods for Stochastic Differential Equations
Kessler, Mathieu; Sorensen, Michael
2012-01-01
The seventh volume in the SemStat series, Statistical Methods for Stochastic Differential Equations presents current research trends and recent developments in statistical methods for stochastic differential equations. Written to be accessible to both new students and seasoned researchers, each self-contained chapter starts with introductions to the topic at hand and builds gradually towards discussing recent research. The book covers Wiener-driven equations as well as stochastic differential equations with jumps, including continuous-time ARMA processes and COGARCH processes. It presents a sp
-
ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY
Enrique Gonzalo Reyes Garcia
2004-01-01
ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY Equations in partial derivatives appeared in the 18th century as essential tools for the analytic study of physical models and, later, they proved to be fundamental for the progress of mathematics. For example, fundamental results of modern differential geometry are based on deep theorems on differential equations. Reciprocally, it is possible to study differential equations through geometrical means just like it was done by o...
-
Fermat type differential and difference equations
Directory of Open Access Journals (Sweden)
Kai Liu
2015-06-01
Full Text Available This article we explore the relationship between the number of differential and difference operators with the existence of meromorphic solutions of Fermat type differential and difference equations. Some Fermat differential and difference equations of certain types are also considered.
-
PARALLEL SOLUTION METHODS OF PARTIAL DIFFERENTIAL EQUATIONS
Directory of Open Access Journals (Sweden)
Korhan KARABULUT
1998-03-01
Full Text Available Partial differential equations arise in almost all fields of science and engineering. Computer time spent in solving partial differential equations is much more than that of in any other problem class. For this reason, partial differential equations are suitable to be solved on parallel computers that offer great computation power. In this study, parallel solution to partial differential equations with Jacobi, Gauss-Siedel, SOR (Succesive OverRelaxation and SSOR (Symmetric SOR algorithms is studied.
-
Introduction to complex theory of differential equations
Savin, Anton
2017-01-01
This book discusses the complex theory of differential equations or more precisely, the theory of differential equations on complex-analytic manifolds. Although the theory of differential equations on real manifolds is well known – it is described in thousands of papers and its usefulness requires no comments or explanations – to date specialists on differential equations have not focused on the complex theory of partial differential equations. However, as well as being remarkably beautiful, this theory can be used to solve a number of problems in real theory, for instance, the Poincaré balayage problem and the mother body problem in geophysics. The monograph does not require readers to be familiar with advanced notions in complex analysis, differential equations, or topology. With its numerous examples and exercises, it appeals to advanced undergraduate and graduate students, and also to researchers wanting to familiarize themselves with the subject.
-
Optimal control problem for the extended Fisher–Kolmogorov equation
Indian Academy of Sciences (India)
by methods of optimal control, such as chemical engineering and vehicle ... ern optimal control theories and applied models are not only represented by .... Obviously, equation (2.5) is an ordinary differential equation and according to ODE.
-
International Nuclear Information System (INIS)
Khuri, S A; Sayfy, Ali
2013-01-01
This study is intended to provide a modified variational algorithm for the numerical solution of the third-order ordinary differential equation y″′ = y -2 which arises in the modeling of thin viscous films with surface tension. The resulting solution is used to solve a problem relevant to Tanner's Law for the speed of a moving three-phase contact line. Numerical results, computational comparisons, suitable error measures and illustrations are provided to testify and demonstrate the convergence and efficiency of the method.
-
Quasi-exact solutions of nonlinear differential equations
Kudryashov, Nikolay A.; Kochanov, Mark B.
2014-01-01
The concept of quasi-exact solutions of nonlinear differential equations is introduced. Quasi-exact solution expands the idea of exact solution for additional values of parameters of differential equation. These solutions are approximate solutions of nonlinear differential equations but they are close to exact solutions. Quasi-exact solutions of the the Kuramoto--Sivashinsky, the Korteweg--de Vries--Burgers and the Kawahara equations are founded.
-
On an improved method for solving evolution equations of higher ...
African Journals Online (AJOL)
In this paper we introduce a new algebraic procedure to compute new classes of solutions of (1+1)-nonlinear partial differential equations (nPDEs) both of physical and technical relevance. The basic assumption is that the unknown solution(s) of the nPDE under consideration satisfy an ordinary differential equation (ODE) of ...
-
XMDS2: Fast, scalable simulation of coupled stochastic partial differential equations
Dennis, Graham R.; Hope, Joseph J.; Johnsson, Mattias T.
2013-01-01
XMDS2 is a cross-platform, GPL-licensed, open source package for numerically integrating initial value problems that range from a single ordinary differential equation up to systems of coupled stochastic partial differential equations. The equations are described in a high-level XML-based script, and the package generates low-level optionally parallelised C++ code for the efficient solution of those equations. It combines the advantages of high-level simulations, namely fast and low-error development, with the speed, portability and scalability of hand-written code. XMDS2 is a complete redesign of the XMDS package, and features support for a much wider problem space while also producing faster code. Program summaryProgram title: XMDS2 Catalogue identifier: AENK_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 2 No. of lines in distributed program, including test data, etc.: 872490 No. of bytes in distributed program, including test data, etc.: 45522370 Distribution format: tar.gz Programming language: Python and C++. Computer: Any computer with a Unix-like system, a C++ compiler and Python. Operating system: Any Unix-like system; developed under Mac OS X and GNU/Linux. RAM: Problem dependent (roughly 50 bytes per grid point) Classification: 4.3, 6.5. External routines: The external libraries required are problem-dependent. Uses FFTW3 Fourier transforms (used only for FFT-based spectral methods), dSFMT random number generation (used only for stochastic problems), MPI message-passing interface (used only for distributed problems), HDF5, GNU Scientific Library (used only for Bessel-based spectral methods) and a BLAS implementation (used only for non-FFT-based spectral methods). Nature of problem: General coupled initial-value stochastic partial differential equations. Solution method: Spectral method
-
Solving delay differential equations in S-ADAPT by method of steps.
Bauer, Robert J; Mo, Gary; Krzyzanski, Wojciech
2013-09-01
S-ADAPT is a version of the ADAPT program that contains additional simulation and optimization abilities such as parametric population analysis. S-ADAPT utilizes LSODA to solve ordinary differential equations (ODEs), an algorithm designed for large dimension non-stiff and stiff problems. However, S-ADAPT does not have a solver for delay differential equations (DDEs). Our objective was to implement in S-ADAPT a DDE solver using the methods of steps. The method of steps allows one to solve virtually any DDE system by transforming it to an ODE system. The solver was validated for scalar linear DDEs with one delay and bolus and infusion inputs for which explicit analytic solutions were derived. Solutions of nonlinear DDE problems coded in S-ADAPT were validated by comparing them with ones obtained by the MATLAB DDE solver dde23. The estimation of parameters was tested on the MATLB simulated population pharmacodynamics data. The comparison of S-ADAPT generated solutions for DDE problems with the explicit solutions as well as MATLAB produced solutions which agreed to at least 7 significant digits. The population parameter estimates from using importance sampling expectation-maximization in S-ADAPT agreed with ones used to generate the data. Published by Elsevier Ireland Ltd.
-
New exact solutions to the generalized KdV equation with ...
Indian Academy of Sciences (India)
is reduced to an ordinary differential equation with constant coefficients ... Application to generalized KdV equation with generalized evolution ..... [12] P F Byrd and M D Friedman, Handbook of elliptic integrals for engineers and physicists.
-
The Landau-Lifshitz equation describes the Ising spin correlation function in the free-fermion model
Rutkevich, S B
1998-01-01
We consider time and space dependence of the Ising spin correlation function in a continuous one-dimensional free-fermion model. By the Ising spin we imply the 'sign' variable, which takes alternating +-1 values in adjacent domains bounded by domain walls (fermionic world paths). The two-point correlation function is expressed in terms of the solution of the Cauchy problem for a nonlinear partial differential equation, which is proved to be equivalent to the exactly solvable Landau-Lifshitz equation. A new zero-curvature representation for this equation is presented. In turn, the initial condition for the Cauchy problem is given by the solution of a nonlinear ordinary differential equation, which has also been derived. In the Ising limit the above-mentioned partial and ordinary differential equations reduce to the sine-Gordon and Painleve III equations, respectively. (author)
-
Exact solutions to two higher order nonlinear Schroedinger equations
International Nuclear Information System (INIS)
Xu Liping; Zhang Jinliang
2007-01-01
Using the homogeneous balance principle and F-expansion method, the exact solutions to two higher order nonlinear Schroedinger equations which describe the propagation of femtosecond pulses in nonlinear fibres are obtained with the aid of a set of subsidiary higher order ordinary differential equations (sub-equations for short)
-
Semiparametric mixed-effects analysis of PK/PD models using differential equations.
Wang, Yi; Eskridge, Kent M; Zhang, Shunpu
2008-08-01
Motivated by the use of semiparametric nonlinear mixed-effects modeling on longitudinal data, we develop a new semiparametric modeling approach to address potential structural model misspecification for population pharmacokinetic/pharmacodynamic (PK/PD) analysis. Specifically, we use a set of ordinary differential equations (ODEs) with form dx/dt = A(t)x + B(t) where B(t) is a nonparametric function that is estimated using penalized splines. The inclusion of a nonparametric function in the ODEs makes identification of structural model misspecification feasible by quantifying the model uncertainty and provides flexibility for accommodating possible structural model deficiencies. The resulting model will be implemented in a nonlinear mixed-effects modeling setup for population analysis. We illustrate the method with an application to cefamandole data and evaluate its performance through simulations.
-
Numerical Analysis of Partial Differential Equations
Lui, S H
2011-01-01
A balanced guide to the essential techniques for solving elliptic partial differential equations Numerical Analysis of Partial Differential Equations provides a comprehensive, self-contained treatment of the quantitative methods used to solve elliptic partial differential equations (PDEs), with a focus on the efficiency as well as the error of the presented methods. The author utilizes coverage of theoretical PDEs, along with the nu merical solution of linear systems and various examples and exercises, to supply readers with an introduction to the essential concepts in the numerical analysis
-
Stochastic differential equation model to Prendiville processes
International Nuclear Information System (INIS)
Granita; Bahar, Arifah
2015-01-01
The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution
-
Stochastic differential equation model to Prendiville processes
Energy Technology Data Exchange (ETDEWEB)
Granita, E-mail: granitafc@gmail.com [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); Bahar, Arifah [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); UTM Center for Industrial & Applied Mathematics (UTM-CIAM) (Malaysia)
2015-10-22
The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution.
-
Introduction to partial differential equations
Borthwick, David
2016-01-01
This modern take on partial differential equations does not require knowledge beyond vector calculus and linear algebra. The author focuses on the most important classical partial differential equations, including conservation equations and their characteristics, the wave equation, the heat equation, function spaces, and Fourier series, drawing on tools from analysis only as they arise.Within each section the author creates a narrative that answers the five questions: (1) What is the scientific problem we are trying to understand? (2) How do we model that with PDE? (3) What techniques can we use to analyze the PDE? (4) How do those techniques apply to this equation? (5) What information or insight did we obtain by developing and analyzing the PDE? The text stresses the interplay between modeling and mathematical analysis, providing a thorough source of problems and an inspiration for the development of methods.
-
New results for exponential synchronization of linearly coupled ordinary differential systems
International Nuclear Information System (INIS)
Tong Ping; Chen Shi-Hua
2017-01-01
This paper investigates the exponential synchronization of linearly coupled ordinary differential systems. The intrinsic nonlinear dynamics may not satisfy the QUAD condition or weak-QUAD condition. First, it gives a new method to analyze the exponential synchronization of the systems. Second, two theorems and their corollaries are proposed for the local or global exponential synchronization of the coupled systems. Finally, an application to the linearly coupled Hopfield neural networks and several simulations are provided for verifying the effectiveness of the theoretical results. (paper)
-
Bethe ansatz and ordinary differential equation correspondence for degenerate Gaudin models
El Araby, Omar; Gritsev, Vladimir; Faribault, Alexandre
2012-03-01
In this work, we generalize the numerical approach to Gaudin models developed earlier by us [Faribault, El Araby, Sträter, and Gritsev, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.83.235124 83, 235124 (2011)] to degenerate systems, showing that their treatment is surprisingly convenient from a numerical point of view. In fact, high degeneracies not only reduce the number of relevant states in the Hilbert space by a non-negligible fraction, they also allow us to write the relevant equations in the form of sparse matrix equations. Moreover, we introduce an inversion method based on a basis of barycentric polynomials that leads to a more stable and efficient root extraction, which most importantly avoids the necessity of working with arbitrary precision. As an example, we show the results of our procedure applied to the Richardson model on a square lattice.
-
Stochastic Differential Equations and Kondratiev Spaces
Energy Technology Data Exchange (ETDEWEB)
Vaage, G.
1995-05-01
The purpose of this mathematical thesis was to improve the understanding of physical processes such as fluid flow in porous media. An example is oil flowing in a reservoir. In the first of five included papers, Hilbert space methods for elliptic boundary value problems are used to prove the existence and uniqueness of a large family of elliptic differential equations with additive noise without using the Hermite transform. The ideas are then extended to the multidimensional case and used to prove existence and uniqueness of solution of the Stokes equations with additive noise. The second paper uses functional analytic methods for partial differential equations and presents a general framework for proving existence and uniqueness of solutions to stochastic partial differential equations with multiplicative noise, for a large family of noises. The methods are applied to equations of elliptic, parabolic as well as hyperbolic type. The framework presented can be extended to the multidimensional case. The third paper shows how the ideas from the second paper can be extended to study the moving boundary value problem associated with the stochastic pressure equation. The fourth paper discusses a set of stochastic differential equations. The fifth paper studies the relationship between the two families of Kondratiev spaces used in the thesis. 102 refs.
-
Spectral theories for linear differential equations
International Nuclear Information System (INIS)
Sell, G.R.
1976-01-01
The use of spectral analysis in the study of linear differential equations with constant coefficients is not only a fundamental technique but also leads to far-reaching consequences in describing the qualitative behaviour of the solutions. The spectral analysis, via the Jordan canonical form, will not only lead to a representation theorem for a basis of solutions, but will also give a rather precise statement of the (exponential) growth rates of various solutions. Various attempts have been made to extend this analysis to linear differential equations with time-varying coefficients. The most complete such extensions is the Floquet theory for equations with periodic coefficients. For time-varying linear differential equations with aperiodic coefficients several authors have attempted to ''extend'' the Foquet theory. The precise meaning of such an extension is itself a problem, and we present here several attempts in this direction that are related to the general problem of extending the spectral analysis of equations with constant coefficients. The main purpose of this paper is to introduce some problems of current research. The primary problem we shall examine occurs in the context of linear differential equations with almost periodic coefficients. We call it ''the Floquet problem''. (author)
-
On implicit abstract neutral nonlinear differential equations
Energy Technology Data Exchange (ETDEWEB)
Hernández, Eduardo, E-mail: lalohm@ffclrp.usp.br [Universidade de São Paulo, Departamento de Computação e Matemática, Faculdade de Filosofia Ciências e Letras de Ribeirão Preto (Brazil); O’Regan, Donal, E-mail: donal.oregan@nuigalway.ie [National University of Ireland, School of Mathematics, Statistics and Applied Mathematics (Ireland)
2016-04-15
In this paper we continue our developments in Hernández and O’Regan (J Funct Anal 261:3457–3481, 2011) on the existence of solutions for abstract neutral differential equations. In particular we extend the results in Hernández and O’Regan (J Funct Anal 261:3457–3481, 2011) for the case of implicit nonlinear neutral equations and we focus on applications to partial “nonlinear” neutral differential equations. Some applications involving partial neutral differential equations are presented.
-
Theory of a higher-order Sturm-Liouville equation
Kozlov, Vladimir
1997-01-01
This book develops a detailed theory of a generalized Sturm-Liouville Equation, which includes conditions of solvability, classes of uniqueness, positivity properties of solutions and Green's functions, asymptotic properties of solutions at infinity. Of independent interest, the higher-order Sturm-Liouville equation also proved to have important applications to differential equations with operator coefficients and elliptic boundary value problems for domains with non-smooth boundaries. The book addresses graduate students and researchers in ordinary and partial differential equations, and is accessible with a standard undergraduate course in real analysis.
-
Analytical exact solution of the non-linear Schroedinger equation
International Nuclear Information System (INIS)
Martins, Alisson Xavier; Rocha Filho, Tarcisio Marciano da
2011-01-01
Full text: In this work we present how to classify and obtain analytical solutions of the Schroedinger equation with a generic non-linearity in 1+1 dimensions. Our approach is based on the determination of Lie symmetry transformation mapping solutions into solutions, and non-classical symmetry transformations, mapping a given solution into itself. From these symmetries it is then possible to reduce the equation to a system of ordinary differential equations which can then be solved using standard methods. The generic non-linearity is handled by considering it as an additional unknown in the determining equations for the symmetry transformations. This results in an over-determined system of non-linear partial differential equations. Its solution can then be determined in some cases by reducing it to the so called involutive (triangular) form, and then solved. This reduction is very tedious and can only performed using a computer algebra system. Once the determining system is solved, we obtain the explicit form for the non-linearity admitting a Lie or non-classical symmetry. The analytical solutions are then derived by solving the reduced ordinary differential equations. The non-linear determining system for the non-classical symmetry transformations and Lie symmetry generators are obtaining using the computer algebra package SADE (symmetry analysis of differential equations), developed at our group. (author)
-
Sparse dynamics for partial differential equations.
Schaeffer, Hayden; Caflisch, Russel; Hauck, Cory D; Osher, Stanley
2013-04-23
We investigate the approximate dynamics of several differential equations when the solutions are restricted to a sparse subset of a given basis. The restriction is enforced at every time step by simply applying soft thresholding to the coefficients of the basis approximation. By reducing or compressing the information needed to represent the solution at every step, only the essential dynamics are represented. In many cases, there are natural bases derived from the differential equations, which promote sparsity. We find that our method successfully reduces the dynamics of convection equations, diffusion equations, weak shocks, and vorticity equations with high-frequency source terms.
-
Yi, Sun; Nelson, Patrick W; Ulsoy, A Galip
2007-04-01
In a turning process modeled using delay differential equations (DDEs), we investigate the stability of the regenerative machine tool chatter problem. An approach using the matrix Lambert W function for the analytical solution to systems of delay differential equations is applied to this problem and compared with the result obtained using a bifurcation analysis. The Lambert W function, known to be useful for solving scalar first-order DDEs, has recently been extended to a matrix Lambert W function approach to solve systems of DDEs. The essential advantages of the matrix Lambert W approach are not only the similarity to the concept of the state transition matrix in lin ear ordinary differential equations, enabling its use for general classes of linear delay differential equations, but also the observation that we need only the principal branch among an infinite number of roots to determine the stability of a system of DDEs. The bifurcation method combined with Sturm sequences provides an algorithm for determining the stability of DDEs without restrictive geometric analysis. With this approach, one can obtain the critical values of delay, which determine the stability of a system and hence the preferred operating spindle speed without chatter. We apply both the matrix Lambert W function and the bifurcation analysis approach to the problem of chatter stability in turning, and compare the results obtained to existing methods. The two new approaches show excellent accuracy and certain other advantages, when compared to traditional graphical, computational and approximate methods.
-
GLOBAL LINEARIZATION OF DIFFERENTIAL EQUATIONS WITH SPECIAL STRUCTURES
Institute of Scientific and Technical Information of China (English)
无
2011-01-01
This paper introduces the global linearization of the differential equations with special structures.The function in the differential equation is unbounded.We prove that the differential equation with unbounded function can be topologically linearlized if it has a special structure.
-
Interactive differential equations modeling program
International Nuclear Information System (INIS)
Rust, B.W.; Mankin, J.B.
1976-01-01
Due to the recent emphasis on mathematical modeling, many ecologists are using mathematics and computers more than ever, and engineers, mathematicians and physical scientists are now included in ecological projects. However, the individual ecologist, with intuitive knowledge of the system, still requires the means to critically examine and adjust system models. An interactive program was developed with the primary goal of allowing an ecologist with minimal experience in either mathematics or computers to develop a system model. It has also been used successfully by systems ecologists, engineers, and mathematicians. This program was written in FORTRAN for the DEC PDP-10, a remote terminal system at Oak Ridge National Laboratory. However, with relatively minor modifications, it can be implemented on any remote terminal system with a FORTRAN IV compiler, or equivalent. This program may be used to simulate any phenomenon which can be described as a system of ordinary differential equations. The program allows the user to interactively change system parameters and/or initial conditions, to interactively select a set of variables to be plotted, and to model discontinuities in the state variables and/or their derivatives. One of the most useful features to the non-computer specialist is the ability to interactively address the system parameters by name and to interactively adjust their values between simulations. These and other features are described in greater detail
-
An approach to rogue waves through the cnoidal equation
Lechuga, Antonio
2014-05-01
Lately it has been realized the importance of rogue waves in some events happening in open seas. Extreme waves and extreme weather could explain some accidents, but not all of them. Every now and then inflicted damages on ships only can be reported to be caused by anomalous and elusive waves, such as rogue waves. That's one of the reason why they continue attracting considerable interest among researchers. In the frame of the Nonlinear Schrödinger equation(NLS), Witham(1974) and Dingemans and Otta (2001)gave asymptotic solutions in moving coordinates that transformed the NLS equation in a ordinary differential equation that is the Duffing or cnoidal wave equation. Applying the Zakharov equation, Stiassnie and Shemer(2004) and Shemer(2010)got also a similar equation. It's well known that this ordinary equation can be solved in elliptic functions. The main aim of this presentation is to sort out the domains of the solutions of this equation, that, of course, are linked to the corresponding solutions of the partial differential equations(PDEs). That being, Lechuga(2007),a simple way to look for anomalous waves as it's the case with some "chaotic" solutions of the Duffing equation.
-
Exact Solution of a Generalized Nonlinear Schrodinger Equation Dimer
DEFF Research Database (Denmark)
Christiansen, Peter Leth; Maniadis, P.; Tsironis, G.P.
1998-01-01
We present exact solutions for a nonlinear dimer system defined throught a discrete nonlinear Schrodinger equation that contains also an integrable Ablowitz-Ladik term. The solutions are obtained throught a transformation that maps the dimer into a double Sine-Gordon like ordinary nonlinear...... differential equation....
-
Lectures on differential Galois theory
Magid, Andy R
1994-01-01
Differential Galois theory studies solutions of differential equations over a differential base field. In much the same way that ordinary Galois theory is the theory of field extensions generated by solutions of (one variable) polynomial equations, differential Galois theory looks at the nature of the differential field extension generated by the solutions of differential equations. An additional feature is that the corresponding differential Galois groups (of automorphisms of the extension fixing the base and commuting with the derivation) are algebraic groups. This book deals with the differential Galois theory of linear homogeneous differential equations, whose differential Galois groups are algebraic matrix groups. In addition to providing a convenient path to Galois theory, this approach also leads to the constructive solution of the inverse problem of differential Galois theory for various classes of algebraic groups. Providing a self-contained development and many explicit examples, this book provides ...
-
Lie algebras and linear differential equations.
Brockett, R. W.; Rahimi, A.
1972-01-01
Certain symmetry properties possessed by the solutions of linear differential equations are examined. For this purpose, some basic ideas from the theory of finite dimensional linear systems are used together with the work of Wei and Norman on the use of Lie algebraic methods in differential equation theory.
-
Introduction to differential equations with dynamical systems
Campbell, Stephen L
2011-01-01
Many textbooks on differential equations are written to be interesting to the teacher rather than the student. Introduction to Differential Equations with Dynamical Systems is directed toward students. This concise and up-to-date textbook addresses the challenges that undergraduate mathematics, engineering, and science students experience during a first course on differential equations. And, while covering all the standard parts of the subject, the book emphasizes linear constant coefficient equations and applications, including the topics essential to engineering students. Stephen Campbell and Richard Haberman--using carefully worded derivations, elementary explanations, and examples, exercises, and figures rather than theorems and proofs--have written a book that makes learning and teaching differential equations easier and more relevant. The book also presents elementary dynamical systems in a unique and flexible way that is suitable for all courses, regardless of length.
-
Low-mode truncation methods in the sine-Gordon equation
International Nuclear Information System (INIS)
Xiong Chuyu.
1991-01-01
In this dissertation, the author studies the chaotic and coherent motions (i.e., low-dimensional chaotic attractor) in some near integrable partial differential equations, particularly the sine-Gordon equation and the nonlinear Schroedinger equation. In order to study the motions, he uses low mode truncation methods to reduce these partial differential equations to some truncated models (low-dimensional ordinary differential equations). By applying many methods available to low-dimensional ordinary differential equations, he can understand the low-dimensional chaotic attractor of PDE's much better. However, there are two important questions one needs to answer: (1) How many modes is good enough for the low mode truncated models to capture the dynamics uniformly? (2) Is the chaotic attractor in a low mode truncated model close to the chaotic attractor in the original PDE? And how close is? He has developed two groups of powerful methods to help to answer these two questions. They are the computation methods of continuation and local bifurcation, and local Lyapunov exponents and Lyapunov exponents. Using these methods, he concludes that the 2N-nls ODE is a good model for the sine-Gordon equation and the nonlinear Schroedinger equation provided one chooses a 'good' basis and uses 'enough' modes (where 'enough' depends on the parameters of the system but is small for the parameter studied here). Therefore, one can use 2N-nls ODE to study the chaos of PDE's in more depth
-
On solutions of variable-order fractional differential equations
Directory of Open Access Journals (Sweden)
Ali Akgül
2017-01-01
solutions to fractional differential equations are compelling to get in real applications, due to the nonlocality and complexity of the fractional differential operators, especially for variable-order fractional differential equations. Therefore, it is significant to enhanced numerical methods for fractional differential equations. In this work, we consider variable-order fractional differential equations by reproducing kernel method. There has been much attention in the use of reproducing kernels for the solutions to many problems in the recent years. We give two examples to demonstrate how efficiently our theory can be implemented in practice.
-
The effect of nonlinearity on unstable zones of Mathieu equation
Indian Academy of Sciences (India)
2017-02-09
Feb 9, 2017 ... Abstract. Mathieu equation is a well-known ordinary differential equation in which the excitation term appears as the non-constant coefficient. The mathematical modelling of many dynamic systems leads to Mathieu equation. The determination of the locus of unstable zone is important for the control of ...
-
Difference and differential equations with applications in queueing theory
Haghighi, Aliakbar Montazer
2013-01-01
A Useful Guide to the Interrelated Areas of Differential Equations, Difference Equations, and Queueing Models Difference and Differential Equations with Applications in Queueing Theory presents the unique connections between the methods and applications of differential equations, difference equations, and Markovian queues. Featuring a comprehensive collection of
-
Tang, Kwong-Tin
2007-01-01
Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.
-
International Nuclear Information System (INIS)
Tang, Bo; He, Yinnian; Wei, Leilei; Zhang, Xindong
2012-01-01
In this Letter, a generalized fractional sub-equation method is proposed for solving fractional differential equations with variable coefficients. Being concise and straightforward, this method is applied to the space–time fractional Gardner equation with variable coefficients. As a result, many exact solutions are obtained including hyperbolic function solutions, trigonometric function solutions and rational solutions. It is shown that the considered method provides a very effective, convenient and powerful mathematical tool for solving many other fractional differential equations in mathematical physics. -- Highlights: ► Study of fractional differential equations with variable coefficients plays a role in applied physical sciences. ► It is shown that the proposed algorithm is effective for solving fractional differential equations with variable coefficients. ► The obtained solutions may give insight into many considerable physical processes.
-
Dynamic Euler-Bernoulli Beam Equation: Classification and Reductions
Directory of Open Access Journals (Sweden)
R. Naz
2015-01-01
Full Text Available We study a dynamic fourth-order Euler-Bernoulli partial differential equation having a constant elastic modulus and area moment of inertia, a variable lineal mass density g(x, and the applied load denoted by f(u, a function of transverse displacement u(t,x. The complete Lie group classification is obtained for different forms of the variable lineal mass density g(x and applied load f(u. The equivalence transformations are constructed to simplify the determining equations for the symmetries. The principal algebra is one-dimensional and it extends to two- and three-dimensional algebras for an arbitrary applied load, general power-law, exponential, and log type of applied loads for different forms of g(x. For the linear applied load case, we obtain an infinite-dimensional Lie algebra. We recover the Lie symmetry classification results discussed in the literature when g(x is constant with variable applied load f(u. For the general power-law and exponential case the group invariant solutions are derived. The similarity transformations reduce the fourth-order partial differential equation to a fourth-order ordinary differential equation. For the power-law applied load case a compatible initial-boundary value problem for the clamped and free end beam cases is formulated. We deduce the fourth-order ordinary differential equation with appropriate initial and boundary conditions.
-
The essential harmony in the classical equations of mathematical physics
Nucci, M C
2001-01-01
The possibility to transform any system of linear ordinary differential equations into a system of constant coefficient equations is demonstrated using Lie theory. Some examples relate the classical equations of mathematical physics to the simple harmonic oscillator. The roles of the third order form of the Ermakov-Pinney equation and of Fleischen-von Weltunter systems are explained.
-
CANONICAL BACKWARD DIFFERENTIATION SCHEMES FOR ...
African Journals Online (AJOL)
This paper describes a new nonlinear backward differentiation schemes for the numerical solution of nonlinear initial value problems of first order ordinary differential equations. The schemes are based on rational interpolation obtained from canonical polynomials. They are A-stable. The test problems show that they give ...
-
First-order partial differential equations
Rhee, Hyun-Ku; Amundson, Neal R
2001-01-01
This first volume of a highly regarded two-volume text is fully usable on its own. After going over some of the preliminaries, the authors discuss mathematical models that yield first-order partial differential equations; motivations, classifications, and some methods of solution; linear and semilinear equations; chromatographic equations with finite rate expressions; homogeneous and nonhomogeneous quasilinear equations; formation and propagation of shocks; conservation equations, weak solutions, and shock layers; nonlinear equations; and variational problems. Exercises appear at the end of mo
-
A Singular Differential Equation Stemming from an Optimal Control Problem in Financial Economics
International Nuclear Information System (INIS)
Brunovský, Pavol; Černý, Aleš; Winkler, Michael
2013-01-01
We consider the ordinary differential equation x 2 u'' = axu'+bu-c(u'-1) 2 , x ∈ (0,x 0 ), with a∈R, b∈R , c>0 and the singular initial condition u(0)=0, which in financial economics describes optimal disposal of an asset in a market with liquidity effects. It is shown in the paper that if a+b 0 then there are infinitely many continuous solutions with indistinguishable asymptotics near 0. Moreover, it is proved that in the latter case there is precisely one solution u corresponding to the choice x 0 =∞ which is such that 0≤u(x)≤x for all x>0, and that this solution is strictly increasing and concave
-
Polygons of differential equations for finding exact solutions
International Nuclear Information System (INIS)
Kudryashov, Nikolai A.; Demina, Maria V.
2007-01-01
A method for finding exact solutions of nonlinear differential equations is presented. Our method is based on the application of polygons corresponding to nonlinear differential equations. It allows one to express exact solutions of the equation studied through solutions of another equation using properties of the basic equation itself. The ideas of power geometry are used and developed. Our approach has a pictorial interpretation, which is illustrative and effective. The method can be also applied for finding transformations between solutions of differential equations. To demonstrate the method application exact solutions of several equations are found. These equations are: the Korteveg-de Vries-Burgers equation, the generalized Kuramoto-Sivashinsky equation, the fourth-order nonlinear evolution equation, the fifth-order Korteveg-de Vries equation, the fifth-order modified Korteveg-de Vries equation and the sixth-order nonlinear evolution equation describing turbulent processes. Some new exact solutions of nonlinear evolution equations are given
-
Reduction operators of Burgers equation.
Pocheketa, Oleksandr A; Popovych, Roman O
2013-02-01
The solution of the problem on reduction operators and nonclassical reductions of the Burgers equation is systematically treated and completed. A new proof of the theorem on the special "no-go" case of regular reduction operators is presented, and the representation of the coefficients of operators in terms of solutions of the initial equation is constructed for this case. All possible nonclassical reductions of the Burgers equation to single ordinary differential equations are exhaustively described. Any Lie reduction of the Burgers equation proves to be equivalent via the Hopf-Cole transformation to a parameterized family of Lie reductions of the linear heat equation.
-
Asymptotic integration of differential and difference equations
Bodine, Sigrun
2015-01-01
This book presents the theory of asymptotic integration for both linear differential and difference equations. This type of asymptotic analysis is based on some fundamental principles by Norman Levinson. While he applied them to a special class of differential equations, subsequent work has shown that the same principles lead to asymptotic results for much wider classes of differential and also difference equations. After discussing asymptotic integration in a unified approach, this book studies how the application of these methods provides several new insights and frequent improvements to results found in earlier literature. It then continues with a brief introduction to the relatively new field of asymptotic integration for dynamic equations on time scales. Asymptotic Integration of Differential and Difference Equations is a self-contained and clearly structured presentation of some of the most important results in asymptotic integration and the techniques used in this field. It will appeal to researchers i...
-
NUMERICAL HOPF BIFURCATION OF DELAY-DIFFERENTIAL EQUATIONS
Institute of Scientific and Technical Information of China (English)
无
2006-01-01
In this paper we consider the numerical solution of some delay differential equations undergoing a Hopf bifurcation. We prove that if the delay differential equations have a Hopf bifurcation point atλ=λ*, then the numerical solution of the equation also has a Hopf bifurcation point atλh =λ* + O(h).
-
Oscillation theory for second order dynamic equations
Agarwal, Ravi P; O''Regan, Donal
2003-01-01
The qualitative theory of dynamic equations is a rapidly developing area of research. In the last 50 years, the Oscillation Theory of ordinary, functional, neutral, partial and impulsive differential equations, and their discrete versions, has inspired many scholars. Hundreds of research papers have been published in every major mathematical journal. Many books deal exclusively with the oscillation of solutions of differential equations, but most of these books appeal only to researchers who already know the subject. In an effort to bring Oscillation Theory to a new and broader audience, the authors present a compact, but thorough, understanding of Oscillation Theory for second order differential equations. They include several examples throughout the text not only to illustrate the theory, but also to provide new direction.
-
Teaching Modeling with Partial Differential Equations: Several Successful Approaches
Myers, Joseph; Trubatch, David; Winkel, Brian
2008-01-01
We discuss the introduction and teaching of partial differential equations (heat and wave equations) via modeling physical phenomena, using a new approach that encompasses constructing difference equations and implementing these in a spreadsheet, numerically solving the partial differential equations using the numerical differential equation…
-
Differential equations I essentials
REA, Editors of
2012-01-01
REA's Essentials provide quick and easy access to critical information in a variety of different fields, ranging from the most basic to the most advanced. As its name implies, these concise, comprehensive study guides summarize the essentials of the field covered. Essentials are helpful when preparing for exams, doing homework and will remain a lasting reference source for students, teachers, and professionals. Differential Equations I covers first- and second-order equations, series solutions, higher-order linear equations, and the Laplace transform.
-
An introduction to differential equations using MATLAB
Butt, Rizwan
2016-01-01
An Introduction to Differential Equations using MATLAB exploits the symbolic, numerical, and graphical capabilitiesof MATLAB to develop a thorough understanding of differential equations algorithms. This book provides the readerwith numerous applications, m-files, and practical examples to problems. Balancing theoretical concepts withcomputational speed and accuracy, the book includes numerous short programs in MATLAB that can be used to solveproblems involving first-and higher-order differential equations, Laplace transforms, linear systems of differentialequations, numerical solutions of differential equations, computer graphics, and more. The author emphasizes thebasic ideas of analytical and numerical techniques and the uses of modern mathematical software (MATLAB) ratherthan relying only on complex mathematical derivations to engineers, mathematician, computer scientists, andphysicists or for use as a textbook in applied or computational courses.A CD-ROM with all the figures, codes, solutions, appendices...
-
Partial differential equations
Agranovich, M S
2002-01-01
Mark Vishik's Partial Differential Equations seminar held at Moscow State University was one of the world's leading seminars in PDEs for over 40 years. This book celebrates Vishik's eightieth birthday. It comprises new results and survey papers written by many renowned specialists who actively participated over the years in Vishik's seminars. Contributions include original developments and methods in PDEs and related fields, such as mathematical physics, tomography, and symplectic geometry. Papers discuss linear and nonlinear equations, particularly linear elliptic problems in angles and gener
-
Optimal overlapping of waveform relaxation method for linear differential equations
International Nuclear Information System (INIS)
Yamada, Susumu; Ozawa, Kazufumi
2000-01-01
Waveform relaxation (WR) method is extremely suitable for solving large systems of ordinary differential equations (ODEs) on parallel computers, but the convergence of the method is generally slow. In order to accelerate the convergence, the methods which decouple the system into many subsystems with overlaps some of the components between the adjacent subsystems have been proposed. The methods, in general, converge much faster than the ones without overlapping, but the computational cost per iteration becomes larger due to the increase of the dimension of each subsystem. In this research, the convergence of the WR method for solving constant coefficients linear ODEs is investigated and the strategy to determine the number of overlapped components which minimizes the cost of the parallel computations is proposed. Numerical experiments on an SR2201 parallel computer show that the estimated number of the overlapped components by the proposed strategy is reasonable. (author)
-
Partial Differential Equations Modeling and Numerical Simulation
Glowinski, Roland
2008-01-01
This book is dedicated to Olivier Pironneau. For more than 250 years partial differential equations have been clearly the most important tool available to mankind in order to understand a large variety of phenomena, natural at first and then those originating from human activity and technological development. Mechanics, physics and their engineering applications were the first to benefit from the impact of partial differential equations on modeling and design, but a little less than a century ago the Schrödinger equation was the key opening the door to the application of partial differential equations to quantum chemistry, for small atomic and molecular systems at first, but then for systems of fast growing complexity. Mathematical modeling methods based on partial differential equations form an important part of contemporary science and are widely used in engineering and scientific applications. In this book several experts in this field present their latest results and discuss trends in the numerical analy...
-
An accurate solution of parabolic equations by expansion in ultraspherical polynomials
International Nuclear Information System (INIS)
Doha, E.H.
1986-11-01
An ultraspherical expansion technique is applied to obtain numerically the solution of the third boundary value problem for linear parabolic partial differential equation in one-space variable. The differential equation with its boundary and initial conditions is reduced to a system of ordinary differential equations for the coefficients of the expansion. This system may be solved analytically or numerically in a step-by-step manner. The method in its present form may be considered as a generalization of that of Dew and Scraton. The extension of the method to the polar-type equations is also considered. (author). 12 refs, 1 tab
-
Hardouin, Charlotte; Minchenko, Andrei; Ovchinnikov, Alexey
2015-01-01
The main motivation of our work is to create an efficient algorithm that decides hypertranscendence of solutions of linear differential equations, via the parameterized differential and Galois theories. To achieve this, we expand the representation theory of linear differential algebraic groups and develop new algorithms that calculate unipotent radicals of parameterized differential Galois groups for differential equations whose coefficients are rational functions. P. Berman and M.F. Singer ...
-
Exact solutions to some modified sine-Gordon equations
International Nuclear Information System (INIS)
Saermark, K.
1983-01-01
Exact, translational solutions to a number of modified sine-Gordon equations are presented. In deriving the equations and the solutions use is made of results from the theory of ordinary differential equations without moving critical points as given by Ince. It is found that kink-like solutions exist also in cases where the coefficients of the trigonometric terms are space- and time-dependent. (Auth.)
-
Lyapunov functionals and stability of stochastic functional differential equations
Shaikhet, Leonid
2013-01-01
Stability conditions for functional differential equations can be obtained using Lyapunov functionals. Lyapunov Functionals and Stability of Stochastic Functional Differential Equations describes the general method of construction of Lyapunov functionals to investigate the stability of differential equations with delays. This work continues and complements the author’s previous book Lyapunov Functionals and Stability of Stochastic Difference Equations, where this method is described for discrete- and continuous-time difference equations. The text begins with a description of the peculiarities of deterministic and stochastic functional differential equations. There follow basic definitions for stability theory of stochastic hereditary systems, and a formal procedure of Lyapunov functionals construction is presented. Stability investigation is conducted for stochastic linear and nonlinear differential equations with constant and distributed delays. The proposed method is used for stability investigation of di...
-
International Nuclear Information System (INIS)
Chowdhury, A.R.; Mukherjee, R.
1984-01-01
The authors have made an exhaustive analysis for an equation introduced by Sabatier (1981) which in the special case reduces to the Harry-Dym equation. First they have deduced the Lie point symmetries and the corresponding ordinary differential equation, through the similarity forms. Next the extended Lie-Backlund type generators are deduced. In the second part the cnoidal wave like solutions are considered. From the Fourier spectrum analysis it is shown that a cnoidal wave breaks into several ordinary solitary waves. (Auth.)
-
A generalized simplest equation method and its application to the Boussinesq-Burgers equation.
Sudao, Bilige; Wang, Xiaomin
2015-01-01
In this paper, a generalized simplest equation method is proposed to seek exact solutions of nonlinear evolution equations (NLEEs). In the method, we chose a solution expression with a variable coefficient and a variable coefficient ordinary differential auxiliary equation. This method can yield a Bäcklund transformation between NLEEs and a related constraint equation. By dealing with the constraint equation, we can derive infinite number of exact solutions for NLEEs. These solutions include the traveling wave solutions, non-traveling wave solutions, multi-soliton solutions, rational solutions, and other types of solutions. As applications, we obtained wide classes of exact solutions for the Boussinesq-Burgers equation by using the generalized simplest equation method.
-
Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em
2017-12-01
Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.