See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/328281723

Dynamical Systems with Applications using Python (Information)

Book · October 2018

CITATION READS

1 4,222

1 author:

Stephen Lynch
Manchester Metropolitan University
67 PUBLICATIONS   1,256 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Dynamical Systems with Applications using Python View project

Dynamical Systems with Applications using Mathematica, 2nd Edition View project

All content following this page was uploaded by Stephen Lynch on 15 October 2018.

The user has requested enhancement of the downloaded file.

 birkhauser-science.de

Stephen Lynch

Dynamical Systems with

Applications using Python
Designed for a broad audience of students in applied mathematics, physics,
and engineering
Represents dynamical systems with popular Python libraries like sympy,
numpy, and matplotlib
Explores a variety of advanced topics in dynamical systems, like neural
networks, fractals, and nonlinear optics, at an undergraduate level
This textbook provides a broad introduction to continuous and discrete dynamical systems. With
its hands-on approach, the text leads the reader from basic theory to recently published
research material in nonlinear ordinary differential equations, nonlinear optics, multifractals,
neural networks, and binary oscillator computing. Dynamical Systems with Applications Using
Python takes advantage of Python’s extensive visualization, simulation, and algorithmic tools to
study those topics in nonlinear dynamical systems through numerical algorithms and
1st ed. 2018, XVI, 665 p. 277 illus., 118
generated diagrams. After a tutorial introduction to Python, the first part of the book deals with
illus. in color.
continuous systems using differential equations, including both ordinary and delay differential
equations. The second part of the book deals with discrete dynamical systems and progresses
Printed book
to the study of both continuous and discrete systems in contexts like chaos control and
Hardcover
synchronization, neural networks, and binary oscillator computing. These later sections are
79,99 € | £69.99 | $89.99
useful reference material for undergraduate student projects. The book is rounded off with
[1]
85,59 € (D) | 87,99 € (A) | CHF
example coursework to challenge students’ programming abilities and Python-based exam
94,50
questions. This book will appeal to advanced undergraduate and graduate students, applied
eBook mathematicians, engineers, and researchers in a range of disciplines, such as biology,
67,82 € | £55.99 | $69.99 chemistry, computing, economics, and physics.
[2]
67,82 € (D) | 67,82 € (A) | CHF
75,50
Available from your library or
springer.com/shop
MyCopy [3]

Printed eBook for just

€ | $ 24.99
springer.com/mycopy

Order online at springer.com / or for the Americas call (toll free) 1-800-SPRINGER /
Lifelong 40% discount for authors or email us at: [email protected]. / For outside the Americas call +49 (0) 6221-345-4301 /
or email us at: [email protected].
The first € price and the £ and $ price are net prices, subject to local VAT. Prices indicated with [1] include VAT for
books; the €(D) includes 7% for Germany, the €(A) includes 10% for Austria. Prices indicated with [2] include VAT for
electronic products; 19% for Germany, 20% for Austria. All prices exclusive of carriage charges. Prices and other
details are subject to change without notice. All errors and omissions excepted. [3] No discount for MyCopy.

Part of
Stephen Lynch

Dynamical Systems with

Applications using Python
Stephen Lynch
Manchester Metropolitan University
Manchester, UK

ISBN 978-3-319-78144-0 ISBN 978-3-319-78145-7 (eBook)

https://doi.org/10.1007/978-3-319-78145-7

Library of Congress Control Number: 2018952351

Mathematics Subject Classification (2010): 00A05, 00A69, 37-01, 34-01, 34Cxx, 34Dxx, 34H10, 34K18,
37-04, 37Nxx, 68U10, 78A60, 92B20

© Springer International Publishing AG, part of Springer Nature 2018

This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of
the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation,
broadcasting, reproduction on microfilms or in any other physical way, and transmission or information
storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology
now known or hereafter developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this publica-
tion does not imply, even in the absence of a specific statement, that such names are exempt from the
relevant protective laws and regulations and therefore free for general use.
The publisher, the authors and the editors are safe to assume that the advice and information in this book
are believed to be true and accurate at the date of publication. Neither the publisher nor the authors or
the editors give a warranty, express or implied, with respect to the material contained herein or for any
errors or omissions that may have been made. The publisher remains neutral with regard to jurisdictional
claims in published maps and institutional affiliations.

This book is published under the imprint Birkhäuser, www.birkhauser-science.com by the registered
company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface

This book provides an introduction to the theory of dynamical systems with

the aid of Python. It is written for both senior undergraduates and graduate
students. Chapter 1 provides a tutorial introduction to Python—new users
should go through this chapter carefully while those moderately familiar and
experienced users will find this chapter a useful source of reference. The first
part of the book deals with continuous systems using differential equations,
including both ordinary and delay differential equations (Chapters 2–12), the
second part is devoted to the study of discrete systems (Chapters 13–17), and
Chapters 18–21 deal with both continuous and discrete systems. Chapter 22
gives examples of coursework and also lists three Python-based examinations
to be sat in a computer laboratory with access to Python. Chapter 23 lists
answers to all of the exercises given in the book. It should be pointed out that
dynamical systems theory is not limited to these topics but also encompasses
partial differential equations, integral and integro-differential equations, and
stochastic systems, for instance. References [1–6] given at the end of the Pref-
ace provide more information for the interested reader. The author has gone
for breadth of coverage rather than fine detail and theorems with proofs are
kept at a minimum. The material is not clouded by functional analytic and
group theoretical definitions, and so is intelligible to readers with a general
mathematical background. Some of the topics covered are scarcely covered
elsewhere. Most of the material in Chapters 9–12 and 16–21 is at postgrad-
uate level and has been influenced by the author’s own research interests.
There is more theory in these chapters than in the rest of the book since it is
not easily accessed anywhere else. It has been found that these chapters are
especially useful as reference material for senior undergraduate project work.
The theory in other chapters of the book is dealt with more comprehensively
in other texts, some of which may be found in the references section of the
corresponding chapter. The book has a very hands-on approach and takes
the reader from the basic theory right through to recently published research
material.

v
vi Preface

Python is extremely popular with a wide range of researchers from all

sorts of disciplines; it has a very user-friendly interface and has extensive
visualization and numerical computation capabilities. It is an ideal package to
adopt for the study of nonlinear dynamical systems; the numerical algorithms
work very quickly, and complex pictures can be plotted within seconds.
The first chapter provides an efficient tutorial introduction to Python.
Simple Python programming is introduced using three basic programming
structures: defining functions, for loops, and if, then, else constructs. New
users will find the tutorials will enable them to become familiar with Python
within a few days. Both engineering and mathematics students appreciate
this method of teaching and I have found that it generally works well with
one staff member to about twenty students in a computer laboratory. In most
cases, I have chosen to list the Python commands at the end of each chapter;
this avoids unnecessary cluttering in the text. The Python programs have
been kept as simple as possible and should run under later versions of the
package. All Python files for the book (including updates and extra files) can
even be downloaded from the Web via GitHub at:

https://github.com/springer-math/dynamical-systems-with-applications-
using-python

Readers will find that they can reproduce the figures given in the text, and
then it is not too difficult to change parameters or equations to investigate
other systems.
Chapters 2–12 deal with continuous dynamical systems. Chapters 2 and 3
cover some theory of ordinary differential equations and applications to mod-
els in the real world are given. The theory of differential equations applied
to chemical kinetics and electric circuits is introduced in some detail. The
memristor is introduced and one of the most remarkable stories in the history
of mathematics is relayed. Chapter 2 ends with the existence and uniqueness
theorem for the solutions of certain types of differential equations. The theory
behind the construction of phase plane portraits for two-dimensional systems
is dealt with in Chapter 3. Applications are taken from chemical kinetics, eco-
nomics, electronics, epidemiology, mechanics, and population dynamics. The
modeling of the populations of interacting species is discussed in some detail
in Chapter 4 and domains of stability are discussed for the first time. Limit
cycles, or isolated periodic solutions, are introduced in Chapter 5. Since we
live in a periodic world, these are the most common type of solution found
when modeling nonlinear dynamical systems. They appear extensively when
modeling both the technological and natural sciences. Hamiltonian, or con-
servative, systems and stability are discussed in Chapter 6, and Chapter 7
is concerned with how planar systems vary depending upon a parameter.
Bifurcation, bistability, multistability, and normal forms are discussed.
Preface vii

