DataMuseum.dk

algol<
begin
_____
integer n,i;real m,s,p;
_______     ____
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);writecr;
n:=read integer;m:=read real;s:=read real;
 begin
 _____
 real array X[1:n];
 ____ _____
 for i:=1 step 1 until n do
 ___      ____   _____   __
 X[i]:=read real;
 for i:=1 step 1 until n do
 ___      ____   _____   __
  begin
  _____
  distrbfunc1(m,s,X[i],p);
  write(|d.dddd|,p);writecr;
        <      >
  end
  ___
 end
 ___
end;t<
___