My guess is that the code now needs to be made into an objective function to be used with a parameter estimation routine.
Also see the documentation on the function to be used to estimate the parameters, since the objective function argument list needs to be compatible with the optimisation function requirements for it. The documentation section on Passing Extra Parameters could also be appropriate.
