Added comments to clarify code
Added overload for simulate to ignore specific parameters and use ParametricNDSolve