The reader is first introduced to the concept of chaos in continuous sys-

tems in Chapters 8 and 9, where three-dimensional systems and Poincaré
maps are investigated. These higher-dimensional systems can exhibit strange
attractors and chaotic dynamics. One can rotate the three-dimensional ob-
jects in Python and plot time series plots to get a better understanding of
the dynamics involved. Once again, the theory can be applied to chemi-
cal kinetics (including stiff systems), electric circuits, and epidemiology; a
simplified model for the weather is also briefly discussed. Chapter 9 deals
with Poincaré first return maps that can be used to untangle complicated
interlacing trajectories in higher-dimensional spaces. A periodically driven
nonlinear pendulum is also investigated by means of a nonautonomous dif-
ferential equation. Both local and global bifurcations are investigated in
Chapter 10. The main results and statement of the famous second part of
David Hilbert’s sixteenth problem are listed in Chapter 11. In order to un-
derstand these results, Poincaré compactification is introduced. The study
of continuous systems ends with one of the authors specialities—limit cycles
of Liénard systems. There is some detail on Liénard systems, in particular,
in this part of the book, but they do have a ubiquity for systems in the
plane. Chapter 12 provides an introduction to delay differential equations
with applications in biology and nonlinear optics.
Chapters 13–17 deal with discrete dynamical systems. Chapter 13 starts
with a general introduction to iteration and linear recurrence (or difference)
equations. The bulk of the chapter is concerned with the Leslie model used
to investigate the population of a single species split into different age classes.
Harvesting and culling policies are then investigated and optimal solutions are
sought. Nonlinear discrete dynamical systems are dealt with in Chapter 14.
Bifurcation diagrams, chaos, intermittency, Lyapunov exponents, periodic-
ity, quasiperiodicity, and universality are some of the topics introduced. The
theory is then applied to real-world problems from a broad range of disci-
plines including population dynamics, biology, economics, nonlinear optics,
and neural networks. Chapter 15 is concerned with complex iterative maps
in the Argand plane, where Julia sets and the now-famous Mandelbrot set
are plotted. Basins of attraction are investigated for these complex systems
and Newton fractals are introduced. As a simple introduction to optics, elec-
tromagnetic waves and Maxwell’s equations are studied at the beginning of
Chapter 16. Complex iterative equations are used to model the propagation
of light waves through nonlinear optical fibers. A brief history of nonlinear
bistable optical resonators is discussed and the simple fiber ring resonator
is analyzed in particular. Chapter 16 is devoted to the study of these op-
tical resonators, and there is discussion on phenomena such as bistability,
chaotic attractors, feedback, hysteresis, instability, linear stability analysis,
multistability, nonlinearity, and steady states. The first and second iterative
methods are defined in this chapter. Some simple fractals may be constructed
viii Preface

using pencil and paper in Chapter 17, and the concept of fractal dimension
is introduced. Fractals may be thought of as identical motifs repeated on
ever-reduced scales. Unfortunately, most of the fractals appearing in nature
are not homogeneous but are more heterogeneous, hence the need for the
multifractal theory given later in the chapter. It has been found that the
distribution of stars and galaxies in our universe is multifractal, and there
is even evidence of multifractals in rainfall, stock markets, and heartbeat
rhythms. Applications in geoscience, materials science, microbiology, and
image processing are briefly discussed. Chapter 18 provides a brief introduc-
tion to image processing which is being used more and more by a diverse
range of scientific disciplines, especially medical imaging. The fast Fourier
transform is introduced and has a wide range of applications throughout the
realms of science.
Chapter 19 is devoted to the new and exciting theory behind chaos con-
trol and synchronization. For most systems, the maxim used by engineers
in the past has been “stability good, chaos bad,” but more and more nowa-
days this is being replaced with “stability good, chaos better.” There are
exciting and novel applications in cardiology, communications, engineering,
laser technology, and space research, for example. A brief introduction to
the enticing field of neural networks is presented in Chapter 20. Imagine
trying to make a computer mimic the human brain. One could ask the ques-
tion: In the future will it be possible for computers to think and even be
conscious? The human brain will always be more powerful than traditional,
sequential, logic-based digital computers and scientists are trying to incor-
porate some features of the brain into modern computing. Neural networks
perform through learning and no underlying equations are required. Mathe-
maticians and computer scientists are attempting to mimic the way neurons
work together via synapses; indeed, a neural network can be thought of as a
crude multidimensional model of the human brain. The expectations are high
for future applications in a broad range of disciplines. Neural networks are
already being used in machine learning and pattern recognition (computer vi-
sion, credit card fraud, prediction and forecasting, disease recognition, facial
and speech recognition), the consumer home entertainment market, psycho-
logical profiling, predicting wave over-topping events, and control problems,
for example. They also provide a parallel architecture allowing for very fast
computational and response times. In recent years, the disciplines of neu-
ral networks and nonlinear dynamics have increasingly coalesced and a new
branch of science called neurodynamics is emerging. Lyapunov functions can
be used to determine the stability of certain types of neural network. There
is also evidence of chaos, feedback, nonlinearity, periodicity, and chaos syn-
chronization in the brain.
Chapter 21 focuses on binary oscillator computing, the subject of UK,
International, and Taiwanese patents. The author and his co-inventor, Jon
Preface ix

Borresen, came up with the idea when modeling connected biological neu-
rons. Binary oscillator technology can be applied to the design of arithmetic
logic units (ALUs), memory, and other basic computing components. It
has the potential to provide revolutionary computational speed-up, energy
saving, and novel applications and may be applicable to a variety of techno-
logical paradigms including biological neurons, complementary metal-oxide-
semiconductor (CMOS), memristors, optical oscillators, and superconducting
materials. The research has the potential for MMU and industrial partners
to develop super fast, low-power computers and may provide an assay for
neuronal degradation for brain malfunctions such as Alzheimer’s, epilepsy,
and Parkinson’s disease!
Examples of coursework and three examination-type papers are listed in
Chapter 22, and a complete set of solutions for the book is listed in Chap-
ter 23.
Both textbooks and research papers are presented in the list of references.
The textbooks can be used to gain more background material, and the re-
search papers have been given to encourage further reading and independent
study.
This book is informed by the research interests of the author, which are
currently nonlinear ordinary differential equations, nonlinear optics, multi-
fractals, neural networks, and binary oscillator computing. Some references
include recently published research articles by the author along with two
patents.
The prerequisites for studying dynamical systems using this book are
undergraduate courses in linear algebra, real and complex analysis, calculus,
and ordinary differential equations; a knowledge of a computer language such
as Basic, C, or Fortran would be beneficial but not essential.

Recommended Textbooks
[1] H.P Langtangen and A. Logg, Solving PDEs in Python: The FEniCS
Tutorial I (Simula SpringerBriefs on Computing), Springer, New York, 2017.

[2] B. Bhattacharya and M. Majumdar, Random Dynamical Systems in Fi-

nance, Chapman & Hall/CRC, New York, 2016.

[3] L.C. de Barros, R.C. Bassanezi and W.A. Lodwick, A First Course in
Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics: Theory and
Applications, Springer, New York, 2016.

[4] V. Volterra, Theory of Functionals and of Integral and Integro-Differential

Equations, Dover Publications, New York, 2005.
x Preface

