Python implementation of DDE solver (one delay) using scipy/numpy