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: |
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
- Milestone changed from sage-9.1 to sage-9.2
comment:2 Changed 15 months ago by
- Milestone changed from sage-9.2 to sage-9.3
comment:3 Changed 9 months ago by
- 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
- Milestone changed from sage-9.4 to sage-9.5
comment:5 Changed 6 weeks ago by
- Milestone changed from sage-9.5 to sage-9.6
Note: See
TracTickets for help on using
tickets.
Moving tickets to milestone sage-9.2 based on a review of last modification date, branch status, and severity.