# amari bump
mh(x)=exp(-abs(x))-c*b*exp(-b*abs(x))
par c=1.05,b=.5
table w % 101 -50 50 mh(dx*t)*dx
par dx=.1
fu[0..199]=heav(u[j])
special z=conv(even,200,50,w,fu0)
u[0..199]'=-u[j]-h+z([j])
par h=.15
@ total=100,meth=euler
init u[90..110]=.8
done