If I have amodel equation:

>>dqpldt=fc*Rc*Nc0*exp(Kgr*t)+(fi_Ri*(exp(-di*t)*(Ni0*di-Nc0*ai+Kgr*Ni0))/(Kgr+di))+(Nc0*ai*exp(Kgr*t))/(Kgr+di)-Ke*qpl;

With parameters' ranges:

fc=[10^(-4),1]; ...

Rc= [4.5*10^(-11),0.45];

Nc0=[1,10^(10)];

Kgr= [3.46*10^(-3),0.578];

and so on for the remaining parameters fi_Ri; ai; di; Ni0; Ke;

For each parameter, how can I divide the range into 21 subintervals using equally spaced nodes and calculate the time (t) at each node, while the other parameters' values are fixed at the baseline?

For the baseline values, I used the following code to obtain the corresponding time (t):

>>syms qpl(t)

>>qpl(t)= ...........;

>>sol = vpasolve(qpl)

Thanks in advance,