These tutorial pages illustrate how to use PyDSTool to study a variety of dynamical systems. It is partly inspired by the XPP tutorial, from which it borrows some of the examples. You are presumed to have read this section on the Getting Started page to get acquainted with some basic steps.
If you have the latest version of the PyDSTool package installed (currently only from the Bazaar repository, see the CodeTopics page), you can follow the tutorial code examples in the PyDSTool/tests/ directory. Unless otherwise specified, all the examples will use integrators that will work "out of the box", but you can adapt the scripts to use the faster C-based integrators if you have the appropriate dependencies installed (see the Getting Started page).
[Examples without links are still in development - come back soon!]
These introductory examples do not necessarily show off the full power of the package, but introduce different aspects in simple steps.
- A long introductory discussion of programmatic elements and Python syntax issues involved in building a very simple linear differential equation model for simple harmonic motion. This is aimed at those new to Python and programmatic environments in general.
- A 1D nonlinear ODE, a Calcium channel model with bifurcation analysis and plotting
- The Lorenz discrete 1D mapping and its Lyapunov exponent
- Using a symbolic Jacobian to speed up integration of a very stiff system
- Adding a zero-crossing event to the Van der Pol oscillator (includes phase plane analysis)
- Importing systems biology models from SBML
- Importing neuroscience models from NineML
- Adding noise to an ODE system
- [More to come...]
- Combining tools: A 4D predator-prey model specified with symbolic objects, using bifurcation analysis and 3D plotting
- Specifiying large systems with macros
- Model reduction and dominant scale analysis: for now, see tests/HH_DSSRT.py and the supplemental material in my PLoS Computational Biology 2012 paper about this.
- Finding periodic orbits using mappings in a reduced 2D biomechanics model
- Parameter estimation example using a Hodgkin-Huxley neuron model
- Hybrid dynamical system models (currently only in PDF format, as the supplemental text to the recent PLoS Compututational Biology article)
- Noisy neuron membranes using symbolic objects and the neural toolbox
Other examples provided in the download:
See the PyDSTool/tests/ directory for all the scripts. Most of these files perform user-understandable tests and often demonstrate features verbosely.
PyDSTool source code is hosted by: