szivacs.mws
>
deq1:=diff(m(t),t)=-alpha*m(t)*(kappa-(kappa0+(m0-m(t))/M));
>
with(DEtools):
>
alpha:=1; M:=40; kappa:=0.5;kappa0:=0.05; m0:=0.1;
>
DEplot(deq1,m(t),t=0..10,
[[m(0)=m0]],scene=[t,m(t)],linecolor=t,stepsize=0.1);
>