Do these integrals really diverge? [closed]

I have the following function which I want to Integrate from zero to infinity:

f[n_] := (-(-1 + Sqrt[1 + y^(2 (n - 2))])^(
1/3) + (1 + Sqrt[1 + y^(2 (n - 2))])^(1/3))^4/ (y^((n - 2)/3))

for $$n=5$$ the integral converges

NIntegrate[f, {y, 0, Infinity}]
5.31289

from literature we expect that it also should converge for $$n>4$$, i.e., $$n=5,6,7,...$$. But Mathematica compute the integral for $$n=5,6,7$$ with errors and it seems that the integral diverges for $$n=6,7$$. My question is that should I trust this outcome? Or maybe there is a way that these errors can be fixed and the integrals become finite?

NIntegrate[f, {y, 0, Infinity}, MaxRecursion -> 100]

During evaluation of In:= NIntegrate::slwcon: Numerical integration converging too slowly; suspect one of the following: singularity, value of the integration is 0, highly oscillatory integrand, or WorkingPrecision too small.

During evaluation of In:= NIntegrate::zeroregion: Integration region {{0.,0}} cannot be further subdivided at the specified working precision. NIntegrate assumes zero integral there and on any further indivisible regions.

During evaluation of In:= NIntegrate::eincr: The global error of the strategy GlobalAdaptive has increased more than 400 times. The global error is expected to decrease monotonically after a number of integrand evaluations. Suspect one of the following: the working precision is insufficient for the specified precision goal; the integrand is highly oscillatory or it is not a (piecewise) smooth function; or the true value of the integral is 0. Increasing the value of the GlobalAdaptive option MaxErrorIncreases might lead to a convergent numerical integration. NIntegrate obtained 442.0378757272312 and 0.0007603109345110766 for the integral and error estimates.

442.038

NIntegrate[f, {y, 0, Infinity}, MaxRecursion -> 100]
3.19445*10^26

NIntegrate[f, {y, 0, Infinity}, MaxRecursion -> 100]

During evaluation of In:= NIntegrate::zeroregion: Integration region {{0.,0}} cannot be further subdivided at the specified working precision. NIntegrate assumes zero integral there and on any further indivisible regions.

During evaluation of In:= NIntegrate::slwcon: Numerical integration converging too slowly; suspect one of the following: singularity, value of the integration is 0, highly oscillatory integrand, or WorkingPrecision too small.

6.72576*10^51
• Series[f, {y, 0, 1}, Assumptions -> y > 0] implies divergence at y == 0. As does Integrate[f, {y, 0, Infinity}], which basically does the same computation. Nov 14 '21 at 17:15
• Is the numerator supposed to vanish as $y\rightarrow0$? It doesn't as coded. Nov 14 '21 at 17:19
• @MichaelE2 Thank you for your comment, but this is also the case for $f$, but integrating $f$ leads to a finite number. Nov 14 '21 at 17:20
• I get that the integral of f converges, $\sim 16^{1/3} y^{-2/3}$ in which the power is $> -1$. Nov 14 '21 at 17:21
• @MichaelE2 you are right. Thank you very much. Nov 14 '21 at 18:01 