Opened 5 years ago

Last modified 2 years ago

#16905 new defect

GSL can't handle really long intervals for integration?

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

Description

f1(x)=1/sqrt(x^3+2)
f2(x)=1/sqrt(x^4+2)
sage: numerical_integral(f2,1,10^8)
(0.8815690504421161, 3.309409685784312e-09)
sage: numerical_integral(f2,1,10^9)
(0.8815690594421439, 2.7280605832086615e-08)
sage: numerical_integral(f2,1,10^10)
(0.8815690603426408, 6.194229607849825e-07)
sage: numerical_integral(f2,1,10^11)
(2.3214916598860602e-07, 4.5569931705290324e-07)

See here for first report.

Change History (2)

comment:1 Changed 5 years ago by kcrisman

Maxima (the .nintegrate method) seems to have a similar problem.

Last edited 5 years ago by kcrisman (previous) (diff)

comment:2 Changed 2 years ago by jakobkroeker

  • Cc jakobkroeker added
Note: See TracTickets for help on using tickets.