p tcp=200
p N=1.0,MaxWt=1.0
p beta1=1.2, beta2=0.2
p Crr=0.9,Cll=0.9,Cinter=0.3,Clr=0.3,Crl=0.3
cint=cinter*heav(tcp-t)
clld=cll-.3*heav(t-tcp)
nl=N*wl/(wl+wr+beta2)
nr=N*wr/(wl+wr+beta2)

sumw=wr+wl

dwr/dt=(Crr*wr+Cint*wl)*(MaxWt-wr)*nr - beta1*sumw*wr
dwl/dt=(Cint*wr+Clld*wl)*(MaxWt-wl)*nl - beta1*sumw*wl



wr(0)=0.03
wl(0)=0.04
@ total=200
@ xp=wl,yp=wr,xlo=0,xhi=1,yhi=1,ylo=0
d





From aehst8+@pitt.edu Fri Jan 15 15:35:23 1999
Date: Fri, 15 Jan 1999 15:29:09 -0500 (EST)
From: Anthony E Harris <aehst8+@pitt.edu>
To: "G. Bard Ermentrout" <bard@math.pitt.edu>
Subject: Re: jcn paper



On Thu, 14 Jan 1999, G. Bard Ermentrout wrote:

> could you also send the actual ODE file - I have over 8000 ode files and
> dont know exactly where to look!
> 
> Bard
> 
> 

bard,

Here is the original odefule. To simulate MD, I just set Cll=0.6,
Cinter=0. here it is: