The question: Let $f$ be a Lebesgue integrable function on $\mathbb{R}$ and let $b\in (0,1)$. Prove that $\int _{0} ^{\infty} \frac{|f(x)|}{|x-a|^b} dx <\infty$ for a.e. $a\in \mathbb{R}$.
My progress: I first consider $E_n := \{x\in [0,\infty): \frac{1}{|x-a|^b} \geq n \}$. Then it is quite obvious that $m(E_n) \leq 2\cdot (\frac{1}{n})^{\frac{1}{b}}$.
Then $\int _{0} ^{\infty} \frac{|f(x)|}{|x-a|^b} dx = \int _{E_n\cap [0,\infty)} \frac{|f(x)|}{|x-a|^b} dx+ \int _{E_n^C\cap [0,\infty)}\frac{|f(x)|}{|x-a|^b} dx$. Here is where I am unsure what to do. $\int _{E_n^C\cap [0,\infty)}\frac{|f(x)|}{|x-a|^b} dx$ is finite for $n<\infty$, but how to deal with $\int _{E_n\cap [0,\infty)} \frac{|f(x)|}{|x-a|^b} dx$ is an issue. Any comments are appreciated!
Let $N$ be a positive integer.
By Tonelli's Theorem we have $\int_{-N} ^{N} \int_{-\infty} ^{\infty}\frac {|f(x)|} {|x-y|^{b}} dxdy=\int_{-\infty} ^{\infty}\int_{-N} ^{N}\frac {|f(x)|} {|x-y|^{b}} dydx$. Now $\int_{-N} ^{N}\frac 1 {|x-y|^{b}} dy=\int_{x-N} ^{x+N}\frac 1 {|z|^{b}} dz$. I will let you check that $x \to \int_{x-N} ^{x+N}\frac 1 {|z|^{b}} dz$ is bounded on the real line for $0<b<1$. Conclude that $\int_{-\infty} ^{\infty}\frac {|f(x)|} {|x-y|^{b}} dx<\infty$ for almost every $y \in [-N,N]$. This is true for each $N$.
[Consider, for example, the limit of $(x+N)^{1-b} -(x-N)^{1-b}$ as $ x \to \infty$. We can write $$(x+N)^{1-b} -(x-N)^{1-b}=x^{1-b}[(1+\frac N x)^{1-b}-(1-\frac N x)^{1-b}]= x^{1-b}[(1+\frac {(1-b)N} x)-(1-\frac {(1-b)N} x)+o(\frac 1 x)]\to 0$$
If a continuous function on $\mathbb R$ has finite limits at $\pm \infty$ then it is bounded].