Opened 2 years ago

Last modified 6 weeks ago

#29142 new defect

numerical_integral sometimes returning (nan,nan) when x is polynomial variable

Reported by: gh-Ninja-Koala Owned by:
Priority: major Milestone: sage-9.6
Component: calculus Keywords:
Cc: Merged in:
Authors: Reviewers:
Report Upstream: N/A Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Status badges

Description

sage:R.<x>=RR[]
sage:f=-15.7724132052927*x^3 + 275.964230612533*x^2 - 1609.26342508113*x + 3127.78392476122
sage:g=2.34170721896176*x^3 - 41.4239781372561*x^2 + 244.135506079109*x - 479.217021253382
sage:t0=5.7642108710894
sage:t=5.764766515597
sage:numerical_integral(sqrt(diff(f)^2+diff(g)^2),t0,t)
(nan,nan)

It works if the first line is changed to

sage:var('x')

then the output is:

sage:numerical_integral(sqrt(diff(f)^2+diff(g)^2),t0,t)
(1.031355273859612e-06, 8.669172196746942e-20)

I'm using sage version 9.0 under Arch Linux My version of gsl is 2.6

Change History (5)

comment:1 Changed 21 months ago by mkoeppe

  • Milestone changed from sage-9.1 to sage-9.2

Moving tickets to milestone sage-9.2 based on a review of last modification date, branch status, and severity.

comment:2 Changed 15 months ago by mkoeppe

  • Milestone changed from sage-9.2 to sage-9.3

comment:3 Changed 9 months ago by mkoeppe

  • Milestone changed from sage-9.3 to sage-9.4

Moving to 9.4, as 9.3 has been released.

comment:4 Changed 5 months ago by mkoeppe

  • Milestone changed from sage-9.4 to sage-9.5

comment:5 Changed 6 weeks ago by mkoeppe

  • Milestone changed from sage-9.5 to sage-9.6
Note: See TracTickets for help on using tickets.