`a_vec=[0.020 0.025 0.030];for i=1:length(a_vec)dP_dt_1(i)=@(P) (a_vec(i).*P)-(b.*P.^2);[t_ode45_vec_1,y_ode45_vec_1] = ode45(@(P,t) dP_dt_1(t),[t_0 t_max],P_0);figure(5)plot(t_ode45_vec_1,y_ode45_vec_1)hold onend`

Hi everyone, when I try to run this piece of code, I get an error :

`Nonscalar arrays of function handles are not allowed; use cell arrays instead.`

How can I resolve this error?

Thanks in advance!

## Best Answer