How to Solve an Improper Integral with High Precision Using Simpson's Rule?

AI Thread Summary
To solve the improper integral from 10 to infinity of (sin(1/x)/(1+x^3))dx with a precision error of e=0.5*10^-4, the discussion emphasizes using Simpson's Rule. Participants suggest breaking the interval into manageable subintervals and estimating the integral on each to ensure the total error remains within the specified limit. A substitution of u=1/x simplifies the integral, allowing for easier evaluation. Numerical approximations using software like Mathematica yield results around 0.000329176 for the integral from 10 to 45 and 0.000332717 for 10 to 142. Ultimately, Simpson's Rule can achieve the desired accuracy when applied correctly.
undefined83
Messages
8
Reaction score
0
Im supposed to solve
integral 10 to +infinity ((sin(1/x)/(1+x^3))dx with error precision of e=0.5*10^-4. Can someone please give me detailed explenation of solving this. (Supposedly by Simpson but i get lost in the way.

P.S. sorry for bad spelling and lack of proper formula notions.
 
Physics news on Phys.org
I process from integral i get
1/3[ln((x+1)/sqrt(x^2-x+1))+sqrt(3)*arctg(2*sqrt(3)*x/3-sqrt(3)/3)] from M to +infinity <=1/4*10^-2

than i get

-ln((M+1)/(sqrt(M^2-M+1))) +sqrt(3)*(pi/2-arctg(2*sqrt(3)*M/3-sqrt(3)/3))<=7.5*10^-3

and i can't find any exact solution to solve that
 
undefined83 said:
Im supposed to solve
integral 10 to +infinity ((sin(1/x)/(1+x^3))dx with error precision of e=0.5*10^-4. Can someone please give me detailed explenation of solving this. (Supposedly by Simpson but i get lost in the way.

P.S. sorry for bad spelling and lack of proper formula notions.
break the interval up
10=x0<x1<x2<...<xn-1<xn=infinity
for
int(10,x1)
int(x1,x2)
.
.
.
int(xn-2,xn-1)
use simpsons rule
for
int(xn-1,infinity)
chose xn-1 large
estimate the integral on each subinterval accurate enough so that total error is within limit
 
Thanks, but i am totaly lost when i need to find int(xn-1,infinity). I try to make it lower than error estimate but i get either logaritmic inverse-trigonometric inequality as above or another insolvable integral.
 
undefined83 said:
Thanks, but i am totaly lost when i need to find int(xn-1,infinity). I try to make it lower than error estimate but i get either logaritmic inverse-trigonometric inequality as above or another insolvable integral.
Don't over think it
|{\int_x^\infty \frac{\sin(\frac{1}{x})}{1+x^3}}|&lt;\int_x^\infty \frac{1}{x^3}dx=\frac{1}{2x^2}
 
Last edited:
Thank you very, very much. You are a life savior.
 
Does someone who has Maple or Mathlab can give me aproximate solution of the integral, posiblly of int(10,45), and int(10,142). Error comarision goes to absurd. I get 3.3*10^-4 on pocket calculator for int(10,35). Simpson with above given boundries gives more than 10 times larger values.
 
undefined83 said:
Does someone who has Maple or Mathlab can give me aproximate solution of the integral, posiblly of int(10,45), and int(10,142). Error comarision goes to absurd. I get 3.3*10^-4 on pocket calculator for int(10,35). Simpson with above given boundries gives more than 10 times larger values.

\int_{10}^{45} \frac{Sin(1/x)}{1+x^3}dx\approx 0.000329176

\int_{10}^{142} \frac{Sin(1/x)}{1+x^3}dx\approx 0.000332717

As per Mathematica's NIntegrate.
 
Last edited:
I don't know if you want to stay with the original integral or not. In any case a substitution will make this much easier.
u=1/x
\int_{10}^\infty \frac{\sin(\frac{1}{x})}{1+x^2}dx=\int_0^\frac{1}{10}\frac{u\sin(u)}{1+u^3}du
then since x(sin(.1)/.1)<sin(x)<x for 0<x<.1
\frac{\sin(.1)}{.3}\log(1.001)=10\sin(1/10)\int_0^\frac{1}{10}\frac{u^2}{1+u^3}du&lt;\int_0^\frac{1}{10}\frac{u\sin(u)}{1+u^3}du&lt;\int_0^\frac{1}{10}\frac{u^2}{1+u^3}du=\frac{1}{3}\log(1.001)
The average of these yeild a good approximation.
simpsons rule will meet the error tolerance
I~(.1-0)/6(f(0)+4f(.05)+f(.1))
where f(x)=x sin(x)/(1+x^3)
 
Last edited:
  • #10
Thank. This afternoon i finaly finished the monster with int(10,150) by Simpson. Got aprox 0.000324. I can't really believe its finished.
 

Similar threads

Back
Top