{VERSION 3 0 "IBM INTEL LINUX" "3.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 
0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }
{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 
0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 
0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 242 "This worksheet computes a
 \"bad\" number, i.e. one where the M-level correlations are not all w
hat they should be, as constructed in Sarnak's paper.  The purpose of \+
this particular worksheet is only to spit out data to be analyzed with
 \"cdf.c\"." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "Digits:=30:v
[-1]:=0: r[-1]:=0: r[0]:=1: v[0]:=1: r[1]:=1: v[1]:=1:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "for i from 2 to 16 do v[i]:=r[i-1]*
v[i-1]^2+v[i-2]: r[i]:=floor(log[2](v[i])) od:" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "v[3];" }
}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"&" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 74 "a[0]:=1: for i from 1 to 14 do a[i]:=r[i-1]^2*v[i-1]^
2+2*r[i-1]*v[i-2] od:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "a[0]+1/(a[1]+1/(a[2]+1/a[3])
);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"#W\"#D" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 6 "44/25;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"#W
\"#D" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "4p[-1]:=0: p[0]:=1:
 q[-1]:=1: q[0]:=0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "for \+
i 6from 1 to 14 do p[i]:=a[i-1]*p[i-1]+p[i-2]: q[i]:=a[i-1]*q[i-1]+q[i
-2]: od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 112 "N:=50000: Digits:=30: alpha:=evalf
(p[14]/q[14]): alpha:=evalf(sqrt(2)); A:=[seq(evalf(frac(i^2*alpha)),i
=1..N)]:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "B:=sort(A):" }}{PARA 0 
"> " 0 "" {MPLTEXT 1 0 55 "C:=[seq([N*(B[i+1]-B[i]),N*(B[i+2]-B[i+1])]
,i=1..N-2)]:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "save C, \"C.txt\";
" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 
1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&alphaG$\"?@C()o,)[]4tBc8
UT\"!#H" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}}}{MARK "10 4 0" 0 }{VIEWOPTS 1 0 0 1 1 1803 }