So I have this ODE to solve (26) with a set of conditions (27). I evaluated the conditions prior to going forward to solve the ODE. I want to solve it symbolically for each mk because I will use the results later on. I wrote a function file for dmk/dt, and then in a loop I want to solve for each mk and save the function. However I am not sure if I am able to save all the functions in a matrix. Moreover, when I input the initial conditions while using ode45, I get an error. This is the code in the loop:
where m0val(x) contains the values for the initial conditions.