Re: How to work out Real part through Kramers-Kronig relation in mathematica?
- To: mathgroup at smc.vnet.net
- Subject: [mg64058] Re: How to work out Real part through Kramers-Kronig relation in mathematica?
- From: "drnedflanders" <hepcj at hotmail.com>
- Date: Sun, 29 Jan 2006 23:10:13 -0500 (EST)
- References: <dri7bs$9sk$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
I use numerical integration and the package CauchyPrincipalValue.
Try
<< NumericalMath`CauchyPrincipalValue`
CauchyPrincipleValue[(x/(x^2-y^2))*eps2(x),{x,0,{b},Infinity}]
I am not sure how your eps2(x) function is defined but place the
point(s) where the singularity(ies) occur where I have indicated {b}.
and give it a go.
Regards,
Carl.
manuelquintana at infosel.net.mx wrote:
> Try this:
> Integrate[(x/(x^2-y^2))*eps2(x),{x,0,Infinity},PrincipalValue->True]
> where eps2(x) is the imaginary part of the dielectric constant.
>
> Link to the forum page for this post:
> http://www.mathematica-users.org/webMathematica/wiki/wiki.jsp?pageName=Special:Forum_ViewTopic&pid=7931#p7931
> Posted through http://www.mathematica-users.org [[postId=7931]]