# mt.ode
# the morse-thule sequence of 0's and 1's 
# the sum mod 2 of the 1's in the binary expansion
# of the integers
# look at the power spectrum for z
# 
init n=0
s0=0
x0=n
%[1..16]
s[j]=s[j-1]+mod(x[j],2)
x[j]=flr(x[j-1]/2)
%
n'=n+1
z'=mod(s16,2)
aux ss=s16
@ meth=discrete,total=10000,maxstor=68000
@ bound=100000
@ xlo=0,xhi=10000,ylo=0,yhi=16,yp=ss
done