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
DO 250 I=1,NK
250	ARG(I)=XM(I)
	SMAX=H
	GOTO 1
500	DO 600 I=1,NK
	DO 600 J=1,NK
	S=0.
	DO 520 K=1,NX
	I3=K+(I-1)*NX
	I4=K+(J-1)*NX
520	S=S+A(I3)*A(I4)
	B2(I,J)=S
600	B2(J,I)=S
C	S MINV FROM PERSSON
	I3=NK-1
	N2=0
	DO 611 L=1,NK
	AUX(NK)=1.0/B2(1,1)
	N2=N2+1
	DO 612 K=1,I3
	KP1=K+1
612	AUX(K)=B2(1,KP1)*AUX(NK)
	DO