ODEs inside a Turing model