#17606: Calculation involving nintegral works with Sage 6.3 but not with Sage
6.4
-------------------------+----------------------------------
Reporter: alyasoff | Owner:
Type: defect | Status: new
Priority: major | Milestone: sage-6.4
Component: numerical | Keywords: nintegral
Merged in: | Authors: Andrew Lyasoff
Reviewers: | Report Upstream: N/A
Work issues: | Branch:
Commit: | Dependencies:
Stopgaps: |
-------------------------+----------------------------------
The following code:
{{{
var('u')
assume(u>0
}}}
{{{
fu=1/8*(I*sqrt(2)*(erf(1/2*sqrt(2)*(I*pi + 2*u + 1)) -
erf(1/2*sqrt(2)*(-I*pi + 2*u + 1)))*e^(u + 1/2) +
sqrt(2)*(-I*erf(1/2*sqrt(2)*(I*pi - 2*u + 1)) + I*erf(1/2*sqrt(2)*(-I*pi -
2*u + 1)))*e^(u + 1/2) - (sqrt(2)*e^(1/8) - sqrt(2)*e^(2*u +
1/8))*(erf(1/4*sqrt(2)*(2*I*pi + 4*u + 1)) + erf(1/4*sqrt(2)*(-2*I*pi +
4*u + 1))) - (sqrt(2)*e^(1/8) - sqrt(2)*e^(2*u +
1/8))*(erf(1/4*sqrt(2)*(2*I*pi - 4*u + 1)) + erf(1/4*sqrt(2)*(-2*I*pi -
4*u + 1))))*e^(-u)/sqrt(pi)
}}}
{{{
fu.nintegral(u,0,2)
}}}
returns:
{{{
(1.4275778409741315, 3.656408002774788e-14, 63, 0)
}}}
with Sage 6.3, but with Sage 6.4 returns:
{{{
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "_sage_input_4.py", line 10, in <module>
exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8
-*-\\n" +
_support_.preparse_worksheet_cell(base64.b64decode("ZnUubmludGVncmFsKHUsMCwyKQ=="),globals())+"\\n");
execfile(os.path.abspath("___code___.py"))
File "", line 1, in <module>
File "/tmp/tmpbwhja8/___code___.py", line 3, in <module>
exec compile(u'fu.nintegral(u,_sage_const_0 ,_sage_const_2 )
File "", line 1, in <module>
File "sage/symbolic/expression.pyx", line 10244, in
sage.symbolic.expression.Expression.nintegral
(build/cythonized/sage/symbolic/expression.cpp:51090)
File "/usr/lib/sagemath/local/lib/python2.7/site-
packages/sage/calculus/calculus.py", line 792, in nintegral
raise TypeError(err)
TypeError: ECL says: Error executing code in Maxima:
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/17606>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" 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/sage-trac.
For more options, visit https://groups.google.com/d/optout.
