Simulating System Dynamics Models in Python

made-with-sphinx-doc JOSS Maintained PyPI license conda package PyPI package PyPI status PyPI pyversions Contributions

This project is a simple library for running System Dynamics models in Python, with the purpose of improving integration of Big Data and Machine Learning into the SD workflow.

PySD translates Vensim or XMILE model files into Python modules, and provides methods to modify, simulate, and observe those translated models. The translation is done throught an intermediate Abstract Synatax Tree representation, which makes it possible to add builders in other languages in a simpler way

Why create a new SD simulation engine?

There are a number of great SD programs out there (Vensim, iThink, AnyLogic, Insight Maker, and others). In order not to waste our effort, or fall victim to the Not-Invented-Here fallacy, we should have a very good reason for starting a new project.

That reason is this: There is a whole world of computational tools being developed in the larger data science community. System dynamicists should directly use the tools that other people are building, instead of replicating their functionality in SD specific software. The best way to do this is to bring specific SD functionality to the domain where those other tools are being developed.

This approach allows SD modelers to take advantage of the most recent developments in data science, and focus our efforts on improving the part of the stack that is unique to System Dynamics modeling.


Currently PySD does not implement all the functions and features of Vensim and XMILE. This may mean that some models cannot be run with the Python translated version, or can only be partially run. In most cases, functions that are not implemented will be translated as pysd.py_backend.functions.not_implemented_function(). However, the most used functions and features are implemented in PySD and most of the models will run properly.

For more information, see the sections on supported Vensim functions, supported Xmile functions, and supported Python builder functions. In case you want to add any new functions, please follow the tips in the development section. The examples of adding functions section may help you.

Additional Resources

PySD Cookbook

A cookbook of simple recipes for advanced data analytics using PySD is available at:

The cookbook includes models, sample data, and code in the form of iPython notebooks that demonstrate a variety of data integration and analysis tasks. These models can be executed on your local machine, and modified to suit your particular analysis requirements.



PySD is currently a community-maintained project, any contribution is welcome.

The code for this package is available at:

If you find any bug, or are interested in a particular feature, see reporting bugs.

If you are interested in contributing to the development of PySD, see the developer documentation listed above.


If you use PySD in any published work, consider citing the PySD Paper (2022):


Martin-Martinez et al., (2022). PySD: System Dynamics Modeling in Python. Journal of Open Source Software, 7(78), 4329,
   doi = {10.21105/joss.04329},
   url = {},
   year = {2022},
   publisher = {The Open Journal},
   volume = {7},
   number = {78},
   pages = {4329},
   author = {Eneko Martin-Martinez and Roger Samsó and James Houghton and Jordi Solé},
   title = {PySD: System Dynamics Modeling in Python},
   journal = {Journal of Open Source Software}

Please, also add the PySD Introductory Paper (2015):

Houghton, J. P., & Siegel, M. (2015). Advanced data analytics for system dynamics models using PySD. Proceedings of the 33rd International Conference of the System Dynamics Society, 2, 1436–1462. ISBN: 9781510815056
   author = {Houghton, James P and Siegel, Michael},
   booktitle = {{Proceedings of the 33rd International Conference of the System Dynamics Society}},
   publisher = {{System Dynamics Society}},
   title = {{Advanced data analytics for system dynamics models using PySD}},
   url = {},
   isbn = {9781510815056},
   volume = {2},
   pages = {1436-1462},
   eventdate = {2015-07-19/2015-07-23},
   location = {Cambridge, Massachusetts, USA},
   year = {2015},
   month = {7},
   keywords = {System Dynamics, Vensim, Python}

You can also cite the library using the DOI provided by Zenodo. It is recomendable to specify the used PySD version and its correspondent DOI. If you want to cite all versions you can use the generic DOI for PySD instead:



For additional help or consulting, contact or