Pymc inference for a high dimensional ODE system

I suspect the unpacking from the vector? Just pass the variables separately. Otherwise there’s nothing I see that could fail