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

C	IVBOLU
C	DETERMINE A,B,ALPH,BETA,K12,K21,KEL AFTER IV ADM
C	OF BOLUS DOSE
C	ALSO CALCULATE AMOUNT OF DRUG IN THE TISSUE COMP.
C	VS TIME
C
C	NORMAL USAGE AS ANY LESQ FUNCTION
CC	IN THIS MODE CALCULATE CP=AE**(-ALPH*T)+B*E**(-BETA*T)
C	4 CONSTANTS IF 5 ENTERED USED AS DT ELSE DT =0
C	5 CONSTANTS DT CALCULATED
C	TO CALCULATE K21,K12,KEL CALL WITH NK=0
C	TO AVOID XT VS TIME PLOT ENTER BLANK LINE ON
C	LOWER,STEP,UPPER
	FUNCTION FUNCT(ARG,X,CP)
	DIMENSION ARG(1)
	REAL KEL,K21,K12
	COMMON IFEL,IX,NK
	A=ARG(1)
	B=ARG(2)
	ALPHA=ARG(3)
	BETA=ARG(4)
	T=X-ARG(5)
	IF (NK.EQ.0)GOTO 100
	E1=ALPHA*T
	E2=BETA*T
	IF(E1.GT.80.OR.E1.LE.0)IFEL=1
	IF(E2.GT.80.OR.E2.LE.0)IFEL=2
	FUNCT=CP-A*EXP(-E1)-B*EXP(-E2)
	RETURN
100	WRITE(0,1000)
	READ(0,2000)FL,FS,FM,DOSE
2000	FORMAT(10F10.0)
	K21=(A*ALPHA+B*BETA)/(A+B)
	KEL=ALPHA*BETA/K21
	K12=ALPHA+BETA-KEL-K21
	WRITE(3,3000)DOSE,K12,K21,KEL
3000	FORMAT('0DOSE=',G15.4,/,' K12  K21  KEL =',3G15.4)
1000	FORMAT(' LOWER,STEP,UPPER ,DOSE:',$)
	IF(FL.EQ.0)RETURN
10	XT=DOSE*K12*(EXP(-BETA*FL)-EXP(-ALPHA*FL))
     1/(ALPHA-BETA)
	WRITE(3,4000)FL,XT
4000	FORMAT(1H ,F15.1,G15.4)
	FL=FL+FS
	IF(FL.LE.FM)GOTO 10
	RETURN
	END