function insects
clear all
tr=[0 ??];
initv=??;
[t,y]=ode45(@f, tr, initv);
plot(t,y)
ylabel('Number of Insects')
xlabel('time')
%
function rk=f(t,y)
rk= ??;