#1




Problem with integration.
Hi, already tried few routines from GSL and it seems it doesn't work.
Function: 1/(x^2  alpha^2) Can anyone tell how to calculate that numerically. Tried to do it by "hand" as well but no good results. Cheers. 
#2




Huh? What doesn't work?
What, exactly, did you try? (What functions?) What, exactly, did you feed the functions? (Integration limits? Value of alpha?) What, exactly, did you expect to get? What, exactly, did the rascally functions give you? What, exactly, do you not understand about the difference between what you expected and what you got? Quote:
What, exactly, is your assignment? Is it to calculate a particular value of integral (with given limits and alpha)? Or is it to derive (or discoverby reading a calculus book) an analytical formula for the integration and apply it to a particular problem? Or is it to discover which numerical methods might be appropriate (and successful) for a given problem by comparing numerical results with the analytical solution? Or what? Regards, Dave 
#3




Hi,
actually the function looks like this [tex] f(y,z) = \int^{\infty}_{\infty} \left( \frac{2y}{x^2y^2} \right) \frac{dx}{e^{xz}+1} [tex] so you have always two singularities at +/y. Already tried to do simple sum getting symmetrically from both sides of singularities  but results are bad. Tried GSL integration: qags, qagp, qawc.  answer I got: bad integral behavior or integral is divergent, or slowly convergent Cheers 
#4




That's a lot different from your first post.
The first part of your original question concerns numerical solutions, but how about looking at the function and seeing what to expect? I mean, even if your goal is an analytical solution, sometimes a numerical approach can give some insight, but you still have to think about things. I see the product of two terms under the integral sign: 2y/(x^2y^2) and 1/(exp(xz) + 1) Let's just pick convenient values of y and z and see what it looks like. For example, let y = 1 and z = 0. You want to integrate from x = minus infinity to infinity. Then the two product terms of the integrand are 2/(x^2  1) and 1/(exp(x) + 1) Now, without trying to come up with an analytic solution, you could note that the second term is pretty civilized (stays finite everywhere; approaches zero at one limit, approaches one at the other I think), but the first term has a nonintegrable singularities at x = 1 and x = 1. (The function in your first post has nonintegrable singularities at x = alpha and x = alpha, right?) What happens with your numeric attempts if you try integrate this specific product (with y = 1, z = 0) as you approach x = 1 from the left? Try it for, say, 0 to 0.9? How about 0 to 0.99? Etc. Is it just possible that the function is not integrable over the interval minus infinity to infinity for these particular values of y and z? Are there any other values of y and z that make it behave differently (and, maybe, allow some additional insight)? Or...maybe I'm missing something. (It wouldn't be the first time.) Regards, Dave 
Thread Tools  
Display Modes  

