DataMuseum.dk

algol<
begin
_____
integer n,p,r,i,j,k,m,t;
_______
real x,sum,P;
____
select(43);
n:=read integer;m:=read integer;
   begin
   _____
   real array A[0:n];integer array B[1:m];
   ____ _____        _______ _____
   A[0]:=x:=0;
   for i:=1 step 1 until n do
   ___      ____   _____   __
   A[i]:=x:=x+ln(i);
   for j:=1 step 1 until m do
   ___      ____   _____   __
   B[j]:=read integer;
   for i:=1 step 1 until m do
   ___      ____   _____   __
      begin
      _____
      p:=B[i];
      for j:=i+1 step 1 until m do
      ___        ____   _____   __
         begin
         _____
         r:=B[j];t:=read integer;
         if -,(n_500∧p_n∧r_p∧t_r) then
         __     <     <   <   <   ____
            begin
            _____
            writetext(|<fejl i datasæt|);writecr;
                      <               >
            end else
            ___ ____
            begin
            _____
            sum:=0;
            for k:=0 step 1 until r do
            ___      ____   _____   __
                  begin
                  _____
                  if n-p-r+k<0 then P:=0 else
                  __           ____      ____
                  P:=exp(A[p]+A[n-p]+A[n-r]+A[r]-A[n]-A[p-k]-A[k]-A[n-p-r+k]
                        -A[r-k]);
                  if k>t then sum:=sum+P else
                  __     ____            ____
                  if k=t then sum:=x:=P
                  __     ____
                  end;writecr;
                  ___
            writetext(|< D(|);
                      <    >
            write(|dd|,i);
                  <  >
            writetext(|<,|);
                      <  >
            write(|dd|,j);
                  <  >
            writetext(|<).  n=|);
                      <       >
            write(|ddd|,n);
                  <   >
            writetext(|<;p=|);
                      <    >
            write(|ddd|,p);
                  <   >
            writetext(|<;r=|);
                      <    >
            write(|ddd|,r);
                  <   >
            writetext(|<;t=|);
                      <    >
            write(|ddd|,t);writecr;writecr;
                  <   >
            writetext(|<  Null hyp.:No pos.sociability.
                      <
  Level of significance:P(X_t)=|);
                           >   >
            write(|__d.ddddd|,sum);writecr;writecr;
                  <,,       >
            writetext(|<  Null hyp.:No neg.sociability.
                      <
  Level of significance:P(X_t)=|);
                           <   >
            write(|__d.ddddd|,1-sum+x)
                  <,,       >
            end;writecr;writecr;writecr
            ___
         end
         ___
      end
      ___
   end
   ___
end;t<
___