First time here? Check out the FAQ!

Ask Your Question
1

numerical_integration returns two numbers. What is the second?

asked 8 years ago

Alan gravatar image

The command numerical_integral returns two numbers. In the documentation, under "numerical integration", the second is called an "error bound", and later on it is called an "error estimate". What is this number, and precisely how is it defined?

Preview: (hide)

2 Answers

Sort by » oldest newest most voted
2

answered 8 years ago

vdelecroix gravatar image

The second number is an estimation of an error bound. You could have seen it looking at the documentation as follows

sage: numerical_integral?
Docstring:
   Returns the numerical integral of the function on the interval from
   a to b and an error bound.
...
Preview: (hide)
link
2

answered 8 years ago

numerical_integral calls routines in the GNU Scientific Library. As described in the relevant introduction, the library returns a numerical integration result and an absolute error that satisfy

|RESULT - I| <= ABSERR <= max(epsabs, epsrel |I|)

where I is the exact value of the integral. The absolute and relative error inputs on the right both default to 1e-6 in SageMath.

Preview: (hide)
link

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 8 years ago

Seen: 417 times

Last updated: Jun 11 '16