That's true! Sorry, I hadn't though it through.
Still puzzled by the concept, but OK, that's just me.
The original problem I was trying to solve had nothing to do with asymptotes.
It's a problem one may encounter in predicting the efficacy of an orally administered drug (I am a medicinal chemist).
If you administer an oral dose D
0 of a drug, it can be shown that under certain hypotheses its concentration in plasma C(t) is given by this function:
C(t)=\frac{D_0 F} {V} \frac{k_a} {k_a-k}(e^{-k t}-e^{-k_a t})
where F, k, k
a and V are constants related to the absorption, elimination and distribution of the drug in the subject.
Here's the plot of the function for a typical case (D
0=1 mg/kg, F=0.8, k=0.5 h
-1, k
a=1 h
-1, V=2 L/kg). The time t is expressed in hours, and C(t) is expressed in mg/L:
To get the desired pharmacological effect, you often want this concentration to be above a given value (let's call it 'EC' or effective concentration) for a sufficiently long time. There are other metrics of efficacy, but this one ('time over EC') is very common.
So, say that in this case EC is 0.1 mg/L. If you add the EC line to the above graph:
you can see that the time over EC is a bit less than 4 h.
Now, this is fine if you have an actual curve and one value of EC. You can use a graphical method to find the intersections between the EC line and the curve, etc.
However, what if you want to study theoretically (and quantitatively) what happens to the 'time over EC' when you change the parameters (k, k
a, V...).
As far as I know, to do that you should find an expression for 'time over EC' as a function of the parameters and EC, differentiate it w.r.t. each parameter and study the derivative.
'Time over EC' is clearly the difference between the two times where C(t) crosses the EC line.
That's where I got stuck. I don't think C(t)=EC can be solved analytically for t in general.
So my next step was to make an Excel file where you could see what happened to the curve when changing the various parameters.
That would be enough for most chemists, but I wanted to find an expression, at least an approximate one, for this flaming 'time over EC'.
So I tried to find an approximation to the C(t) curve that I could intersect analytically with EC.
The approach worked. I used the Pade approximation of the Taylor series centred on t=0 for the 'left' hand side of C(t) (before the max).
For the right hand side, I couldn't find a nice Taylor approximation, so I considered that Ln(C(t)) has a right asymptote that is reached quite early, and decided to intersect Ln(EC) with this asymptote. Here's the plot:
Obviously there is a small error in using the asymptote rather than the actual curve, but nothing I couldn't live with.
And here we come to the problem I presented here.
When k
a=k, C(t) can't be calculated with the usual function, it gives 0/0. If you take the limit of C(t) for k
a→k, you find:
C(t)=\frac{D_0 F} {V} k t e^{-k t}
which is fine with Taylor on the left hand side, but on the right hand side gives the mess we just discussed.
What I find utterly crazy is that if I 'cheat' and use the original C(t) with k
a≈k (e.g. 1.01 and 1, respectively), I get a curve that is almost identical to the one obtained with the new equation in k only, even in log scale:
Not unexpected, OK; but then I'm even more gobsmacked that the new equation should not allow asymptotes...
There you go, live and learn.
In practice, to approximate the 'time over EC' I think I'll use the curve where the constants are different, and in the cases where they are the same I will use two very close values for them. Quite distasteful, but what can you do...
Thanks again for your input, by the way.