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
~} |g3D.Db@{g%(/1@/z)|5v$`/'X.yo 5bx)XyipwrvXyrxu@eztf|sY$b(ry
Hoq)+&()&'(&&'&%&&5h/po9dsi5@/]5"`n4/`4&sU4" oy3b nyml9m=<#a
Sm\2L7WL:qDW@DW@-
K]VTE