Computer Software 153
PROGRAM:
NORMAL.FOR
c
************************************************************
C
FITTING
A
NORMAL
DISTRIBUTIO
N
TO
LOGS
OF
OBSER
VATIONS
AND
C
TESTING:
NORMAL.FOR
C
FOR
PROCESSI
NG
ACCUMULATION
/THICKNESS
/GRADE
AS
THE
CASE
MAY
BE
C
THIS
PROGRAM
READS
ONE
VA
LU
E
PER
RECORD
AND
PROCESSES
.
C
************************************************************
CHARACTER
*1
PNE
W
DIMENSIO
N A(31) ,AN(104) ,FR(100) ,AV(100) ,FL(100) ,EL(100)
DIMENSION
PCR(100) ,ALE
VE
L(6)
DIMENSION
R(2) , T(2 ) , P (32) ,PP(2) ,Q(52) ,E(52)
C
C
C
************************************************************
DISCRETE
VALUES
FOR
COMPUTING
THE
NORMAL
PROBABILITY
************************************************************
1
2
3
DATAA
/1
.,6.0,6.6666
,8 .4
,10
.2857,12.2222,14.1818
,16.1538
,
18.1333,2 0.1176,22 .1 052,24 .0952,26 . 0869,28 . 08,3 0. 0740,32 . 0689,
34.0645,36.0606,38 .0571,40.0540,42 .0512,44 .
048
7,46 .0465,48.
0444
,
50.0425,52
.0408,54.0392,56.0377,58
.0363,60.0350,62.03/
C
************************************************************
1111
C
C
C
5
51
52
111
PNE
W
=CHAR
(12 )
OPE
N
(UNIT
=6,FILE =
'NORMAL
.DAT'
,STATUS
='OLD')
OPE
N
(UNIT
=2,FILE =
'NORMAL.
RES
, , STATUS='
UN
KN
O
WN
' )
FORMAT
(A)
************************************************************
READS
CLASS
LIMITS
AN
(1)
AND
AN
(2)
FOR
ACCUMULATION
************************************************************
READ(6
,51)
(ALEVEL(I)
, 1=1, 6)
FORMAT
(6A4)
WRI TE(2, 52)
(ALEVEL(I)
,1 =1,6)
FORMAT
(5X,
6A4
/ )
READ(6,111)N1
,KEY1 ,KEY2
FORMAT
(312)
READ(6
,*)AN(1) ,AN(2) ,W
N2
=N1*2
************************************************************
************************************************************
************************************************************
N2
STANDS
FOR
THE
NUMBER
OF
CLASS
INTERVALS
C IS A
CO
NSTANT
TO
BE
ADDED
IN
CASE
THE
DISTRIBUTION
IS A 3-
PARAMETER
LOG
NORMAL
DISTRIBUTION
. C
CAN
BE
0
ALSO.
CASES
CASES
USUALLY
AN(l)
=O
.O,
AN(2)
=0.5
MAY
SUFFICE
FOR
ALL
OF
ACCUMULATIO
N
AND
AN(l)
=O.O
AND
AN(2)
=0.2
FOR
ALL
OF
GRADE/THICKNESS
.
DO
100 K=3,N2 ,2
AN(K)
=AN(K
-1)
AN(K+1)
=AN(K)+
W
CONTINUE
C
C
C
C
C
200
100
C
C
C
C