Searching for help: Warning: Explicit integral could not be found.

1 view (last 30 days)
Hello,
I'm trying to find an integral for a function, but there is always an Error Message: Warning: Explicit integral could not be found.
The function include two variables:
syms H positive;
syms x;
and looks like this:
f = (40000*(5962725554385875/(9223372036854775808*(40000*((5000*x - 75)^2 + 81/(256*H^2))^(1/2) + 1)) + 5690802099995323/18446744073709551616)*((5000*x - 75)^2 + 81/(256*H^2))^(1/2) - (3162897844476069375*1000^((235758936893218125*((5000*x - 75)^2 + 81/(256*H^2))^(1/2))/1125899906842624 + 10588771210153143/576460752303423488)*((5000*x - 75)^2 + 81/(256*H^2))^(1/2))/549755813888)
I tried:
F = int(f,x,0,B) % with B=0.03 (previously set)
but there is always the Error Message. If anybody have a clue please answer. Btw, the same problem occurs at another function:
g = -(40000*((5000*x - 75)^2 + 81/(256*H^2))^(1/2))/((3330505768082351*(40000*((5000*x - 75)^2 + 81/(256*H^2))^(1/2))^(95251/50000))/9007199254740992 - 811793349164925/274877906944)
G = int(g,x,0,B)

Answers (1)

Walter Roberson
Walter Roberson on 9 Sep 2013
You cannot realistically integrate that f, unless perhaps H is fairly large. If you evaluate f at x=0, H=17, then you will get a value on the order of -6E47122 . The largest f(x) is at x=3/200; with H=17, it is still -1.2E26. The function is symmetric, so by x=3/100 it is back down to -6E47122
I have poked around but I see no reason to expect that there might be a closed form for the symbolic integral.
  2 Comments
Andreas
Andreas on 10 Sep 2013
Thank you, the value of H should be 0.0016, but with this help too, i can't integrate the function.
Walter Roberson
Walter Roberson on 10 Sep 2013
With H = 0.0016, f(0) is about -2 * 10^225826 and f(3/200) is -5 * 10^220856
What sort of range are you expecting that the integral would come out as?

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!