load(foo.ode) loads the reaction system but then infers reactions, which might not be wanted
tutorial notebooks and the documentation for load/1 say that
load_ode on it (but that command does not exist).
Moreover, what is done is basically
import_ode(foo.ode), load_reactions_from_ode_system. when one might only want to do the first part…
In other words
load(foo.ode) should probably just call
load/1 doc and tutorials might need to be updated