I am pleased that you found my code helpful!
The code fits only one vector, so it only has one output. To get all the results, the easiest way would be to copy the differential equation function, and in the second version of it (with a slightly different name), set:
S = Sv;
and then run it with the fitted parameters. Use a different ‘tspan’ as well, with a higher resolution to get a smoother curve (use linspace to create it), if desired. Another option would be to save it and the fitted parameters to a .mat file.
.
