Added a method to generate a numerical rhs function for the Park model.