I am trying to solve the following equation numerically:
Equation to solve for $y$:
$$na(ay-y)^{n}\int_{ay}^{+\infty} \frac{(x-y)^{-n}}{x-ay}dx=b$$
with for example $a=4$, $n=1.25$ and $b=1.6$.
I have adapted an answer to a similar question that I found on StackExchange:
Pdf[x_, y_, a_, n_] = n (x - y)^(-1 - n) ((-1 + a) y)^n
eps[x_, y_, a_] = (x - y)*a/(x - y*a)
Func[x_, y_, a_, n_] = FullSimplify[eps[x, y, a]*Pdf[x, y, a, n]]
f[y_?NumericQ] :=
NIntegrate[Func[x, y, 4, 1.25], {x, 4*y, Infinity}];
FindRoot[f[y] == 1.6, {y, 0}, WorkingPrecision -> 10000]
But I am getting a list of more than 10 error messages and the solution I get does not seem to be correct. My guess (from looking at plots of the integral NIntegrate[Func[x, y, 4, 1.25], {x, 4*y, Infinity}]
which also generate loads of error messages) is that the solution is very small, that's why I am trying to search around 0. But I have no idea how to find it. Can anyone help please?