[5] J. Mallet-Paret (Editor), J. Wu (Editor), H. Zhu (Editor), Y. Yi (Editor),

Infinite Dimensional Dynamical Systems (Fields Institute Communications),
Springer, New York, 2013.

[6] C. Bernido, M.V. Carpio-Bernido, M. Grothaus et al., Stochastic and In-

finite Dimensional Analysis, Birkhäuser, New York, 2016.

Special thanks go to Ben Nuttall (Python guru), Community Manager,

the Raspberry Pi Foundation, Cambridge, UK (www.raspberrypi.org), for
reviewing this book. I would also like to express my sincere thanks to all of the
reviewers of this book and the other editions of my books. As always, thanks
also go to Birkhäuser and Springer, especially Samuel DiBella (Assistant
Editor, Springer Nature). Finally, thanks to my family and especially my
wife Gaynor, and our children, Sebastian and Thalia, for their continuing
love, inspiration, and support.

Manchester, UK Stephen Lynch FIMA SFHEA

Contents

1 A Tutorial Introduction to Python 1

1.1 The IDLE Integrated Development Environment for Python 2
1.1.1 Tutorial One: Using Python as a Powerful Calculator 4
1.1.2 Tutorial Two: Simple Programming with Python . . . 6
1.1.3 Tutorial Three: Simple Plotting Using
the Turtle Module . . . . . . . . . . . . . . . . . . . . 9
1.2 Anaconda, Spyder and the Libraries,
Sympy, Numpy, and Matplotlib . . . . . . . . . . . . . . . . 14
1.2.1 Tutorial One: A Tutorial Introduction to Sympy . . . 15
1.2.2 Tutorial Two: A Tutorial Introduction to Numpy and
Matplotlib . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.2.3 Tutorial Three: Simple Programming, Solving
ODEs, and More Detailed Plots . . . . . . . . . . . . . 20
1.3 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

2 Differential Equations 33
2.1 Simple Differential Equations and Applications . . . . . . . 34
2.1.1 Linear Differential Equations . . . . . . . . . . . . . . 34
2.1.2 Separable Differential Equations . . . . . . . . . . . . 35
2.1.3 Exact Differential Equations . . . . . . . . . . . . . . 39
2.1.4 Homogeneous Differential Equations . . . . . . . . . . 40
2.2 Applications to Chemical Kinetics . . . . . . . . . . . . . . . 44
2.3 Applications to Electric Circuits . . . . . . . . . . . . . . . . 48
2.4 Existence and Uniqueness Theorem . . . . . . . . . . . . . . 53
2.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 57
2.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59

3 Planar Systems 65
3.1 Canonical Forms . . . . . . . . . . . . . . . . . . . . . . . . . 66
3.1.1 Real Distinct Eigenvalues . . . . . . . . . . . . . . . . 68
3.1.2 Complex Eigenvalues (λ = α ± iβ) . . . . . . . . . . . 69
3.1.3 Repeated Real Eigenvalues . . . . . . . . . . . . . . . 70

xi
xii Contents

3.2 Eigenvectors Defining Stable and Unstable Manifolds . . . . 72

3.3 Phase Portraits of Linear Systems in the Plane . . . . . . . 74
3.4 Linearization and Hartman’s Theorem . . . . . . . . . . . . 78
3.5 Constructing Phase Plane Diagrams . . . . . . . . . . . . . . 79
3.6 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 87
3.7 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90

4 Interacting Species 95
4.1 Competing Species . . . . . . . . . . . . . . . . . . . . . . . 96
4.2 Predator-Prey Models . . . . . . . . . . . . . . . . . . . . . . 99
4.3 Other Characteristics Affecting Interacting Species . . . . . 104
4.4 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 107
4.5 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108

5 Limit Cycles 113

5.1 Historical Background . . . . . . . . . . . . . . . . . . . . . 114
5.2 Existence and Uniqueness of Limit Cycles in the Plane . . . 117
5.3 Nonexistence of Limit Cycles in the Plane . . . . . . . . . . 123
5.4 Perturbation Methods . . . . . . . . . . . . . . . . . . . . . 127
5.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 136
5.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139

6 Hamiltonian Systems, Lyapunov Functions, and Stability 145

6.1 Hamiltonian Systems in the Plane . . . . . . . . . . . . . . . 146
6.2 Lyapunov Functions and Stability . . . . . . . . . . . . . . . 151
6.3 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 157
6.4 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 158

7 Bifurcation Theory 163

7.1 Bifurcations of Nonlinear Systems in the Plane . . . . . . . . 164
7.1.1 A Saddle-Node Bifurcation . . . . . . . . . . . . . . . 165
7.1.2 A Transcritical Bifurcation . . . . . . . . . . . . . . . 167
7.1.3 A Pitchfork Bifurcation . . . . . . . . . . . . . . . . . 167
7.1.4 A Hopf Bifurcation . . . . . . . . . . . . . . . . . . . . 170
7.2 Normal Forms . . . . . . . . . . . . . . . . . . . . . . . . . . 170
7.3 Multistability and Bistability . . . . . . . . . . . . . . . . . . 174
7.4 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 178
7.5 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180

8 Three-Dimensional Autonomous Systems and Chaos 185

8.1 Linear Systems and Canonical Forms . . . . . . . . . . . . . 186
8.2 Nonlinear Systems and Stability . . . . . . . . . . . . . . . . 190
8.3 The Rössler System and Chaos . . . . . . . . . . . . . . . . 194
Contents xiii

8.3.1 The Rössler Attractor . . . . . . . . . . . . . . . . . . 194

8.3.2 Chaos . . . . . . . . . . . . . . . . . . . . . . . . . . . 196
8.4 The Lorenz Equations, Chua’s Circuit, and the
Belousov-Zhabotinski Reaction . . . . . . . . . . . . . . . . . 199
8.4.1 The Lorenz Equations . . . . . . . . . . . . . . . . . . 199
8.4.2 Chua’s Circuit . . . . . . . . . . . . . . . . . . . . . . 201
8.4.3 The Belousov-Zhabotinski (BZ) Reaction . . . . . . . 204
8.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 207
8.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 211

9 Poincaré Maps and Nonautonomous Systems in the Plane 215

9.1 Poincaré Maps . . . . . . . . . . . . . . . . . . . . . . . . . . 216
9.2 Hamiltonian Systems with Two Degrees of Freedom . . . . . 221
9.3 Nonautonomous Systems in the Plane . . . . . . . . . . . . . 227
9.4 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 235
9.5 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 239

10 Local and Global Bifurcations 245

10.1 Small-Amplitude Limit Cycle
Bifurcations . . . . . . . . . . . . . . . . . . . . . . . . . . . 246
10.2 Gröbner Bases . . . . . . . . . . . . . . . . . . . . . . . . . . 252
10.3 Melnikov Integrals and Bifurcating Limit Cycles
from a Center . . . . . . . . . . . . . . . . . . . . . . . . . . 258
10.4 Bifurcations Involving Homoclinic Loops . . . . . . . . . . . 260
10.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 262
10.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 266

11 The Second Part of Hilbert’s Sixteenth Problem 271

11.1 Statement of Problem and Main Results . . . . . . . . . . . 272
11.2 Poincaré Compactification . . . . . . . . . . . . . . . . . . . 275
11.3 Global Results for Liénard Systems . . . . . . . . . . . . . . 281
11.4 Local Results for Liénard Systems . . . . . . . . . . . . . . . 290
11.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 291
11.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 292

12 Delay Differential Equations 297

12.1 Introduction and the Method of Steps . . . . . . . . . . . . . 298
12.2 Applications in Biology . . . . . . . . . . . . . . . . . . . . . 304
12.3 Applications in Nonlinear Optics . . . . . . . . . . . . . . . 310
12.4 Other Applications . . . . . . . . . . . . . . . . . . . . . . . 313
12.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 315
12.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 320
xiv Contents

