Mathematica Definite integral with some unknown variables

Click For Summary
The discussion centers on evaluating an integral in Mathematica involving unknown variables a, b, and c. The integral is expressed as ∫₀¹ dy (y² (1 - b³y³)^(1/2)) / ((1 - a⁴c²y⁴)^(1/2)). Users highlight that the computation may take a long time due to potential singularities in the integrand, which could lead to the integral being improper or non-convergent. To address this, it's suggested to define a new variable d = a⁴c² and impose the assumption d < 1 to help the software handle the integral more effectively. A code snippet is provided, but it fails to return an output, indicating that further refinement or assumptions may be necessary to achieve convergence and improve computation time.
shinobi20
Messages
277
Reaction score
20
TL;DR
I am evaluating a definite integral with unknown variables in Mathematica.
I am trying to evaluate an integral with unknown variables ##a, b, c## in Mathematica, but I am not sure why it takes so long for it to give an output, so I just decided to cancel the running. The integral is given by,

##\int_0^1 dy \frac{ y^2 (1 - b^3 y^3)^{1/2} }{ (1 - a^4 c^2 y^4)^{1/2} }##
 
Physics news on Phys.org
Depending on the values of ##a, b, c## you can expect singularities in the integrand that may make the integral improper or even non-convergent. So, if the problem allows it, you should probably tell the software explicitly about assumptions on ##a, b, c##.

For example, I would begin by defining ##d := a^4 c^2## and passing the assumption that ##d < 1##.
 
Last edited:
  • Like
Likes Dale
S.G. Janssens said:
Depending on the values of a,b,c you can expect singularities in the integrand that may make the integral improper or even non-convergent. So, if the problem allows it, you should probably tell the software explicitly about assumptions on a,b,c.

For example, I would begin by defining d:=a4c2 and passing the assumption that d<1.
This is my code

Code:
d=2;
z=1;
b=a/z;
h=c^2 a^2d;
Integrate[ ( y^d (1 - (b y)^(d+1))^1/2 )/(1 - h^2d y^2d)^1/2 , {y,0,1}, Assumptions -> {y>0, h<1} ]

It does not return anything.
 

Similar threads

  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
3
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K