Boolean Delay Equations (BDEs) can simulate surprisingly complex behaviour, despite their relative simplicity. In addition to steady state dynamics, BDEs can also generate periodic and quasiperiodic oscillations, m:n frequency locking and even chaos. Furthermore, the enumerability of Boolean update functions and their compact parametrisation means that BDEs can be leveraged to generate low-level descriptions of biological networks, from which more detailed formulations (e.g. differential equation models) can be constructed. However, although several studies have demonstrated the utility of BDE modelling in computational biology, a current barrier to the wider adoption of the BDE approach is the absence of freely-available simulation software.In this work, we present BDEtools -an open-source MATLAB package for numerically solving BDE models. After giving a brief introduction to BDE modelling, we describe the package's solver algorithms, together with several utility functions that can be used to provide solver inputs and to process solver outputs. We also demonstrate the functionality of BDEtools by illustrating its application to an established model of a gene regulatory network that controls circadian rhythms.BDEtools makes it straightforward for researchers to quickly build reliable BDE models of biological networks. We hope that its ease of use and free availability will encourage more researchers to explore BDE formulations of their systems of interest. Through the continued use of BDEs by the computational biology community, we will no doubt discover their potential applicability to a broader class of biological networks.