DataMuseum.dk

algol<
begininteger n,i,j,k;real m,s,p,r,t,D;
____________         ____
procedure parameter(n,X,MEAN,STDEVIATION);
_________
value n;
_____
integer n;real array X;real MEAN,STDEVIATION;
_______   ____ _____   ____
begin
_____
integer i;real sum,squaresum;
_______   ____
sum:=0;squaresum:=0;
for i:=1 step 1 until n do
___      ____   _____   __
 begin
 _____
 sum:=sum+X[i];squaresum:=squaresum+X[i]∧2;
                                        |
 end;
 ___
MEAN:=sum/n;
STDEVIATION:=sqrt((squaresum-sum∧2/n)/(n-1));
                                |
end procedure parameter;
___
procedure distrbfunc1(m,s,x,PROBABILITY);
_________
value m,s,x;
_____
real m,s,x,PROBABILITY;
____
begin
_____
real A;
____
 real procedure psi(y);
 ____ _________
 value y;real y;
 _____   ____
 begin
 _____
 real a;a:=1/sqrt(2);
 ____
 psi:=1-1/(1+0.07052308⨯a⨯y+0.04228201⨯(a⨯y)∧2+0.00927053⨯(a⨯y)∧3
                                            |                  |
     +0.00015201⨯(a⨯y)∧4+0.00027657⨯(a⨯y)∧5+0.00004306⨯(a⨯y)∧6)∧16;
                      |                  |                  |  |
 end real procedure psi;
 ___
A:=(x-m)/s;
if A<0 then PROBABILITY:=0.5-0.5⨯psi(-A) else
__     ____                              ____
PROBABILITY:=0.5+0.5⨯psi(A);
end procedure distrbfunc1;
___
select(19);
k:=0;
nextdatagroup: for n:=read integer while n|0 do
               ___                 _____  =  __
 begin
 _____
 real array X[1:n];k:=k+1;p:=0;m:=0;s:=0;D:=0;t:=0;
 ____ _____
 writetext(|< sample number:|);write(|ddd|,k);writecr;
           <                >        <   >
 writetext(|< N=|);write(|ddd|,n);writecr;
           <    >        <   >
 for i:=1 step 1 until n do
 ___      ____   _____   __
  begin
  _____
  X[i]:=read real;
  if X[i]=0 then
  __        ____
   begin
   _____
   writetext(|< too few elements in the datagroup.|);
             <                                    >
   writecr;writecr;writecr;go to nextdatagroup;
                           __ __
   end;
   ___
  end;
  ___
 i:=0;
 for r:=read real while r|0 do
 ___              _____  =  __
 i:=i+1;
 if i|0 then
 __  =  ____
  begin
  _____
  write(|ddd|,i);
        <   >
  writetext(|<  elements too many in the datagroup.|);
            <                                      >
  writecr;writecr;writecr;go to nextdatagroup;
                          __ __
  end;
  ___
 for i:=n-1 step -1 until 1 do
 ___        ____    _____   __
 for j:=1 step 1 until i do
 ___      ____   _____   __
 if X[j]>X[j+1] then
 __             ____
  begin
  _____
  r:=X[j];X[j]:=X[j+1];X[j+1]:=r;
  end;comment the elements of X have been ordered;
  ___ _______
 for i:=1 step 1 until n do
 ___      ____   _____   __
  begin
  _____
  write(|ddd.dddd|,X[i]);writecr;
        <        >
  end;
  ___
 parameter(n,X,m,s);
 for i:=1 step 1 until n do
 ___      ____   _____   __
  begin
  _____
  distrbfunc1(m,s,X[i],p);
  if abs(p-(i-1)/n)_abs(p-i/n) then
  __               >           ____
  t:=abs(p-(i-1)/n) else t:=abs(p-i/n);
                    ____
  if D<t then D:=t;
  __     ____
  end;
  ___
writetext(|< observed value of DN:|);
          <                       >
write(|ddd.dddd|,D);writecr;
      <        >
writetext(|< level of significance:|);writecr;
          <                        >
writetext(|<mean value:|);
          <            >
write(|ddd.ddddd|,m);writecr;
      <         >
writetext(|<standard deviation:|);
          <                    >
write(|ddd.ddddd|,s);writecr;writecr;writecr;
      <         >
end
___
end;t<
___