Managed to get my model to run just by messing with some with some parameter values however for this problem that is not really relevant.
When the problem is generated using ODAE the error of MAXITERS appears. However, when my problem is generated using ODE problem, using my above code, the parameter samples for which the solution diverges the output is a matrix of zeros as specified by my code.
So I think the issue is when I generate the problem using ODAE. Do you know of a work round for this?