C	IVCON1
C	THIS PROGRAMME DETERMINES A,B,ALPHA,BETA ,K12,K21 & KEL
C	FROM THE POST INFUSION PLASMA LEVELS
C	AFTER CONSTANT INFUSION OVER A SHORT PERIOD, TAU.
C
C	NORMAL USE AS ANY LESQ FUNCTION
C	4 CONSTANTS
C
	FUNCTION FUNCT(ARG,X,CP)
	DIMENSION ARG(1)
	REAL K12,K21,KEL
	COMMON IFEL,IX,NK
	A=ARG(1)
	B=ARG(2)
	ALPHA=ARG(3)
	BETA=ARG(4)
	TAU=ARG(5)
	IF (NK.EQ.0)GOTO 100
	E1=ALPHA*TAU
	E2=BETA*TAU
	E3=ALPHA*(T-TAU)
	E4=BETA*(T-TAU)
	IF(E1.GT.80.OR.E1.LE.0)IFEL=1
	IF(E2.GT.80.OR.E2.LE.0)IFEL=2
	IF(E3.GT.80.OR.E3.LE.0)IFEL=3
	IF(E4.GT.80.OR.E4.LE.0)IFEL=4
	IF(IFEL.NE.0)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