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(foo.ode)
call 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 import_ode
, the load/1
doc and tutorials might need to be updated