This paper shows how Python and Scipy can be used to simulate the time-dependent and steady-state behaviour of reaction networks, and introduces Pysces, a Python modelling toolkit.