\\ 20894 Menezes/, 182.
t_autocorr (v,d) = {my (s=0,n=#v);
for (i=1,n-d,s+=(v[i]!=v[i+d])); s}

t_erf (x) = 1-erfc(x)

\\ Kolmogorov-Smirnov test function, summing for t*k<=up.
\\ Checked vs. 1-Q(\lambda) in 2031 Fisz, 726.
t_ks (t,up=20,maxt=200) = {my (s,sgn,k,t2t2);
s=0; sgn=1; k=1; t2t2=2*t*t;
while ((k<=maxt)&&(t*k<=up),\
s+=(sgn*exp(-t2t2*k*k)); k+=1; sgn=-sgn); 2*s}

\\ Maximum of forward partial |sums|.
t_maxsum (v) = {my (s=0,m=0);
for (i=1,#v, s+=v[i]; sa=abs(s); if (sa>m,m=sa)); m}

\\ N(0,1).
t_Phi (x) = 0.5*(1+t_erf(x/sqrt(2)))