Hello Aaron and Chris and thanks for your replies. @Aaron: Perhaps you didn't see the script since I've already used SymPy's limit in it. The question was how to do it *without* limit (i.e. for implementation in C++ where SymPy doesn't exist). And I read up a bit on the Gruntz method, but probably won't require such a sledgehammer to tackle my limited problem.
@Chris: That's a nice trick, but I'm still getting a NaN because both the numeration and denominator of the curvature function evaluate to zero. And even by trying the same trick for the successively differentiated numerators and denominators (as per L'Hopital) I ran up against the same problem as before, i.e. when the numerator becomes non-zero, the denominator becomes a NaN. Perhaps I should sub-apply L'Hopital to the denominator to resolve the NaN status. Thanks anyhow... -- Shriramana Sharma ஶ்ரீரமணஶர்மா श्रीरमणशर्मा -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. For more options, visit https://groups.google.com/groups/opt_out.
