On Sun, Jun 16, 2013 at 10:29 AM, Shriramana Sharma <[email protected]> wrote: > 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.
Ah, I scanned through the code, but all I saw was subs. But I see the call to limit now. I guess I would note the structure of your expression when factored (look at factor(kappa.subs(mybez1))). Be careful about applying certain power identities (see http://docs.sympy.org/tutorial/tutorial/simplification.html#powers). Aaron Meurer > > @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. > > -- 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.