13 Linear Discrete Dynamical Systems 327

13.1 Recurrence Relations . . . . . . . . . . . . . . . . . . . . . . 328
13.2 The Leslie Model . . . . . . . . . . . . . . . . . . . . . . . . 333
13.3 Harvesting and Culling Policies . . . . . . . . . . . . . . . . 337
13.4 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 342
13.5 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 343

14 Nonlinear Discrete Dynamical Systems 347

14.1 The Tent Map and Graphical Iterations . . . . . . . . . . . . 348
14.2 Fixed Points and Periodic Orbits . . . . . . . . . . . . . . . 353
14.3 The Logistic Map, Bifurcation Diagram,
and Feigenbaum Number . . . . . . . . . . . . . . . . . . . . 360
14.4 Gaussian and Hénon Maps . . . . . . . . . . . . . . . . . . . 368
14.5 Applications . . . . . . . . . . . . . . . . . . . . . . . . . . . 373
14.6 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 376
14.7 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 380

15 Complex Iterative Maps 385

15.1 Julia Sets and the Mandelbrot Set . . . . . . . . . . . . . . . 386
15.2 Boundaries of Periodic Orbits . . . . . . . . . . . . . . . . . 391
15.3 The Newton Fractal . . . . . . . . . . . . . . . . . . . . . . . 395
15.4 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 396
15.5 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 399

16 Electromagnetic Waves and Optical Resonators 403

16.1 Maxwell’s Equations and Electromagnetic Waves . . . . . . 404
16.2 Historical Background . . . . . . . . . . . . . . . . . . . . . 406
16.3 The Nonlinear SFR Resonator . . . . . . . . . . . . . . . . . 412
16.4 Chaotic Attractors and Bistability . . . . . . . . . . . . . . . 413
16.5 Linear Stability Analysis . . . . . . . . . . . . . . . . . . . . 417
16.6 Instabilities and Bistability . . . . . . . . . . . . . . . . . . . 420
16.7 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 424
16.8 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 428

17 Fractals and Multifractals 433

17.1 Construction of Simple Examples . . . . . . . . . . . . . . . 434
17.2 Calculating Fractal Dimensions . . . . . . . . . . . . . . . . 441
17.3 A Multifractal Formalism . . . . . . . . . . . . . . . . . . . . 448
17.4 Multifractals in the Real World and Some Simple Examples 452
17.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 459
17.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 464
Contents xv

18 Image Processing with Python 471

18.1 Image Processing and Matrices . . . . . . . . . . . . . . . . 472
18.2 The Fast Fourier Transform . . . . . . . . . . . . . . . . . . 477
18.3 The Fast Fourier Transform on Images . . . . . . . . . . . . 484
18.4 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 487

19 Chaos Control and Synchronization 491

19.1 Historical Background . . . . . . . . . . . . . . . . . . . . . 492
19.2 Controlling Chaos in the Logistic Map . . . . . . . . . . . . 497
19.3 Controlling Chaos in the Hénon Map . . . . . . . . . . . . . 498
19.4 Chaos Synchronization . . . . . . . . . . . . . . . . . . . . . 505
19.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 509
19.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 513

20 Neural Networks 519

20.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . 520
20.2 The Delta Learning Rule and Backpropagation . . . . . . . 526
20.3 The Hopfield Network and Lyapunov Stability . . . . . . . . 531
20.4 Neurodynamics . . . . . . . . . . . . . . . . . . . . . . . . . 541
20.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 545
20.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 550

21 Binary Oscillator Computing 557

21.1 Brain Inspired Computing . . . . . . . . . . . . . . . . . . . 558
21.2 Oscillatory Threshold Logic . . . . . . . . . . . . . . . . . . 563
21.3 Applications and Future Work . . . . . . . . . . . . . . . . . 571
21.4 An Assay for Neuronal Degradation . . . . . . . . . . . . . . 577
21.5 Python Programs . . . . . . . . . . . . . . . . . . . . . . . . 579
21.6 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 584

22 Coursework and Examination-Type Questions 591

22.1 Examples of Coursework Questions . . . . . . . . . . . . . . 592
22.2 Examination 1 . . . . . . . . . . . . . . . . . . . . . . . . . . 607
22.3 Examination 2 . . . . . . . . . . . . . . . . . . . . . . . . . . 610
22.4 Examination 3 . . . . . . . . . . . . . . . . . . . . . . . . . . 613

23 Solutions to Exercises 619

23.1 Chapter 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 619
23.2 Chapter 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 622
23.3 Chapter 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 623
23.4 Chapter 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 625
23.5 Chapter 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 626
23.6 Chapter 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 628
23.7 Chapter 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 628
xvi Contents

23.8 Chapter 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 630

23.9 Chapter 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 631
23.10 Chapter 10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 631
23.11 Chapter 11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 632
23.12 Chapter 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 634
23.13 Chapter 13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 634
23.14 Chapter 14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 636
23.15 Chapter 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 637
23.16 Chapter 16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 638
23.17 Chapter 17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 638
23.18 Chapter 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 639
23.19 Chapter 19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 639
23.20 Chapter 20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 640
23.21 Chapter 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 640
23.22 Chapter 22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 641

Appendix A Index of Python Programs 645

A.1 IDLE Python Programs . . . . . . . . . . . . . . . . . . . . 645
A.2 Anaconda Python Programs . . . . . . . . . . . . . . . . . . 646

Index 651
Appendix A

Index of Python Programs

Readers can download the Python program files via GitHub:

https://github.com/springer-math/dynamical-systems-with-applications-
using-python

These files will be kept up-to-date and extra files will be added in the
forthcoming years.

A.1 IDLE Python Programs

These files include solutions to the Exercises listed in Chapter 1.

euclid_algorithm.py --- See Exercise 10.

F2C.py --- See Exercise 1(a).
F2K.py --- Converts degrees Fahrenheit to Kelvin.
fibonacci.py --- Lists first n terms of the Fibonacci sequence.
fmu.py --- The logistic function.
fractal_tree.py --- Plots a fractal tree.
fractal_tree_color.py --- Plots a color fractal tree.
grade.py --- Converts a score to a grade.
guess_number.py --- Guess the number game.
koch_snowflake.py --- See Exercise 1(d).
koch_square.py --- Plots a Koch square fractal.
Pythag_Triples.py --- See Exercise 1(c).
sierpinski.py --- Plots a Sierpinski triangle fractal.
sierpinski_square.py --- Plots a Sierpinski square fractal.

© Springer International Publishing AG, part of Springer Nature 2018 645

S. Lynch, Dynamical Systems with Applications using Python,
https://doi.org/10.1007/978-3-319-78145-7
646 Appendix: Springer
c

sum_primes.py --- See Exercise 1(b).

sum_n.py --- Sums the natural numbers to n.

A.2 Anaconda Python Programs

If you have difficulty with the animation programs in Spyder, you have to
change the backend to run an animation in the IPython console. You can do
that by running

In[1]: %matplotlib qt5

before the animation. If you don’t want to use this command every time,
you can go to: Tools, Preferences, IPython Console, Graphics, Backend, and
change it from “Inline” to “Automatic.”

Program_01a.py --- Solve a simple ODE.

Program_01b.py --- Solve a second order ODE.
Program_01c.py --- Plot two curves on one graph.
Program_01d.py --- Subplots.
Program_01e.py --- Surface and contour plot in 3D.
Program_01f.py --- A parametric curve in 3D.
Program_01g.py --- Animation of a simple curve.

Program_02a.py ---
Solve a separable ODE.
Program_02b.py ---
Solve the logistic ODE.
Program_02c.py ---
Power series solution.
Program_02d.py ---
Power series solution for van der Pol.
Program_02e.py ---
Plot series solution against numerical
solution.
Program_02f.py --- Solve a linear first order ODE.
Program_02g.py --- Solve a linear second order ODE.

