File LILI.FT (FORTRAN source file)

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

C TEILPROGRAMM LILI VON KURVEN
C KOORDINATEN : ABSZISSE : LINEAR / ORDINATE : LINEAR
C FORTRAN II / OPTION (I)
C SR : ERASE.FT / ALPHA.FT / HOME.FT / FDIS.FT / KAKO.FT
C CH : MRVA.SV / NP.SV
	COMMON NO,N1,NG,ITEXT,XA,XE,YA,YE
	DIMENSION NG(99),ITEXT(3,15),ABZW(99),M(4)
	READ(1,100)ID
100	FORMAT('MESSWERTE (1) ODER STATISTIK-WERTE (2) : 'I1)
	IF(ID-1)4,4,2
2	WRITE(1,101)
101	FORMAT(/'ABSZISSEN-WERTE DER MESSREIHEN')
	DO 3 I=1,N1
	WRITE(1,102)NG(I),LEER,
102	FORMAT('  MR # 'I2' : 'I0)
	READ(1,103)ABZW(I)
103	FORMAT(F15.6)
3	CONTINUE
	WRITE(1,104)
104	FORMAT(/' AUSWAHL STATISTIK-WERTE (JA=1,NEIN=0) :')
	READ(1,105)(M(J),J=1,4)
105	FORMAT('  ARITHMETISCHER MITTELWERT : 'I1/
	1'        MAXIMUM UND MINIMUM : 'I1/
	2'          XM+2*S UND XM-2*S : 'I1/
	3'          XM+3*S UND XM-3*S : 'I1)
4	ANZ=0.
	CALL KAKO(XA,XE,YA,YE)
	PAUSE
	DO 21 I=1,N1
	IK=1
	IG=NG(I)
S	JMS NTD
C NAME DER DATEI : TEMUXX
	IN1=1285
	IN2=0853
	IN3=IE
5	CALL IOPEN('P1',IN1)
	READ(4,109)IOPT,KZ
109	FORMAT(A2)
	READ(4,110)FIZ,EF,EF
110	FORMAT(A6)
	READ(4,109)(IL,K=1,KZ)
	GOTO(20,6),ID
6	GOTO(7,10),IK
7	CONTINUE
S	JMS AMWB
	IK=2
	GOTO 5
10	CONTINUE
S	JMS STRB
	WRITE(1,111)ABZW(I),YM,SS,YMI,YMA
111	FORMAT(5E12.5)
	GOTO 21
20	ANZ=ANZ+FIZ
S	JMS PKTE
21	CONTINUE
	CALL HOME
	WRITE(1,190)((ITEXT(I,J),J=1,15),I=1,3),ANZ
190	FORMAT(2(15A2,/),15A2,'   N = 'F9.0)
	PAUSE
	CALL ERASE
	CALL CHAIN('MRVA')
99	CALL CHAIN('NP')
SNTD,	0
C 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 (K6060
S	DCA \IE
S	JMP I NTD
SPKTE,	0
C MESSWERTE ALS PUNKTE DARSTELLEN
	DX=1000./(XE-XA)
	DY=700./(YE-YA)
	X=0.
1000	READ(4,110)Y
	IF(XA-X)1002,1002,1006
1002	IF(XE-X)1006,1003,1003
1003	IF(YA-Y)1004,1004,1006
1004	IF(YE-Y)1006,1005,1005
1005	FX=DX*(X-XA)
	FY=DY*(Y-YA)
	IX=FX
	IY=FY
	CALL FDIS(-1,IX,IY)
1006	X=X+1.
	FIZ=FIZ-1.
	IF(FIZ)1007,1007,1000
1007	CONTINUE
S	JMP I PKTE
SAMWB,	0
C ARITHMETISCHEN MITTELWERT BERECHNEN
	AZ=FIZ
	YM=0.
1010	READ(4,110)Y
	YM=YM+Y
	FIZ=FIZ-1.
	IF(FIZ)1011,1011,1010
1011	YM=YM/AZ
S	JMP I AMWB
SSTRB,	0
C STREUUNG UND YMAX,YMIN
	FIZ=FIZ-1.
	AZ=FIZ
	READ(4,110)Y
	SS=(Y-YM)**2
	YMI=Y
	YMA=Y
1012	READ(4,110)Y
	SS=SS+(Y-YM)**2
	IF(Y-YMI)1013,1014,1014
1013	YMI=Y
1014	IF(Y-YMA)1016,1016,1015
1015	YMA=Y
1016	FIZ=FIZ-1.
	IF(FIZ)1017,1017,1012
1017	SS=SQRT(SS/AZ)
S	JMP I STRB
	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