For each $0 \leq c \leq \frac{1}{\sqrt{2}}$ define the following region
$$R=\{(x,y): x^2+y^2 \leq 1\} \setminus ([-c,c] \times [-c,c])$$
Compute
$$\lim_{c \rightarrow 0} \iint_{R} \frac{1}{(x^2+y^2)^{3/4}}\,dA$$
Now, if we remove a circle instead of a square then using polar coordinates, we can easily find the answer which is $4\pi$. If we consider R as a type one or two region, then I do not know how to calculate the integral because it seems that we cannot make use of the polar coordinate anymore. Any hint would be appreciated.
Let $I_c$ be the integral of the function on the unit disk with $[-c,c]\times [-c,c]$ square removed. Let $J_c$ be the integral on the unit disk with the radius $c$ smaller disk centered on the origin removed.
As you’ve noticed using polar coordinates, $J_c=\int_c^1 2\pi r/r^{3/2}dr=[4\pi\sqrt{r}]^1_c=4\pi (1-\sqrt{c})$.
Since the function is positive, increasing the covered area will only increase the integral. Since the square in $I_c$ covers a circle of radius $c$ and lies within a circle of radius $2c$, we have that the region of $J_{2c}$ is within that of $I_c$ which is within that of $J_{c}$.
Thus,
$$ 4\pi(1-\sqrt{2c})=J_{2c}\leq I_c\leq J_c = 4\pi(1-\sqrt{c})$$
Applying the Squeeze Theorem gives that $\lim_{c \rightarrow 0^+} I_c=4\pi$ as desired.