Program_03a.py --- Plot the phase portrait of a linear system.

Program_03b.py --- Plot the phase portrait of a nonlinear
system.
Program_03c.py --- Finding critical points.

Program_04a.py --- Phase portrait and time series of Holling-

Tanner model.
Appendix 647

Program_05a.py --- Limit cycle of a Fitzhugh-Nagumo system.

Program_05b.py --- Approximate and numerical solutions to ODEs.
Program_05c.py --- Error between one-term and numerical
solution.
Program_05d.py --- Lindstedt-Poincare technique.

Program_06a.py --- Contour plot.

Program_06b.py --- Surface plot.

Program_07a.py --- Animation of a simple curve.

Program_07b.py --- Animation of a subcritical Hopf bifurcation.
Program_07c.py --- Animation of a SNIC bifurcation.

Program_08a.py --- The Rossler attractor.

Program_08b.py --- The Lorenz Attractor.
Program_08c.py --- The Belousov-Zhabotinsky reaction.
Program_08d.py --- Animation of a Chua circuit bifurcation.

Program_09a.py --- Simple Poincare return map.

Program_09b.py --- Hamiltonian with two degrees of freedom plot.
Program_09c.py --- Phase portrait and Poincare map for the
Duffing system.
Program_09d.py --- Bifurcation diagram of Duffing equation.

Program_10a.py --- Computing Lyapunov quantities.

Program_10b.py --- Division algorithm for multivariate
polynomials.
Program_10c.py --- S-polynomial.
Program_10d.py --- Computing the Groebner basis.
Program_10e.py --- Computing Groebner basis of Lyapunov
quantities.
Program_10f.py --- Animation of a homoclinic limit cycle
bifurcation.
Program_10g.py --- Animation of a homoclinic limit cycle
bifurcation.

Program_11a.py --- Animation of a Lienard limit cycle.

Program_12a.py --- The method of steps.

Program_12b.py --- Plot of solution by method of steps.
Program_12c.py --- The Mackey-Glass DDE.
Program_12d.py --- The Lang-Kobayashi DDEs.
648 Appendix: Springer
c

Program_13a.py --- Computing bank interest.

Program_13b.py --- Solving a second order recurrence relation.
Program_13c.py --- The Leslie matrix, eigenvalues and
eigenvectors.

Program_14a.py --- Graphical iteration of the tent map.

Program_14b.py --- Bifurcation diagram of the logistic map.
Program_14c.py --- Computing Lyapunov exponents for the
logistic map.
Program_14d.py --- Iteration of the Henon map.
Program_14e.py --- Lyapunov exponents of the Henon map.

Program_15a.py --- Point plot for a Julia set.

Program_15b.py --- Colormap of a Julia set.
Program_15c.py --- Color Mandelbrot set.
Program_15d.py --- Color Newton fractal Julia set.

Program_16a.py --- Intersection of implicit curves.

Program_16b.py --- Chaotic Attractor of the Ikeda map
Program_16c.py --- Bifurcation diagram of the Ikeda map.

Program_17a.py --- The Koch curve.

Program_17b.py --- Chaos game and the Sierpnski triangle.
Program_17c.py --- Barnsley’s fern.
Program_17d.py --- Subplots of tau, D_q and f(alpha)
multifractal spectra.

Program_18a.py --- Generating a multifractal image

Program_18b.py --- Counting pixels in a color image.
Program_18c.py --- Image and statistical analysis on
microbes.png
Program_18d.py --- Fast Fourier transform of a noisy signal
Program_18e.py --- Iterative map and power spectra
Program_18f.py --- Fast Fourier transform of Lena image
Program_18g.py --- Edge detection in Lena image

Program_19a.py --- Chaos control in the logistic map.

Program_19b.py --- Chaos control in the Henon map.
Program_19c.py --- Chaos synchronization between two Lorenz
systems.
Program_19d.py --- Generalized synchronization.

Program_20a.py --- The generalized delta learning rule.

Appendix 649

Program_20b.py --- The discrete Hopfield network.

Program_20c.py --- Iteration of a minimal chaotic neuromodule.
Program_20d.py --- Bifurcation diagram of neuromodule.

Program_21a.py --- The Hodgkin-Huxley equations.

Program_21b.py --- The Fitzhugh-Nagumo half-adder.
Program_21c.py --- Phase portrait Josephson junction limit
cycle.
Program_21d.py --- Animated Josephson junction limit cycle.
Program_21e.py --- Pinched hysteresis of a memristor.
Index

absorptive nonlinearity, 407 ArtistAnimation, 178

action potential, 114, 559 assay for neuronal degradation,
activation 577
function, 521, 550 associative memory, 524, 531
level, 532 asymptotic expansion, 127
potential, 522 asymptotically stable
ADALINE network, 524 critical point, 152
affine linear transformation, 440 asynchronous updating, 536
age class, 106, 333 attractor, 193
Airy equation, 61 attributes, 529
algebraicity of limit cycles, 283 autocatalysis, 204
All or None principle, 560 autonomous differential equation,
Alzheimer’s disease, 577 54
ampere, 48 autonomous system, 186
Ampere’s law, 405 auxiliary system approach, 507
Anaconda, 14 average Lyapunov exponent, 366
Anaconda Python programs, 646 ax.set title, 235
anemia, 374 Axes3D(fig), 235
angiogenesis, 308 axial flow compressors, 175
angular frequency of the wave, axon, 521, 558
406
animation, 26 Baby computer, 571
Spyder, 646 backpropagation, 526
ants and termites, 96 algorithm, 528
aperiodic, 201, 360 backward training, 388
behavior, 196 ballistic propagation, 571
append, 262 bandwidth, 407
applying a damper, 493 Barnsley’s fern, 441
arrhythmic, 495 basin of attraction, 159, 193, 372,
Artificial Intelligence Group, 526 388
artificial neural networks, 520 basis, 253

© Springer International Publishing AG, part of Springer Nature 2018 651

S. Lynch, Dynamical Systems with Applications using Python,
https://doi.org/10.1007/978-3-319-78145-7
652 Index

batch data, 530 blowflies, 360

Belousov-Zhabotinski reaction, bluegill sunfish, 104
204, 212 Boston housing data, 529
Bendixson’s criterion, 123 boundaries of periodic orbits, 391
bias, 521 box-counting dimension, 443, 448
bifurcating limit cycles from a brain functions, 521
center, 258 Briggs-Rauscher, 204
bifurcation bursting, 310
curve, 167 butterfly effect, 200
diagram, 166 BZ reaction, 204
at infinity, 273
point, 362
canny edge detector, 477
value, 164
canonical form, 66, 67, 186
bifurcation diagram
Cantor
CR resonator, 409
multifractal, 454
DDE
field components, 313 set, 235, 435
Mackey-Glass, 307 capacitance, 49
Duffing equation, 232 capacitor, 49
Gaussian map, 370 cardioid, 393
Ikeda map, 425 cardiology, 492
Josephson junction, 575 carrying capacity, 38
limit cycles, 177 cavity ring (CR) resonator, 408
logistic map, 365 cavity round-trip time, 408
neuromodule, 544 cell body, 558
periodically forced center, 69, 148, 246
pendulum, 233 manifold
SFR resonator, 422 theorem, 191
SNIC, 178 changing the system parameters,
binarize, 475 493
binary half adder, 565 chaologist, 492
biology, 374 chaos, 194, 196, 350, 353, 409
bipolar activation function, 522 control, 541, 552
bistability, 174, 175, 368, 421 synchronization, 505
bistable, 175, 206, 370, 408 chaos control
cycle, 203 OGY method, 493
device, 407, 421 periodic proportional pulses,
neuromodule, 543 513
optical resonator, 495 chaos game, 439
region, 233, 241, 409, 414 chaotic
solution, 176 attractor, 197, 229, 415
bistable region, 544 dynamics, 199
blit, 178 phenomena, 348
Index 653

