PROGRAM d13r12(input,output);
(* driver for routine PROBKS *)
CONST
   npts=20;
   eps=0.1;
   iscal=40;
VAR
   alam,aval : real;
   i,j : integer;
   txt : PACKED ARRAY [1..iscal] OF char;

(*$I MODFILE.PAS *)
(*$I PROBKS.PAS *)

BEGIN
   writeln('probability function for kolmogorov-smirnov statistic');
   writeln;
   writeln('lambda':7,'value:':10,'graph:':13);
   FOR i := 1 to npts DO BEGIN
      alam := i*eps;
      aval := probks(alam);
      FOR j := 1 to iscal DO BEGIN
         IF (j <= round(iscal*aval)) THEN BEGIN
            txt[j] := '*'
         END ELSE BEGIN
            txt[j] := ' '
         END
      END;
      writeln(alam:8:6,aval:10:6,' ':5,txt)
   END
END.