szivacs.mws

>    deq1:=diff(m(t),t)=-alpha*m(t)*(kappa-(kappa0+(m0-m(t))/M));

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;

alpha := 1

M := 40

kappa := .5

kappa0 := .5e-1

m0 := .1

>    DEplot(deq1,m(t),t=0..10,
[[m(0)=m0]],scene=[t,m(t)],linecolor=t,stepsize=0.1);

[Maple Plot]

>