File PROINH.FT (FORTRAN source file)

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

C	PROINH
C	WITH INIBITOR
C		Y=SUM(<VIKT*X1*10***K>/<1+X1*10***K+X2*10**K2>
C
C	K=BINDING CONST=VIKT(NK+I)
C	K2=INHIB CONST TO BE D=ARG(I)ETERMINED
C	X=CONS DRUG
C	X2=CONS INHIBITOR
C	VIKT=RELATIVE WEIGHTS OF THE TERMS=VIKT(I)
C
C	MUST BE INITIALIZED<GIVES ERROR 1 OTHERWAYS>
C	INITILIZING<NK=0>
C	INPUT NK-1 WEIGHTS
C	INPUT NUMBER IN RESTTERM
C	INPUT NK VALUES FOR K
C	IF UNINHIBITED TERM IS WANTED ADD 2 KONSTANTS
C	WEIGHT OF TERM AND BINDING CONSTANT
	FUNCTION FUNCT(ARG,X,Y)
	COMMON IFEL,IX,NK,NX,DUMM(60),X2(30)
	DIMENSION ARG(1)
	DIMENSION VIKT(12)
	DATA VIKT(1)/0.0/
	Z2(U)=10.**VIKT(NK+I)*X/(1+10.**VIKT(NK+I)*X+10.**U*X2(IX))
	IF (NK.EQ.0)GOTO 100
	IF (VIKT(1).EQ.0)IFEL=1
	A=0
	NK1=NK-1
	IF (NK.EQ.1)GOTO 11
	DO 10 I=1,NK1
10	A=A+Z2(ARG(I))*VIKT(I)
11	I=NK
	A=A+VIKT(NK)*Z2(ARG(NK))
	IF (VIKT(2*NK+1).NE.0)A=A+VIKT(NK+2)*(10.**VIKT(NK+1)
     1*X)/(1.+10.**VIKT(NK+1)*X)
	FUNCT=Y-A
	RETURN
100	WRITE(0,1000)
	READ(4,2000)VIKT
	WRITE(3,3000)VIKT
	RETURN
1000	FORMAT(' WEIGHTS,NUMBER,BIND CONST <F6.0>:')
2000	FORMAT(12F6.0)
3000	FORMAT(1H0,12F8.2)
	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