odesolver plot "intermediate" values

Hi!

I am trying to save "intermediate" values that are calculated with an odesolver. I have also asked the question on StackOverflow , but with no response as of yet.

Here is my m file as it is currently.

You can see one of my old attempts to plotting one of those values, effFactor, in lines 1, 34 and 261. I tried returning the value alongside the differentials, but that did not work. The plot of effFactor was just a positive-sloped line.

I did some reading and it seems that outputFcn is the way to go. However, I have never done this and am quite new to Matlab. In the m file, you can also see my attempt to get started with the outputFcn method (lines 249-259), but that did not work out so well.

Please, could someone give a hand?