Finding Error Bounds for Simpson's Method in Approximate Integration

Telemachus
Messages
820
Reaction score
30

Homework Statement


Use the Simpson method to estimate \displaystyle\int_{0}^{1}\cos(x^2)dx with an approximation error less than 0.001.

Well, I have a problem. Actually I'm looking for a bound for the error of approximate method integration by using Simpson's method.

I have to bring \displaystyle\int_{0}^{1}\cos(x^2)dx with an error less than 0.001.

I started looking for the fourth order derivative, and got:

f^4(x)=48x^2\sin(x^2)-12\cos(x^2)+16x^4\cos(x^2)

Now, I have to find a bound K for this derivative in the interval [0,1]. What I did do was watch for if they had maximum and minimum on the interval, then calculate the derivative of order five:
f^5(x)\in{[0,1]}\Rightarrow{f^5(x)=0\Longleftrightarrow{x=0}}

From here I did was ask:

If f^5(x)\in{[0,1]}\Rightarrow{f^5(x)=0\Longleftrightarrow{x=0}}

So here I know is that zero is a maximum or minimum. Then I wanted to look for on the concavity of the curve, and I thought the easiest thing would be to look at the sixth derivative, to know how it would behave the fourth derivative of the original function.

f^6(x)=720x^2\cos(x^2)-480x^4\sen(x^2)-64x^6\cos(x^2)+120\sen(x^2)

The problem is that when evaluated

f^6(0)=120\sen(0^2)=0

So, I get zero the derivative sixth, and I do not say whether it is concave upwards or downwards is concave.

What should I do? there may be a less cumbersome to work this, if so I would know.

Greetings.
 
Last edited:
Physics news on Phys.org
You could just plot the function and see it hits a maximum on that interval at x=1. Alternately, you can use arguments like x2sin x2≤x2 to find an upper bound for the function. In these types of problems, it's usually not terribly important to find the actual maximum; an upper bound for the maximum is good enough.
 
I complicated it too much. Just by evaluating the 6th derivative in any other point of the interval I could get the concavity of the curve. An other way was just seeing for the value of the 5th derivative, as the concavity changes only in zero then I would know if it creases or decreases.

Bye, and thanks.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top