Can I use nlparci with the solve function?

Hello

I solved a set of ODEs using the problem based approach from this example and the solve function:

https://www.mathworks.com/help/optim/ug/fit-ode-problem-based-least-squares.html

My ODEs are derived from chemical reactions like in the example. Now I want confidence intervals for the rate constants. FYI: it looks like the solve function is choosing to use fmincon to solve the problem.

I found someone who used nlparci with the lsqcurvefit function

https://www.mathworks.com/support/search.html/answers/548961-how-to-solve-for-a-kinetic-model-system-with-confidence-intervals.html?fq=asset_type_name:answer%20category:curvefit/fit-postprocessing&page=1

My question is, can I use nlparci with the solve function to get confidence intervals? I don't see how to get the output for the jacobian from fmincon, which I think nlparci needs as an input.

Edit: I added the code and data needed to run the code. Be sure to change the filepath as needed.

Also here is the result, The graph shows the fit being performed on four different data sets. Each one has different rate constants for the same set of reactions. There is only measured data for three of the species.