DataMuseum.dk

algol<
begin
_____
integer n,i,j;
_______
real my,sigma,NFTT,p,t,alfa1,alfa2,u,v,w;
____
 boolean ran;
 _______
 real random;
 ____
 core code ran,random;
 ____ ____
 2,46 ; boolean ran
        _______
 2,45 ; real random
        ____
 e:  pm re1,mln re2
     gr re1 X
     sr re1,gr re1
     hv re LKA
    tl -6,ga re3
     pm re1,mln re2
     gr re1 X
     sr re1,gr re1
     it re4
 e3: pm _ X
        0
     gm (r-1),nkf
     hv re LZ
     grf a2,hr s1
 e1: bsn (p210) DLPB t -451
 e2: m 9901
     _
 e4: qq
     qq
     qq
     qq
     qq
     qq
     qq
     qq
     e;
     _
select(35);writecr;writecr;
n:=read integer;
begin
_____
real array X[1:2;1:n];
____ _____
for j:=1 step 1 until n do
___      ____   _____   __
X[1,j]:=j;
for j:=1 step 1 until n do
___      ____   _____   __
X[2,j]:=read real;p:=read real;
if p|0 then
__  =  ____
 begin
 _____
 writetext(|<  Wrong number of elements in the datagroup|);
           <                                            >
 writecr;writecr;
 go to end;
 __ __
 end;
 ___
alfa1:=read real;alfa2:=read real;u:=X[2,1];v:=X[2,2];writecr;
my:=read real;sigma:=read real;j:=0;
writetext(|<  TIDSSERIEANALYSE.(Yule-serier)|);writecr;writecr;writecr;
          <                                 >
writetext(|<  Test for randomness|);writecr;writecr;
          <                      >
writetext(|<  alfa1= |);write(|-dd.dddd|,alfa1);writecr;
          <          >        <        >
writetext(|<  alfa2= |);write(|-dd.dddd|,alfa2);writecr;
          <          >        <        >
writetext(|<  sigma= |);write(|_dd.dddd|,sigma);writecr;writecr;
          <          >        <,       >
for i:=1 step 1 until n-2  do
___      ____   _____      __
 begin
 _____
 gier(ran);p:=random;
if p_.5 then
__  <   ____
  begin
  _____
  t:=sqrt(-2⨯ln(p));
  NFTT:=(t-(2.515517+.802853⨯t+.010328⨯t∧2)/
                                        |
        (1+1.432788⨯t+.189269⨯t∧2+.001308⨯t∧3))⨯sigma+my;
                               |           |
  end else
  ___ ____
  begin
  _____
  p:=1-p;
  t:=sqrt(-2⨯ln(p));
  NFTT:=-(t-(2.515517+.802853⨯t+.010328⨯t∧2)/
                                         |
         (1+1.432788⨯t+.189269⨯t∧2+.001308⨯t∧3))⨯sigma+my;
                                |           |
  end;
  ___
 j:=j+1;
 w:=-alfa1⨯u-alfa2⨯v+NFTT;X[2,j]:=w-X[2,j+2];
 write(|____-dddd.ddddd|,X[2,j]);
       <,,,,           >
 if entier(j/5)=j/5 then writecr;
 __                 ____
 u:=w;v:=u;
 end;
 ___
end:
end;123 121 end;t<
___         ___