sage: numerical_integral(lambda t:imaginary(gamma(1-I*t)),-15,15)
(0.0, 5.3925521851144085e-15)
On 18 Kwi, 00:02, Ian Petrow <ianpet...@gmail.com> 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 <ianpet...@gmail.com> 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 sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