chaotic attractor completely

Hénon map, 501 integrable, 221
neuromodule, 542 reduced, 254
Sierpiński, 440 complex eigenvalues, 69, 330
Chapman cycle, 211 complex iterative equation, 417
characteristic complex(x,y), 396
equation, 302, 330 compound interest, 329
exponent, 283 computer algebra, 252
multiplier, 219 concentrations, 44
charge density, 405 conditional Lyapunov exponents,
chemical 506
kinetics, 44, 86, 117, 204, conductivity, 406
505, 561 conformal mapping, 386
reaction, 59 conservation of energy, 146
signals, 521 conservation of mass, 86
substance, 61 conservative, 146
chemical law of mass action, 44 contact rate, 212
Chua’s circuit, 51, 117, 201, 493, content-addressable memory, 532
514 continuation lines, 20
circle map, 221 continuous Hopfield model, 532
circular frequency of light, 410 contour plot, 24
classical symmetry argument, 248 control curves, 499
classification of critical points, 71 control engineering, 526
climate change, 314 control parameter, 495
clipping problem, 451 control region, 495
clockwise bistable cycle, 632 controlling chaos
clockwise hysteresis, 241, 410 Hénon map, 498
cluster, 524 logistic map, 497
cmap, 396 conversational agents, 526
CMOS oscillators, 571 convex closed curve, 123
coarse Hölder exponent, 450 convoluted surfaces, 191
codimension-1 bifurcation, 182 coordinates of an image, 474
codimension-2 bifurcation, 182 core area of the fiber, 413
coexistence, 96 corollary to Poincaré-Bendixson
coexisting chaotic attractors, 372 theorem, 120
col, 69 correlation dimension, 450
collect, 128 coulomb, 48
comb(k,s), 459 counterclockwise hysteresis, 410
common typing errors, 20 coupler, 412
commutative ring, 252 critical point, 55, 72, 186, 191,
competing species, 96, 110, 140 300
complete synchronization, 506 at infinity, 276
654 Index

culling, 106 dimension, 448

policy, 337 direction
current, 48 field, 66
density, 405 vector, 66
cusp, 84 discrete Fourier transform, 479
cylindrical polar coordinates, 193 discrete Hopfield model, 536
dispersive nonlinearity, 407
damping, 116 displacement function, 247
damping coefficient, 282 distributive laws, 252
dangerous bifurcation, 174 divergence test, 248
Daphnia dentifera, 104 domain of stability, 97, 193, 388
dashed curve, 22 double-coupler fiber ring
data, 545 resonator, 410, 428
data mining, 520 double Hopf bifurcation, 308, 322
databases, 529 double-scroll attractor, 203
DDE, 298 double-well potential, 151
dde23, 315 Dq , 448
ddeint, 298 driver system, 506, 507
def, 6 dsolve, 57
defibrillator, 495 Duffing
defraction, 61 equation, 129, 227
degenerate system, 493, 611
critical point, 148 Dulac’s criterion, 122
node, 70 Dulac’s theorem, 273
degree, 273
degree lexicographical order, 253 EC , 191
delay differential equation, 298 economic model, 117, 322, 332
deleted neighborhood, 121 economics, 92, 374, 382, 453, 505
delta learning rule, 524, 526 edge detection
dendrites, 521, 558 Roberts, 485
depolarization, 559 Sobel, 485
depolarize, 566 eig, 342
derivative of the Poincaré map eigenvector, 72
test, 220 El Niño, 313
desired vector, 527 electric
deterministic chaos, 195, 492 circuit, 48, 92, 117, 201, 532
deterministic system, 520 displacement, 405
Df , 441 displacement vector, 405
dielectric, 406 field, 412, 417, 495
difference equation, 328, 551 field strength, 405
differential amplifier, 407 flux density, 405
diffusion limited aggregates electromotive force (EMF), 49
(DLA), 453 elementary steps, 44
Index 655

elliptic integral, 259 Fibonacci sequence, 343

EMF, 49 field, 253
energy level, 227 figsize, 107
enrichment of prey, 106 fine focus, 246
ENSO model, 313 first integral, 146
Enthought Canopy, 14 first iterative method, 419, 422,
environmental effects, 106 552
environmental model, 313 first return map, 216
epidemic, 61, 85, 106, 117, 212 first-order difference equation,
epilepsy, 577 328
epoch, 524 fish population, 38, 181, 345
equilibrium point, 55 Fitzhugh-Nagumo
ergodicity, 366, 495 equations, 116
error backpropagation rule, 528 oscillator, 114
error function, 527 system, 563
erythrocytes, 374 fixed point, 55, 355
ES , 72, 186, 191 period m, 221
EU , 72, 186, 191 period N, 357
Euclidean dimension, 448 period one, 217, 414, 497
Euclid’s algorithm, 29 period two, 498
exact, 39 fixed size box-counting algorithm,
exact differential equation, 39 451
excitatory, 521, 558 fixed weight box-counting
existence and uniqueness algorithm, 451
limit cycle, 117 flow, 118
existence theorem, 53 focal values, 247
extinct, 96 fold bifurcation, 182
for loop, 6
Fabry-Perot forced system, 227
interferometer, 407 forward rate constant, 46
resonator, 407 fossil dating, 59
farad, 49 Fourier spectrum, 478
Faraday’s law, 49 Fourier transform, 477
of induction, 404 fractal, 434, 441
fast Fourier transform, 480 attractor, 197, 441
feedback, 175, 407, 420 dimension, 441
feedback mechanism, 543 Cantor set, 442
feedforward single layer network, Koch curve, 442
523 Koch square, 442
Feigenbaum constant, 365 Sierpiński triangle, 442
Ferranti Mark 1, 571 geometry, 434
FFT, 480 structure, 196, 201, 386
fiber parameters, 423 fragmentation ratios, 449
656 Index

f (α) spectrum, 448 Hamiltonian systems

FuncAnimation, 178, 207 with two degrees of freedom,
Function, 57 221
function approximators, 526 handcrafted patterns, 541
fundamental memory, 536 hard bifurcation, 174
fuzzy discs, 451 Hartman’s theorem, 79
harvesting, 106, 181
Gauss’s law policy, 337
electricity, 405 Hausdorff dimension, 448
magnetism, 405 Hausdorff index, 441
Gauss-Newton method, 528 Hausdorff-Besicovich dimension,
Gaussian input pulse, 421 450
Gaussian map, 368 Heaviside function, 523
Gaussian pulse, 638 Hebb’s learning law, 522
generalized delta rule, 528 Hebb’s postulate of learning, 536
generalized fractal dimensions, help command, 20
448 Hénon map, 370, 445, 451, 609,
generalized mixed Rayleigh 612
Liénard equations, 267 henry, 49
generalized synchronization, 507 heteroclinic
gestation period, 304 bifurcation, 234
GitHub, 2 orbit, 120, 150, 234, 274
global bifurcation, 260, 273 tangle, 234
global warming, 314 heterogeneous, 448
globally asymptotically stable, hidden layer, 524, 528, 603
194, 211 high pass filter, 484
glucose in blood, 60 Hilbert numbers, 273
Gröbner bases, 252 Hints for programming, 20
gradient, 66 history, 176
gradient vector, 527 history function, DDEs, 298
graphene nano-ribbon, 573 Hodgkin-Huxley equations, 114,
graphic, 120 560
graphical method, 351, 417 Holling-Tanner model, 101, 139,
gray scale, 458, 473 164
Green’s theorem, 122 homoclinic
Gross National Product (GNP), bifurcation, 200, 234, 261
374 loop, 260, 261, 267
Guido van Rossum, 2 orbit, 150, 234, 274
tangle, 234
Hénon-Heiles Hamiltonian, 225 homogeneous, 448
haematopoiesis, 306 homogeneous differential
Hamiltonian, 146, 611 equation, 40
Index 657

