CC ORAL 1 C THIS PROGRAM DETERMINES KA,KE,DOS/VC(=DVC) & DT C FOR THE ONE COMPARTMENT MODEL. C DT=0 WHEN NO LOG TIME EXISTS. C C NORMAL USE AS ANY LESQ FUNCTION. C 4 CONSTANTS FUNCTION FUNCT(ARG,X,CP) DIMENSION ARG(1) REAL KA,KE COMMON IFEL,IX,NK KA=ARG(1) KE=ARG(2) DVC=ARG(3) T=X-ARG(4) E1=KE*T E2=KA*T IF (E1.GT.80.OR.E1.LE.0)IFEL=1 IF (E2.GT.80.OR.E2.LE.0)IFEL=2 IF (KA.EQ.KE)IFEL=3 IF (IFEL.NE.0)RETURN FUNCT=CP-KA*DVC*(EXP(-E1)-EXP(-E2))/(KA-KE) RETURN END A E1 FNEG FSTA #TMP+0003 JSR EXP JA .+0004 JA #TMP+0003 FSTA #TMP+0003 FLDA E2 FNEG FSTA #TMPX+0000 JSR EXP JA .+0004 JA #TMPX+0000 FNEG FADD #TMP+0003 FMUL #TMP+0000 FSTA #TMP+0000 FLDA KA FSUB KE FSTA #TMP+0003 FLDA