Solving Ordinary Differential Equations [chapter]

Svein Linge, Hans Petter Langtangen
2016 Programming for Computations - MATLAB/Octave  
Differential equations constitute one of the most powerful mathematical tools to understand and predict the behavior of dynamical systems in nature, engineering, and society. A dynamical system is some system with some state, usually expressed by a set of variables, that evolves in time. For example, an oscillating pendulum, the spreading of a disease, and the weather are examples of dynamical systems. We can use basic laws of physics, or plain intuition, to express mathematical rules that
more » ... n the evolution of the system in time. These rules take the form of differential equations. You are probably well experienced with equations, at least equations like ax Cb D 0 or ax 2 Cbx Cc D 0. Such equations are known as algebraic equations, and the unknown is a number. The unknown in a differential equation is a function, and a differential equation will almost always involve this function and one or more derivatives of the function. For example, f 0 .x/ D f .x/ is a simple differential equation (asking if there is any function f such that it equals its derivative -you might remember that e x is a candidate). The present chapter starts with explaining how easy it is to solve both single (scalar) first-order ordinary differential equations and systems of first-order differential equations by the Forward Euler method. We demonstrate all the mathematical and programming details through two specific applications: population growth and spreading of diseases. Then we turn to a physical application: oscillating mechanical systems, which arise in a wide range of engineering situations. The differential equation is now of second order, and the Forward Euler method does not perform well. This observation motivates the need for other solution methods, and we derive the Euler-Cromer 87
doi:10.1007/978-3-319-32452-4_4 fatcat:f45zmkhfkvd43esoasukvghkki