File BARTLE.FT (FORTRAN source file)

Directory of image this file is from
This file as a plain text file

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



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search