tutorials on using vector variables programattically

[isaacsas]

And note we currently need to have users call collect(vec_symbolic_variable) before passing to ReactionSystem.