I have been trying to solve this problem for sometime and so far I have received an answer that requires the symbolic toolbox which I dont have - I would appreciate any help!
I need to solve these 2 differential equations simultaneously.
But I dont know how to code the dr*^3/dt. My code is below:
y(1) = r* and
y(2) = C*
dydt(1) = dr*/dt and
dydt(2) = dC*/dt
In my case dydt(1) needs to be replaced with something that would solve dr*^3/dt and not dr*/dt. The rest of the code is below.
[t,y]=ode45(@(t,y) odefcnNY_v3(t,y,D,Cs,rho,r0,Af,N,V), tspan, y0);