C TEILPROGRAMM BARTLE VON VAR1
C MULTIPLER T-TEST UND BARTLETT-TEST
C FORTRAN II / OPTION (I)
C FS : PFT.FT / PCHI.FT
C CH : DUNCAN.SV
	COMMON NG,N1,FM,VS,NZ,FKW,FMW
	DIMENSION NG(49),FM(49),VS(49),NZ(49),NC(49)
	SPACE=' '
C MULTIPLER T-TEST
	DO 2 I=1,N1
2	NC(I)=NG(I)
	WRITE(1,100)LEER,
100	FORMAT('MULTIPLER T-TEST ? (JA=1,NEIN=0) : 'I0)
	READ(1,103)ID
103	FORMAT(I1)
	IF(ID)12,12,9
9	WRITE(1,104)LEER
104	FORMAT(I0)
	FMW=SQRT(FMW)
	IF(FKW-2047.)6,7,7
6	KW=IFIX(FKW)
	GOTO 8
7	KW=2047
8	DO 11 I=1,N1
	DO 10 J=1,N1
	IF(NC(J))10,10,3
3	IF(NC(J)-NG(I))5,4,5
4	NC(J)=0
	GOTO 10
5	IG=NG(I)
S	JMS NTD
C NAME TEMUXX
	IN1=1285
	IN2=0853
	IN3=IE
	CALL IOPEN('P1',IN1)
	READ(4,101)IOPT,KZ
101	FORMAT(A2)
	READ(4,102)EFA,EFA,PARA
102	FORMAT(A6)
	IG=NG(J)
S	JMS NTD
C NAME TEMUXX
	IN1=1285
	IN2=0853
	IN3=IE
	CALL IOPEN('P1',IN1)
	READ(4,101)KZ,KZ
	READ(4,102)EFB,EFB,PARB
	FNZI=FLOAT(NZ(I))
	FNZJ=FLOAT(NZ(J))
	T=SQRT(FNZI*FNZJ/(FNZI+FNZJ))
	T=T*(FM(I)-FM(J))/FMW
	TW=ABS(T)
	TT=100.*PFT(1,1,KW,TW)
	WRITE(1,105)NG(I),EFA,PARA,NG(J),EFB,PARB,TT
105	FORMAT(' #'I2' ('A6'/'A6') = #'I2' ('A6'/'A6
	1')    P ='F7.2' %')
10	CONTINUE
	WRITE(1,107)SPACE
107	FORMAT(A1)
11	CONTINUE
12	WRITE(1,107)SPACE
	FN=0.
	DO 15 I=1,N1
15	FN=FN+FLOAT(NZ(I))
	FT=0.
	DO 16 I=1,N1
16	FT=FT+1./FLOAT(NZ(I)-1)
	FT=FT-1./(FN-FLOAT(N1))
	C=1./(3.*FLOAT(N1-1))
	C=1.+C*FT
	UM=2.3025851
	C=UM/C
	FT=0.
	DO 17 I=1,N1
17	FT=FT+(FLOAT(NZ(I)-1))*ALOG(VS(I))/UM
	FMW=FMW*FMW
	FT=C*((FN-FLOAT(N1))*ALOG(FMW)/UM-FT)
	WRITE(1,111)
111	FORMAT('BARTLETT-TEST'/13('=')/)
	M=N1-1
	C=PCHI(M,FT)*100.
	WRITE(1,112)C
112	FORMAT('   SIGNIFIKANZ-NIVEAU FUER GLEICHHEIT DER STREU
	1UNGEN : 'F6.2' %'/)
	CALL CHAIN('DUNCAN')
SNTD,	0		/NUMMER DER DATEI
	IZ=(IG/10)*10
	IE=IG-IZ
	IZ=IZ/10
S	CLA CLL
S	TAD \IZ
S	CLL RTL;RTL;RTL
S	TAD \IE
S	TAD (6060
S	DCA \IE
S	JMP I NTD
	END