{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 "Text Output" -1 2 1 
{CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 }1 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 
1 0 0 255 1 0 0 0 0 0 0 1 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 }{PSTYLE "" 11 12 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 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 37 "p[-1]:=0: p[0]:=1: q[-1]:=1: q[0]:=0:" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "for i from 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 90 "N:=100000: Digits:=30: alpha:=evalf(p[14]/q[14]); 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 35 "C:=[seq(N*(B[i
+1]-B[i]),i=1..N-1)]:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "E:=sort(C)
:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#>%&alphaG$\"?S$z%QQbU8E52_)*f<!#H" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 111 "disc:=proc() tempvar:=0: for j from 1 to N-1 do temp
var:=max(tempvar,abs(j/N-(1-exp(-E[j])))) od: tempvar; end;" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 7 "" 1 "" {TEXT -1 47 "Warning, \+
`tempvar` is implicitly declared local" }}{PARA 7 "" 1 "" {TEXT -1 41 
"Warning, `j` is implicitly declared local" }}{PARA 12 "" 1 "" 
{XPPMATH 20 "6#>%%discGR6\"6$%(tempvarG%\"jGF&F&C%>8$\"\"!?(8%\"\"\"F0
,&%\"NGF0!\"\"F0%%trueG>F,-%$maxG6$F,-%$absG6#,(*&F/\"\"\"F2!\"\"F0F3F
0-%$expG6#,$&%\"EG6#F/F3F0F,F&F&F&" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 7 "disc();" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\">>fz&Gc'
4%*\\m>GU6#!#I" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}
{MARK "0 0 0" 242 }{VIEWOPTS 1 0 0 1 1 1803 }