Numerical analysis (composite numerical integration)

AI Thread Summary
The discussion focuses on calculating the error bound for the composite trapezoidal rule when integrating the function tan(x) from 0 to π/4. Participants clarify that the second derivative of tan(x) is 2sec²(x)tan(x), and the maximum value occurs at π/4, yielding a bound of 4 for the error calculation. The error formula is applied using h = π/16, leading to the expression for error being less than (π/4)/(12) * (π²/16²) * 4. Additionally, for the function exp(x)sin(x) over the interval from 1 to 3, the maximum of the second derivative is determined to be at x = 3, which is used for error estimation. The conversation emphasizes the importance of identifying maximum values of derivatives for accurate error bounds in numerical integration.
willbe
Messages
7
Reaction score
0
using composite trapezoidal rule with n=4 how can i get a bound for the error of I=integration tan(x) from x=0 to x=pi/2

i know that the term of error in comp trapezoidal rule is (b-a)/12 h^2 f''(eita)
i got the second derevative of tanx to be 2sec^2 x tanx then i don't know with what value exactly i need to substitute in this function to get the bound.
 
Physics news on Phys.org
Since tan(pi/2) is infinite, you cannot perform this integration.
 
sorry it is to x=pi/4
should i substitute by sec^2x by (1+tan^2x)
 
willbe said:
sorry it is to x=pi/4
should i substitute by sec^2x by (1+tan^2x)

You must have 0 < \eta < \pi/4

If you want a worst case error, differentiate again and find the value of \eta that maximizes f''.
 
ok i did and the second derivative is 2sec^2 x tanx now is the maximum value of this function is 4 as sec^2(x)=tan^2(x) + 1 ?
 
Taking the derivative of f'' and setting it to zero gives an immaginary root. So, there is no point of maximum.
The value of f'' at \pi/4 is 1, so we take it as the worst case.
h = \pi/4/4 = \pi/16
Then your error will be less then \frac{\pi/4}{12}\cdot \frac{\pi^2}{16^2}.
 
why we set the value of f'' to zero?
we set it to a value that gives us the maximum of the function so we substitute in 2sec^2(x)tan(x) which is equal to 2(1+tan^2(x))tan(x) by PI/4 so we get 2(1+1) which is 4.
right??
 
willbe said:
why we set the value of f'' to zero?
we set it to a value that gives us the maximum of the function so we substitute in 2sec^2(x)tan(x) which is equal to 2(1+tan^2(x))tan(x) by PI/4 so we get 2(1+1) which is 4.
right??

I have set the value of f''' to zero. The maximum of f'' happens when its derivative is zero.
Since f''' has no real roots the maximum of f'' must be at one of its boundary points 0 or \pi/4. The value of the funtion at 0 is 0 and at \pi/4 is 4 (and not 1 as I wrote previously). So in the expression for the error, you replace f''(\eta) by 4, getting \frac{\pi/4}{12}\cdot \frac{\pi^2}{16^2}\cdot 4
 
ok what about f(x)=exp x sinx using composite trapezoidal rule what is the bound of error for integration from 1 to 3

the second derevative will be 2exp x cos x then to get the bound of error which value should i substitute with??
 
  • #10
willbe said:
ok what about f(x)=exp x sinx using composite trapezoidal rule what is the bound of error for integration from 1 to 3

the second derevative will be 2exp x cos x then to get the bound of error which value should i substitute with??

Since exp x is monotone and grows very fast, the maximum value of f'' is at x = 3. This is the value to substitute in the expression.
 

Similar threads

Replies
7
Views
2K
Replies
3
Views
2K
Replies
11
Views
4K
Replies
1
Views
3K
Replies
3
Views
1K
Replies
13
Views
2K
Back
Top