File PT.FT (FORTRAN source file)

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

	FUNCTION PT(N,X)
C WAHRSCHEINLICHKEIT VON ST(N) = X
C ST = STUDENT-VERTEILUNG
C  N = FREIHEITSGRADE
C NACH DECUS NO. FOCAL 8-308
C FORTRAN II
	X=X*X+.1E-15
	ID=2*(N/2)-N+2
	W=X/FLOAT(N)
	Z=1./(1.+W)
	IF(ID-1)1000,1001,1000
1000	PT=SQRT(W*Z)
	D=.5*PT*Z/W
	GOTO 1002
1001	PT=SQRT(W)
	Y=.31830989
	D=Y*Z/PT
	PT=2.*Y*ATAN(PT)
1002	Y=2.*W/Z
	J=ID+2
1003	IF(J-N)1004,1004,1005
1004	D=D*Z*(1.+1./FLOAT(J-2))
	PT=PT+D*Y/FLOAT(J-1)
	J=J+2
	GOTO 1003
1005	RETURN
	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