File WI.F4

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

	DIMENSION DY(7),PHID(7),GAM(7),XX(2),CC(2)
C------DATA VALUES
	PI=3.14159
	DTR=PI/180.
	RTD=1./DTR
	AL=3.289
	HEL=AL*.5
	WD=.014
	PSI=PI/2.
	GS=.5*(SQRT(5.0)-1.0)
	DO 45 I=1,7
45	DY(I)=.045-(I-1)*.015
	GAM(1)=31.5781
	GAM(2)=330.9301
	GAM(3)=63.4400
	GAM(4)=0.0
	GAM(5)=299.1013
	GAM(6)=91.0496
	GAM(7)=270.
	PHID(1)=14.2601
	PHID(2)=PHID(1)
	PHID(3)=24.2763
	PHID(4)=0.0
	PHID(5)=PHID(3)
	PHID(6)=PHID(1)
	PHID(7)=PHID(1)
C------FIND THE X WHERE THE SPACING IS A MINIMUM
	DO 16 II=1,6
	JJ=II+1
	KNT=1
	XS=0.0
	XF=1.5
	F1=PHID(II)*DTR*.5
	F2=PHID(JJ)*DTR*.5
	SN1=SIN(GAM(II)*DTR)
	SN2=SIN(GAM(JJ)*DTR)
	CS1=COS(GAM(II)*DTR)
	CS2=COS(GAM(JJ)*DTR)
C------GOLDEN SECTION SEARCH
	JN=3
11	DX=GS*(XF-XS)
	XX(1)=XF-DX
	XX(2)=XS+DX
	DO 7 J=1,2
	IF(J.EQ.JN) GO TO 7
	X=XX(J)
	CALL BN2(X,Y1,A1,EL,F1,HEL,PSI)
	CALL BN2(X,Y2,A2,EL,F2,HEL,PSI)
	DYY=Y1*CS1+DY(II)-Y2*CS2-DY(JJ)
	DZZ=Y1*SN1-Y2*SN2
	CC(J)=SQRT(DYY**2+DZZ**2)
7	CONTINUE
C-----INCREMENT COUNT AND QUIT IF OVER 20
	KNT=KNT+1
	IF(KNT.GT.20) GO TO 1
	IF(CC(1)-CC(2)) 8,9,10
C------C(1)<C(2) : LEFT SIDE OF INTERVAL
8	CC(2)=CC(1)
	JN=2
	XF=XX(2)
	GO TO 11
C-------C(1)=C(2) : CENTER OF INTERVAL
9	XS=XX(1)
	XF=XX(2)
	JN=3
	GO TO 11
C------C(1)>C(2) : RIGHT SIDE OF INTERVAL
10	CC(1)=CC(2)
	JN=1
	XS=XX(1)
	GO TO 11
1	C=.5*(CC(1)+CC(2))-.007*(1.0/COS(A1)+1.0/COS(A2))
16	WRITE(5,4) II,JJ,X,Y1,Y2,C
4	FORMAT(//,' WIRES ',2I1,'  X : ',F4.2,' Y1 : ',F8.6,
	2' Y2 : ',F8.6,'  C : ',F8.6)
	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