EXAMPLE PROGRAM "PRICLAND-2" EXAMPLE PROGRAM "PRICLAN -2" 
  
  
  
IX DIMENSION ABSVAR(200),SVVAR( 200) DIMENSION K2(25),K4(25),IXA$900,24) , 
MIR * GUBX ( 300) , TRENX ( 300) ,11 (77) ,£ (300) *18(900) ,Ix (9005 ,1R(900) ,1c(900), 
40 C DO 44 KK=1,50 +15 (900) 11¢1(9003 ,1¢2(900) ,X1 (900) , 
rere READ(S,20)N, TT * AK2(900, 24) ,K2(900, 24.) 
78,05 20 FORMAT(I3,77A1) READ(5,1)N,M,N2,N3 
20.67 PRINT 20,N,TT PRINT 4,N,M,N2,N2 
119465 READ 200,(X(1),1=1,N) 1 FORMAT(4I3) 
4^ 5 O4, PRINT 200,(X(1),1-1,N) DO 2 1=1,8 
1453 200 FORMAT(16F5.0) READ (5, 33 IE(T) 3 (IXA(T J) ,J24,M) 
1 «99 SYMzO, 2 CONTINUE 
2,54 XMIN=X(1) 3 FORMAT(I4,24#I2) 
5,97 XMAX=XMIN DO 4 I=1,M 
3427 DO 97 I=1,N WRITEGO,2)IE(I),(IXA(I,J),J21,M) 
26,02 ICE EE 4 CONTINUE 
0,54. IF(X(I).GT.XMAX)XMAX-X(I) DO 66 J=1,M 
0,29 77 CONTINUE MIN=IXA(1,J 
1.04 WRITE(6,71) MAX =MIN 
0.57 PRINT 88,XMIN,XMAX DO 67 J=1,M 
4437 88 FORMAT(1H ‚ XMIN=",F10.2, RA D T ETAT) 
20,52 + "XMAX? ,840.2) IF(IXA(I,J).GT.MAX )MAXzIXA(I,J) 
5415 DO 4 I=1,N 67 CONTINUE 
momie SYM=SYM+X (I) K3(J)=MIN 
4 CONTINUE K4(J)=MAX 
tion of AN=N 66 CONTINUE 
AER SR=SYM/AN WRITE(6,42) 
BE WRITE(6,71) WRITE (6,667) 
tories WRITE(6,400)SYM,SR WRITE(6,96)(K3(J),J=1,M) 
LOOFORMAT(1H, , SYM=",F11.4, MIDDIE °, WRITE(6,42) 
* ARTTHMETICAL SR= ',FA1.4) WRITE (6,666) 
AVAR=0, RO EE 1 132300 
AKVAR=0, WRITE (6,42 
AINFOR=04 WRITE(6,777) 
CUBSUM=0. 1000 DO 5 J=1,M 
TREND=O, DO 6 IzA,N 
DO 23 I=1,N IX(1)=IXA(I,J) 
DOLJA=X (19/SYM 6 CONTINUE 
AINFOR=ATNFOR+3. 3222* ALOGTO(DOLJA)* DT2=N2 
*DOLJA DO 8 I=1,N 
ABSVAR(I)-X(I)-SR IR(I)=LT2 
AVAR=AVAR+ABSVAR(I)**2 8 CONTINUE 
SVVAR(I)-ABSVAR* $2 ISWA=ISWA+IR(I) 
CUBX (1 )=ABSVAR** 3 IS(I)=ISWA 
TRENX (1)=ABSVAR**4 IF(I.GT.4)GO TO 62 
CUBSUM=CUBSUM+CUBX (I) IC(I)=K3(J) 
TREND=TREND+TRENX (1) GO TO 9 
AKVAR=AKVAR+SVVAR(I) 62 I1=I-1 
23 CONTINUE IC(I)=K3(J)*IS(I1) 
DOTK=SQRT(AKVAR/ (AN=-1,0)) 9 CONTINUE 
ASIM=AN* CUBSUM/ (AN-=1 5 0)* (AN-2. )* DO 60 I=1,N 
— *DOTK** 3 ICA(I)=IC(I) 
ution of EXES=TREND/AN* DOTK* * 4 -5. IC2(I)=IC(I+41) 
ation of HDAszSQRT(6.* (AN-1, )/ (AN41, )* (AN425)) 60 CONTINUE 
andscape HDE-ZSQRT(24.,* AN* (AN-2, )* AR )/ DO 220 I=1,N3 
* ((AN=10)**2)* (AN+3. )* (AN+5. ) KT=0 
TA=ASIM/HDA DO 55 K=1,N 
TE=EXES /HDE IF(IX(K).GE,ICA(I).^ AND,IX(K) LT, 
AMAXINz2, 2222* AIOGAO(AN) *IC2(I1)))80 TO 50 
OTINFZ(AMAXIN-ABS (AINPFOR))/AMAXIN NX=0 
PRINT 75,AINFOR,AMAXIN,OTINF GO TO 40 
75 FORMAT(1H ,°GENERAL ENTROPY =°,F10. 50 NX=1 
*F10,4,10x, "MAXIMAL ENTROPY =", 40 K1(K)=NX 
*P10.4,5X, RELATIVE ENTROPY =,F10.4) KT=KT+K1(K) 
CV=DOTK/SR 55 CONTINUE 
WRITE(6,888)ASIM,EXES ,HDA,HDE,DOTK, K1(1)=KT 
*CV | : K2(1,J)=K1(1) 
888 FORMAT(1H , “COEFFICIENT ASIMMETRY”, 220 CONTINUE 
** KA=",F10,4/14 , “COEFFICIENT ^,. 5 CONTINUE 
* "EXCESS KEz^,P40,A/4H , REPRESEN’, DO 116 J=1,M 
* “PATIVITY COEFFICIENT ASSIMMETRY ^, K2(1,J)=K2(1,J)=1 
** HDAz^,FA0.4/4H ,"REPRESENTATI^, 416 CONTINUE 
* 'VITY COEFFICIENT EXCESS HDEz^ DO 40 Iz1 + 
F F10,4/1H , "MIDDLE SQUARE VARI&®, WRITE (6,402 ICAQG),GO2G,0) ,4921,M), 
* *'TION DOTK-",F45.4/4H ,'COEFFI", *7C(1) 
* CIENT OF VARIATION CV=",F10.4) 40 CONTINUE 
END WRITE(6,777) 
49