SUBROUTINE AMEAN(X.N»XBAR) 



1 SUBROUTINE AMEAN(X»N»XBAR> AMEAN 



C AMEAN 



C CALCULATES THE ARITHMETIC MEAN OF A SEQUENCE OF DATA POINTS AMEAN 

 C ' — — 'AMEAN 



2 DIMENSION X(l) AMEAN 



3 XBAR=0. AMEAN 

 <♦ DOI 1 = 1. N AMEAN 



5 1 XBAR=XBAR+X(I> AMEAN 



6 XBAR = XBAR/FLOATCN) AMEAN 



7 RETURN AMEAN 



8 END AMEAN 



AMEAN 



- 



1 





FLOAT 



- 



6 





I 



- 



<♦ 



5 



N 



- 



1AG 



l» 



RETURN 



- 



7 





X 



- 



1AG 



2DI 



XBAR 



- 



1A6 



3= 



7 



8 

 9 

 10 

 11 

 12 

 13 



APPROXIMATION TO INVERSE NORMAL DISTRIBUTION 



6= 



FUNCTION TINORM(ALPHArS) TINORM 



DIMENSION A<3).bl3) TINORM 



DATAUI I). 1 = 1 1 3 )/.0 10328 r. 802853* 2. 515517/ ► IB< I > » I=l»3)/.0010308r TINORM 

 1.189269.1.132788/ TINORM 



-T I NORM 

 TINORM 

 -TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 

 TINORM 



IF(. NOT. (ALPHA. GT.O.. AND. ALPHA. LT.l.)) GO TO 1 



X=ALPHA 



IFIX.GT..5) X=l.-X 



X=SQKT<-2.*AL0G(X> > 



TINORM=X-(A(3)+X*(A(2)+X*A(l)))/«l.+X*(B(3)+X»(B(2)+X* B(l)))> 



CALL OVERFL(I) 



IF(I.EQ.l) RETURN 2 



IFIALPHA.LT..5) TINORM=-TINORM 



RETURN 



RETURN 2 



END 



1 



- 



4 



13* 



A 



- 



2DI 



3D A 



ALOG 



- 



7 





ALPHA 



- 



1AG 



4 



B 



- 



201 



3D A 



I 



- 



3DA 



9AG 



OVERFL 



- 



9 





RETURN 



- 



10 



12 



SORT 



- 



7 





TINORM 



- 



1 



8= 



X 



- 



5= 



6 



-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦-♦—f— ••-■»—♦-♦-♦-+-♦-♦-♦-♦-♦-♦— ♦-+-♦-+-♦-♦-■•— ♦-♦-♦-+ 



45 



