Wow. Great. Thanks. Ian
On Apr 17, 3:12 pm, achrzesz <[email protected]> wrote: > sage: numerical_integral(lambda t:imaginary(gamma(1-I*t)),-15,15) > (0.0, 5.3925521851144085e-15) > > On 18 Kwi, 00:02, Ian Petrow <[email protected]> wrote: > > > > > > > > > Also, imaginary(gamma(1.+I*5)), say, works fine as long as it's not > > inside a numerical integral > > > Ian > > > On Apr 17, 2:57 pm, Ian Petrow <[email protected]> wrote: > > > > Hello Everyone, > > > > I know numerical integration isn't great in SAGE, but I think the > > > following isn't too complicated. I try > > > > sage: numerical_integral(imaginary(gamma(1-I*t)),-15,15) > > > > which fails, and give a LONG error message. I won't put the whole > > > thing here, but it starts: > > > > ERROR: An unexpected error occurred while tokenizing input > > > The following traceback may be corrupted or invalid > > > The error message is: ('EOF in multi-line statement', (35, 0)) > > > > --------------------------------------------------------------------------- > > > TypeError Traceback (most recent call > > > last) > > > > /Users/ianpetrow/<ipython console> in <module>() > > > > /Applications/sage/local/lib/python2.6/site-packages/sage/gsl/ > > > integration.so in sage.gsl.integration.numerical_integral (sage/gsl/ > > > integration.c:1622)() > > > > /Applications/sage/local/lib/python2.6/site-packages/sage/symbolic/ > > > expression.so in sage.symbolic.expression.Expression._fast_float_ > > > (sage/symbolic/expression.cpp:29478)() > > > > /Applications/sage/local/lib/python2.6/site-packages/sage/symbolic/ > > > expression_conversions.pyc in fast_float(ex, *vars) > > > 1212 1.4142135623730951 > > > 1213 """ > > > -> 1214 return FastFloatConverter(ex, *vars)() > > > 1215 > > > 1216 ################# > > > > etc etc etc, and says at the end > > > > TypeError: cannot convert I to real number > > > > Whereas on the other hand, > > > > numerical_integral(imaginary(1+I*t),-15,15) > > > (0.0, 2.4968976137148349e-12) > > > > works just fine. What gives? I'm running version 4.6 in OS X. > > > > Thanks, > > > > Ian -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
