@ChrisRackauckas I’ve seen that the issue has been worked on ! Thanks for the support !
In the meantime, I’ve been working on the code without integrating it into a package for the moment. And I have another question that may be related to broadcasting.
I saw that it is possible to define initial conditions with dictionaries associating a key made from a problem variable and an initial condition value such as
prob = ODEProblem(fol, [fol.x => 0.0], (0.0, 10.0), [fol.τ => 3.0])
from Getting Started with ModelingToolkit.jl · ModelingToolkit.jl
Should this work with vector-valued initial conditions ?