Dear Folks, I find this interesting. I tried this integral using the CauchyPrincipalValue command and got one result. I also tried doing the principal value manually (integrating up to the pole then from the pole to infinity and adding the two) and ended up with something different. Any thoughts? Could somebody try it on Mathematica (I don't have it on my computer)? Int(omega*sqrt(x)/(x*(omega-x)),x=0..infinity)=int(omega*sqrt(x)/(x*(omega-x)),x=0..infinity,CauchyPrincipalValue); int(sqrt(x)*omega/(x*(omega-x)), x = 0 .. omega-epsilon, assume(omega > 0, epsilon > 0))+int(sqrt(x)*omega/(x*(omega-x)), x = omega+epsilon .. infinity, assume(omega > 0, epsilon > 0)); Drew

Please Wait...