View source for Kolmogorov test/Fig1

<center><asy>

srand(2014011);

import stats;

int size = 13;
real [] sample = new real[size+1];
real lambda = 1.3/size;
real width = 2.0;

for (int k=0; k<size; ++k) {
  sample[k] = Gaussrand();
}
sample[size] = 10;

sample = sort(sample);

// for (real x : sample ) {
//   write(x);
// }

real x0 = -10;
int k = 0;
for (real x : sample ) {
  filldraw( box( (x0,k/size-lambda), (x,k/size+lambda) ), rgb(0.8,0.8,0.8) );
  draw( (x0,k/size-lambda)..(x,k/size-lambda), currentpen+1.5 );
  draw( (x0,k/size)..(x,k/size), currentpen+1.5 );
  draw( (x0,k/size+lambda)..(x,k/size+lambda), currentpen+1.5 );
  k += 1;
  x0 = x;
  draw( (x,(k-1)/size-lambda)..(x,k/size+lambda) );
}

clip( box((-width,-0.005),(width,1.005)) );

draw ((-width,0)--(width,0),Arrow);
draw ((0,-0.1)--(0,1.3),Arrow);
draw ((-width,1)--(width,1));

draw ((sample[2],0)..(sample[2],2/size));
draw ((sample[size-1],0)..(sample[size-1],0.48), dashed);
draw ((sample[size-1],0.7)..(sample[size-1],1-1/size), dashed);

label("$x$",(width,0),S);
label("$y$",(0,1.3),W);
label("$0$",(0,0),SW);
label("$1$",(0,1),NE);

label("$X_{(1)}$",(sample[0],0),S);
label("$X_{(2)}$",(sample[1],0),S);
label("$X_{(3)}$",(sample[2],0),S);
label("$X_{(n)}$",(sample[size-1],0),S);

label("$F_n(x)+\lambda_n(\alpha)$",(-1.55,0.35));
draw ((-1.35,0.25)..(-1.2,1/size+lambda));
dot((-1.2,1/size+lambda));

label("$F_n(x)$",(0.4,0.3));
draw ((0.4,0.4)..(0.3,8/size));
dot((0.3,8/size));

label("$F_n(x)-\lambda_n(\alpha)$",(1.5,0.6));
draw ((1.6,0.7)..(1.7,1-lambda));
dot((1.7,1-lambda));

shipout(scale(100,100)*currentpicture);
</asy></center>