PROGRAM sequential
	PARAMETER(l = 960,lp1 = l + 1)
	INTEGER l,lp1,max,iseed,ip,ibm,i,itime,mag
	LOGICAL n(0:lp1)
	REAL*4 p
	WRITE(6,*)'enter p,max,iseed'
	WRITE(6,*) 'possible inputs: p=0.001,max=900,iseed=123456789'
	read(5,*) p,max,iseed
	ip = (2.0*p - 1.0)*2147483648.0
	ibm = 2*iseed - 1
*	initialize lattice
	do 10 i = 1,l
	   ibm = ibm*16807
	   n(i) = ibm .lt. ip
10	continue
	do 100 itime = 1,max
	   mag = 0
*	   accumulate data
	   do 20 i = 1,l
	      if (n(i)) mag = mag + 1
20	   continue
	   WRITE(5,*) itime,mag
	   if (mag .eq. l .or. mag .eq. 0) stop
	   n(0) = n(l)
	   n(lp1) = n(1)
	   do 30 i = 1,l
*	      execute cellular automata rule
	      n(i) = n(i-1) .or. n(i + 1)
30	   continue
100	continue
	stop
	END