# swindale.ode
# creates occular dominance stripes
# make a weight matrix - mexican hat 
table wa % 31 -15 15 exp(-a*t*t)-c*exp(-b*t*t)
par  mu=.5,nu=1.3
par a=.1,b=.025,c=.035
# convolve it
special k=conv(periodic,101,15,wa,n0)
# random initial data
n[0..100](0)=ran(1)
# ODEs
n[0..100]'=n[j]*(1-n[j])*(mu*k([j])-nu)
@ total=50
done