Hopf information dimension, 450

bifurcation, 170, 182, 303 inhibitory, 521, 558
singularity, 182 initial value problem, 38
Hopfield network, 156, 531, 551, input vector, 521
610, 613 insect population, 109, 345, 616
Hopfield neural network, 525 instability, 421
horseshoe dynamics, 235 instant physician, 526
host-parasite system, 104 integrable, 223
human population, 85, 343 integrate and fire neuron, 116
hyperbolic Integrated Development
attracting, 283 Environment, 2, 14
critical point, 79 integrating factor, 34
fixed point, 220, 371 intensity, 412
iterated function system, 440 interacting species, 95, 626
repelling, 283 intermittency, 203, 212, 363
stable limit cycle, 220 route to chaos, 365
unstable limit cycle, 220 invariant, 118, 201, 416
hyperpolarize, 566 axes, 82, 98
hyperpolarized, 559 inverse discrete Fourier
hysteresis, 175, 371, 421 transform, 479
curves, 313 inverted Koch snowflake, 609
Josephen junction, 574 inverted Koch square, 438
inverted pendulum, 322
IDE, 14 io.imsave, 474
ideal, 252 isoclines, 67
IDLE, 2 isolated periodic solution, 114
IDLE Python programs, 645
isothermal chemical reaction, 86
if, elif, else, 6
iterated function system (IFS),
Ikeda
440
DDE, 310
iteration, 328
map, 375, 414, 428, 513
im, 396
image analysis, 453 Jacobian, 171
image compression, 434, 484 Jacobian matrix, 78, 192, 205,
incident, 407 371, 502
indentation level, 20 Jaynes-Cummings model, 595
index, 125 Jordan curve, 121, 286
inductance, 49 Josephson junction
infected population, 105 mathematical model, 573
infectives, 85 Josephson junction (JJ), 571
inflation unemployment model, jth point of period i, 357
382 Julia set, 386, 389, 434, 612
infodict, 107 color, 390
658 Index

KAM system, 116, 123, 139, 260,

theorem, 227 281
tori, 227 large parameter, 286
kernel machines, 525 local results, 290
Kerr theorem, 293
effect, 407, 413 limit cycle, 106, 114, 118, 206,
type, 413 595, 611
kinetic energy, 146 hyperbolic, 259
Kirchhoff’s neuron, 116
current law, 49 nonexistence, 608
laws, 532 3-D, 195
voltage law, 49 Lindstedt-Poincaré technique, 131
Koch linear differential equation, 34
curve, 436 linear phase shift, 412, 423
snowflake, 464 linear stability analysis, 56, 300,
square, 436 417
linear transformation, 187
linearization, 78
ladybirds and aphids, 99 linearized system, 78
lambdas, 315 Lipschitz
laminarize, 494 condition, 53
Lang-Kobayashi equations, 310 continuous, 53
Laplace transform, 50 local bifurcation, 273
large-amplitude limit cycle, 175 log-log plot, 445
bifurcation, 175 logic gates, 407
laser, 182, 374, 410, 494 logic operations, 565
model, 310 logistic
LaTeX, 24 equation, 38, 302
law of mass action, 86 function, 6
learning process, 521 growth, 101
learning rate, 527 map, 360, 497, 612
least mean squared (LMS) Lorenz
algorithm, 524 attractor, 201
legend, 23, 107 equations, 199, 494
Legendre transformation, 450 loss in the fiber, 412
Leslie Lotka-Volterra model, 99, 164,
matrix, 334 212, 304
model, 333 low-gain saturation function, 523
lexicographical order, 253 low pass filter, 484
lie detector, 526 lowest common multiple, 256
Liénard Lyapunov
equation, 248 quantity, 247
plane, 282 stability, 531
Index 659

Lyapunov domain of stability, 155 mathematical model, 574

Lyapunov exponent, 197, 366, 612 meshgrid, 157
Lorenz system, 603 meteorology, 199
Lyapunov function, 151, 154, 194, method of multiple scales, 134
247, 284, 550, 607 method of steepest descent, 527
Hopfield network, 532 method of steps, 299
Lyapunov quantities, 290 mgrid, 425
Lyapunov stability theorem, 152 micro-parasite—zooplankton—
lynx and snowshoe hares, 99 fish system,
104
Mac OS, 2 minimal chaotic neuromodule,
Mackey-Glass model, 306 542
magnetic field vector, 405 minimal Gröbner basis, 257
magnetic flux, 405 mixed fundamental memories, 537
magnetostrictive ribbon, 494 mixing, 353
Mandelbrot, 443 modulo, 254
Mandelbrot set, 389, 391, 434 monomial, 253
manifold, 72 ordering, 253
Maple, 26 mortgage assessment, 526
math module, 4 motif, 434
Mathematica, 26 mplot3d, 207
MATLAB, 26 multi-electrode array, 577
MATLAB code to Python, 603 multidegree, 254
matplotlib, 18 multifractal, 447, 472, 602
maximal interval of existence, 54, formalism, 448
62, 118 Hénon map, 459
Maxwell’s equations, 404 Sierpiński triangle, 459
Maxwell-Bloch equations, 406 spectra, 448
Maxwell-Debye equations, 406 multistability, 174
McCulloch-Pitts neuron, 522 multistable, 151, 175, 206, 241,
MEA, 577 543
mean, 545 murder, 60
infectivity period, 212 muscle model, 60
latency period, 212 mutual exclusion, 96
mechanical oscillator, 91 myimages, 292
DDE, 321 mylein sheath, 558
mechanical system, 117, 176
Melnikov national income, 332
function, 259 negative
integral, 258 limit set, 118
memory devices, 407 semiorbit, 118
memristance, 52 negatively, invariant, 118
memristor, 51, 572, 573 net reproduction rate, 339
660 Index

network architecture, 521 nullclines, 67

neural network, 375, 505, 521 numerical solutions, 42
DDE, 315 numpy, 18
neuristor, 572
neurodynamics, 541 occasional proportional feedback
neurological assay, 577 (OPF), 494
neuromodule, 541 ODE, 34
neuron odeint, 42, 57
module, 375 OGY method, 495
neuron(s), 114, 521, 551, 558 ohm, 49
neuronal model, 521 Ohm’s law, 48
neurotransmitters, 558 optical
Newton fractal, 395, 602 bistability, 407
Newton’s law of cooling, 60 computer, 407
fiber, 410
Newton’s law of motion, 146
fiber double ring, 410
Newton’s method, 395, 528
memories, 407
noise, 497
oscillators, 572
NOLM, 404
resonator, 176
with feedback, 410
sensor, 408
nonautonomous system, 116, 227
optimal sustainable, 340
nonconvex closed curve, 124
optogenetics, 578
nondegenerate
orbit, 66, 118
critical point, 148, 246
ordinary differential equation, 34
nondeterministic chaos, 195, 492
oscillation of a violin string, 114
nondeterministic system, 520
oscillatory threshold logic, 563
nonexistence of limit cycles, 123
output vector, 521
nonhyperbolic
ozone production, 211
critical point, 79, 151, 607
fixed point, 371 parasitic infection, 111
nonlinear Parkinson’s disease, 577
center, 247 partial differential equations, 34
optics, 374 partition function, 448
phase shift, 412 Pascal’s triangle, 464
refractive index coefficient, passive circuit, 51
413 Peixoto’s theorem in the plane,
nonlinearity, 175, 407 164
nonperiodic behavior, 196 pendulum, 147, 159, 241
nonsimple canonical system, 67 double, 593
normal form, 164, 170 perceptron, 522
normalized eigenvector, 340 perihelion, 593
not robust, 101 period, 259
notebook, 14 bubblings, 368
np.mgrid, 88 limit cycle, 103, 119, 205
Index 661

undoublings, 368 positive

period-doubling, 203 limit set, 118
period-doubling bifurcations to semiorbit, 118
chaos, 363 positively, invariant, 118
period-n cycle, 195 potato man, 394
period-one behavior, 350 potential difference, 48
period-two, 196 potential energy, 146, 151
behavior, 350 potential function, 151
period-three behavior, 351 pow(x,y), 315
periodic power, 412
behavior, 115 law, 443
orbit, 259 spectra, 199, 203
windows, 363 of a waterwheel, 92
periodicity, 348, 353 power-splitting ratio, 412
permittivity of free space, 405 pprint, 57
perturbation methods, 127 Prandtl number, 200
phase portrait, 66 preallocate, 20
phase shift, 412 predation, 104
physiology, 505 rate, 101
piecewise, 300, 315 predator-prey, 117
piecewise linear function, 523 DDE, 304
pinched hysteresis, 52, 574 models, 99
pitchfork bifurcation, 167 system, 109
pixels, 451 probe vector, 536
planar manifold, 187 propagation, 412
plastics, 453 psychological profiling, 526
plt.axes, 292 PyDDE, 298
Poincaré pydelay, 298, 311
compactification, 275 Pyragas’s method, 493
map, 119, 199, 216, 259, 371, Python
495 based exam, 607, 613
section, 216, 611 files download, 2
Poincaré-Bendixson theorem, 120,
227, 281, 287 qth moment, 448
Poisson brackets, 223 qualitative behavior, 66
polar coordinates, 69, 276 qualitatively equivalent, 74
pole placement technique, 496 quasi-periodicity, 221
pollution, 106 quasi-polynomials, 302
polymer, 453 quasiperiodic, 544, 552
population, 92 route to chaos, 203
of rabbits, 86 quasiperiodic forcing, 228
population model, 614 quiver, 88
662 Index

Rössler reversed fundamental memories,

attractor, 194 537
system, 194 RGB image, 474
radioactive decay, 610 rgb2gray, 477
randint(a,b), 459 ring, 252
random behavior, 195 ringing, 241
Raspberry Pi, 2 RLC circuit, 51, 117
rate constant, 45 roach:fish population, 338
rate-determining step, 45 robust, 103
Rational(1,2), 376 Rotating Wave Approximation
rationally independent, 221 (RWA), 312
ravel, 157 rsolve, 342
Rayleigh number, 200 rubbers, 453
Rayleigh system, 115
re, 396 S-polynomial, 256
reaction rate equation, 45 saddle point, 69, 148
real distinct eigenvalues, 68 saddle-node bifurcation, 165
recurrence relation, 328 saddle-node on an invariant cycle
recurrent neural network, 524, bifurcation, 176
531 safe bifurcation, 174
red blood cells, 374 save image, 472
red and grey squirrels, 96 savefig, 22
reduced, 262 scaling, 443, 448
reduced Gröbner basis, 257 scatter, 425
reflected, 407 sea lions and penguins, 96
refractive index, 407 seasonal effects, 106
refractive nonlinearity, 407 seasonality, 212
refuge, 106 second iterative method, 420,
regionprops, 477 422, 552
regulator poles, 496 second order linear difference
relative permeabilities, 406 equation, 330
relative permittivities, 406 second part of Hilbert’s sixteenth
repeated real eigenvalues, 70 problem, 272
repolarization, 559 second-order differential equation,
resistance, 49 50
resonance terms, 173 secular term, 131
resonant, 173 sedimentary rocks, 453
response system, 506, 508 self-similar, 448
restoring coefficient, 282 self-similar fractal, 441
restoring force, 116 self-similarity, 434
restrictions in programming, 291 semistable
return map, 247, 607 critical point, 56
reverse rate constant, 46 limit cycle, 118, 220, 286
Index 663

sensitivity to initial conditions, manifold, 72, 79, 186, 191,

196, 348, 353, 492 495
separable differential equation, 35 node, 69
separation of variables, 35 staircases, 356
separatrix, 151 stationary point, 55
cycle, 261 std, 545
series solutions, 42 steady state, 51, 103
SFR, 403 stem cell, 578
resonator, 409, 412 stiff system, 47, 212
sharks and fish, 99 stiffness, 116
Sierpiński triangle, 439 stochastic methods, 520
sigmoid function, 523 stock market analysis, 453
signal processing, 453, 482 stoichiometric equations, 45
simple canonical system, 68 Stokes’s theorem, 405
simple nonlinear pendulum, 146 strange attractor, 197
simply connected domain, 123 stretching and folding, 348
singlet, 212 strictly dominant, 336
singular node, 70
structurally
Smale horseshoe map, 234, 373
stable, 103, 164
Smale-Birkhoff theorem, 235
unstable, 101, 164
small perturbation, 56, 417
subcritical Hopf bifurcation, 174,
small-amplitude limit cycle, 246
200
soft bifurcation, 174
subharmonic oscillations, 229
solar system, 492
subplots, 22, 107
solution curves, 36
summing junction, 521
solve, 88, 258
soma, 521, 558 superconductor, 573
spatial vector, 405 supercritical Hopf bifurcation,
spectrum of Lyapunov exponents, 174
197 supervised learning, 524
speed of light, 406 surface plot, 24
spike train, 559 susceptible population, 105
spin-glass states, 537 susceptibles, 85
spirals, 355 sustainable, 338
spurious steady state, 537 switches, 407
SR flip-flop, 569 symbols, 57
stability, 151, 190 sympy, 15
diagram, 420 synaptic
stable cleft, 558
critical point, 55, 152 gap, 558
fixed point, 361, 388 vesicles, 558
focus, 69 weights, 521
limit cycle, 103, 118 synchronization, 494, 505, 541
664 Index

synchronization of chaos, 505 uint8, 476

synchronous updating, 537 unconstrained optimization
problem, 527
target vector, 524, 527 uncoupled, 187
targeting, 497 uniform asymptotic expansion,
τ (q), 449 127
Taylor series expansion, 56, 78, uniform harvesting, 341
371, 417 unipolar activation function, 522
tent map, 348, 608 uniqueness theorem, 53
3D plot, 23 universality, 365
three-dimensional system, 186 Unix, 2
threshold, 559 unstable
logic, 564 critical point, 55, 153
value, 85 fixed point, 361, 388
time series, 105, 364 focus, 69
chaos detection, 198 limit cycle, 118
plot, 198 manifold, 72, 79, 186, 191
Tinkerbell map, 601 node, 68
Toda Hamiltonian, 240 unsupervised learning, 524
topological dimension, 448
topologically equivalent, 74 vacuum, 405
torus, 228 vacuum tube oscillator, 571
total degree, 253 value of homes in Boston, 529
totally van der Pol equation, 130
connected, 388 van der Pol system, 114, 259
disconnected, 388 vector field, 66
training, 524 plot, 533
trajectory, 66, 118 velocity of light, 413
transcritical bifurcation, 167 Verhulst’s equation, 38, 96
transfer function, 521, 551 virus, mobile phone, 592
transient, 51 viscosity, 200
transmitted, 407 viscous fingering, 434
transversal, 247 volt, 48
transversely, 216 voltage drop, 48
travelling salesman problem, 531
triangular pulse, 421 wave equations, 404
trigsimp, 130 wave vector, 406
trivial fixed point, 356 wavelength, 406
turbulence, 453, 492 light, 413
2D plot, 21 WC , 191
two-neuron module, 533 while loop, 6, 29
 Index 665

Windows, 2 XOR gate, 522

wing rock, 175
WinPython, 14 You Tube, 386
WS , 79, 191 youngest class harvesting, 339
WU , 79, 191
X-ray spectroscopy, 453 Zq , 448

View publication stats