In Matlab you can design a set of loops where the parameters and the initial values can vary at a specific range. Then you just store the results or do an animation plot of the solutions. That's all in general, in practise it needs work.
Carter uses Mathematica (using Mathematica 7) for his solutions. A very helpful part of Carter's notes is that he has provided Mathematics notebooks with detailed code with comments. These notebook can be ported to a local copy of Mathematica for experiments. Another good place for a discussion about simulating HODEs is given in