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;
alfa1:=read real;alfa2:=read real;u:=read real;v:=read real;writecr;
my:=read real;sigma:=read real;n:=read integer;j:=0;
writetext(|<  TIDSSERIEANALYSE.(Yule-serier)|);writecr;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 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;
 write(|____-dddd.ddddd|,w);
       <,,,,           >
 if entier(j/5)=j/5 then writecr;
 __                 ____
 u:=w;v:=u;
 end;
 ___
end;t<
___