I want to evaluate the numerical integral:
I1[y_] := NIntegrate[PolyGamma[(3*((0.61^2) - x (1 - x)*(y)^2))/(4*0.001) + 1], {x, 0, 1}]
The plot of the function defined by the integral above is:
The result is highly oscillatory for y>2*0.61. This happens because the PolyGamma function has poles if the argument is a negative integer or zero. My idea is to use Principal Value and exclude the points which the argument is negative and integer.
I reckon that the best approach is:
I2[y_] := NIntegrate[PolyGamma[(3*((0.61^2) - x (1 - x)*(y)^2))/(4*0.001) + 1], {x, 0, 1}, Method ->PrincipalValue, Exclusions -> {(3*((0.61^2) - x (1 - x)*(y)^2))/(4*0.001) + 1 == n}]
But I don’t know how to define "n" as a negative integer. I have two questions:
How can I define n as a negative integer and use it to exclude the problematic points?
There is another way to compute the principal value of this integral?
Mathematica Asked by Victor Alencar on December 28, 2020
1 AnswersAnother way.
Using primary definition of PolyGamma
function as a series:
Integrate[(z/(k (k + z)) - 1/z - EulerGamma) /.
z -> Rationalize[(3*((0.61^2) - x (1 - x)*(y)^2))/(4*0.001) + 1,
0], {x, 0, 1}, Assumptions -> y >= 0]
(*-EulerGamma + 1/k - (
4 ArcTan[(50 Sqrt[3] y)/Sqrt[11203 + 40 k - 7500 y^2]])/(
5 Sqrt[3] y Sqrt[11203 + 40 k - 7500 y^2]) - (
4 ArcTan[(50 y)/Sqrt[11203/3 - 2500 y^2]])/(
5 y Sqrt[33609 - 22500 y^2])*)
and we have:
F[y_, M_] := -EulerGamma +
Sum[1/k - (4 ArcTan[(50 Sqrt[3] y)/Sqrt[11203 + 40 k - 7500 y^2]])/(5 Sqrt[3] y
Sqrt[11203 + 40 k - 7500 y^2]), {k, 1, M}]
- (4 ArcTan[(50 y)/Sqrt[11203/3 - 2500 y^2]])/(5 y Sqrt[33609 - 22500 y^2])
Plot[F[y, 2000] // Re, {y, 0, 2}, PlotRange -> {Automatic, {2, 6}}]
Plot[F[y, 5000] // Re, {y, 0, 2},
PlotRange -> {{1.225, 1.5}, {2, 6}}](*In range 1.225 < x < 1.5*)
Answered by Mariusz Iwaniuk on December 28, 2020
1 Asked on January 31, 2021
1 Asked on January 30, 2021 by markokram
1 Asked on January 30, 2021 by john-hall
0 Asked on January 30, 2021 by anon248
2 Asked on January 30, 2021
1 Asked on January 30, 2021 by jakubl
1 Asked on January 29, 2021 by chikitin
9 Asked on January 29, 2021 by nasser
1 Asked on January 29, 2021
1 Asked on January 29, 2021
0 Asked on January 28, 2021
differential equations equation solving plotting warning messages
0 Asked on January 28, 2021 by kirma
0 Asked on January 28, 2021
functions optional arguments options output formatting simplifying expressions
0 Asked on January 28, 2021
1 Asked on January 28, 2021 by shrocat
4 Asked on January 28, 2021 by literal
1 Asked on January 27, 2021 by user34018
Get help from others!
Recent Answers
Recent Questions
© 2022 AnswerBun.com. All rights reserved.