Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Difficulty in solving an integral in Mathematica

  1. Dec 8, 2015 #1
    Hello,

    I have to solve the following integral in Mathematica but although I don't get any mistakes the value I get is <10^-700 which can't be right!!

    For the following values

    b=0.052868
    c=0.998047
    g=0.00762712
    s=0.902416

    I have to solve the following:

    Fx = Table[NIntegrate[(Exp[-t*x^2/g]/x^3)*(Sin[x] - x*Cos[x])*(b*x*Sin[x]*Cos[s*x] - (c*Sin[x] - x*Cos[x])*Sin[s*x])/((c*Sin[x] - x*Cos[x])^2 + (b^2*x^2*Sin[x]^2)),{x, 0, Infinity}],{t, 13, 100, 1}];

    Could you possibly help me???I guess the problem is the exponential part which is really huge.This problem has been bothering me for a couple of weeks now!

    Thank you guys

    k0st123
     
  2. jcsd
  3. Dec 9, 2015 #2

    DrClaude

    User Avatar

    Staff: Mentor

    It works fine up to t = 42. For t=43, the integral doesn't converge to the prescribed accuracy. You need to tweak the desired accuracy.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Difficulty in solving an integral in Mathematica
Loading...