Opened 5 years ago

Last modified 4 months ago

#17606 new defect

Calculation involving nintegral works with Sage 6.3 but not with Sage 6.4

Reported by: alyasoff Owned by:
Priority: major Milestone: sage-6.4
Component: calculus Keywords: nintegral
Cc: Merged in:
Authors: Andrew Lyasoff Reviewers:
Report Upstream: N/A Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Description (last modified by jdemeyer)

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: 

Attachments (1)

Sage-6.3-6.4-diff.sws (7.9 KB) - added by alyasoff 5 years ago.

Download all attachments as: .zip

Change History (4)

Changed 5 years ago by alyasoff

comment:1 Changed 5 years ago by jdemeyer

  • Description modified (diff)

comment:2 Changed 5 years ago by rws

  • Component changed from numerical to calculus

comment:3 Changed 4 months ago by chapoton

works fine in 8.9.b7

Note: See TracTickets for help on using tickets.