File TWCOM1.FT (FORTRAN source file)

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

C	TWCOM1
C	4 CONST<T=X-DT>
C	OR
C	3 CONST<DT=0>
C	FUNCTION:
C
C	CP=[KA*DOS/<VC(KA-KE)>]*[EXP(-KE*<T-DT>)-EXP(-KA*<T-DT>)]
C	INITIAL ESTIMATES FOR :
C	KA  K21  D/VC  AND DT


	FUNCTION FUNCT(ARG,X,Y)
	DIMENSION ARG(1)
	COMMON IFEL
C	ERRORTESTS FOR 1)KA=KE
C		       2)EXPONENTS > 100
	FKA=ARG(1)
	FKE=ARG(2)
	IF(FKA.EQ.FKE)GOTO 1
10	T=X-ARG(4)
	E1=FKA*T
	IF(E1.GT.85)IFEL=2
20	E2=FKE*T
	IF(E2.GT.100)IFEL=3
	IF(E1.LE.0.OR.E2.LE.0)GOTO 4
	FUNCT=FKA*ARG(3)/(FKA-FKE)
	FUNCT=FUNCT*(EXP(-E2)-EXP(-E1))
	FUNCT=Y-FUNCT
	RETURN
4	IFEL=3
1	IFEL=IFEL+1
	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