Hi,
The error is because of the following lines,
t = t0:dt:tf; %%from time discretization section
[mlSol(t), rpSol(t), TpSol(t)] = dsolve(odes, conds); %%from initial conditions section
"dsolve" returns the symbolic variables as the output. As the time discretization section is written above the initial conditions section, 't' is no longer a symbolic variable, it is now a vector of double. So, mlSol(t) is trying to index the mlSol symbolic variable with the vector 't' and assign the output value, which is not possible as the values of t are 0,0.01,0.02,...etc. It is an invalid indexing.
It can be solved by, moving the time discretization section next to the initial condtions section, so 't' would be a symbolic variable during the execution of the initial conditions section and it would assign the results to mlSol(t) (i.e to a symbolic variable). Refer to this documentation for more information on the output arguments of the "dsolve" function.
%Initial conditions
cond1 = ml(0) == 1.93*10^-6;
cond2 = rp(0) == 0.00076;
cond3 = Tp(0) == 298.15;
conds = [cond1 cond2 cond3];
[mlSol(t), rpSol(t), TpSol(t)] = dsolve(odes, conds);
% time discretization
dt = 0.01;
t0 = 0;
tf = 100;
t = t0:dt:tf;
tspan=[t0 tf];
fplot(mlSol, 'r', 'linewidth', 2)
title('Droplet mass versus Time')
xlabel('Time,t')
ylabel('Mass, mg')
grid on
fplot(rpSol, 'g', 'linewidth', 2)
title('Droplet size versus time')
xlabel('Time,t')
ylabel('Droplet size, mm')
grid on
fplot(TpSol, 'b', 'linewidth', 2)
title('Temperature versus time')
xlabel('time,t')
ylabel('Temperature, K')
grid on
