# ca100.ode
# 101 cell CA model
# here is the function for the new value
table f ca5.tab 
# add up with equal weight 5 cells
table wgt % 5 -2 2 1
# The convolution is created
special k=conv(periodic,100,2,wgt,u0)
# here is the update
u[0..100]'=f(k([j]))
# make sure XPP knows it is discrete time
@ total=200,meth=discrete
# integrate u[0..100]=heav(ran(1)-.5)
# view the array of values
#
# make sure you have the table ca5.tab
# 6
# 0
# 5
# 0
# 0
# 1
# 0
# 1
# 1
#
done