I have a matrix of the form $A = bI - J$ where
1. $b$ is a large positive constant so that $A$ is positive definite
2. $J_{ij} = 0$ for $i=j$ and follows a power law off-diagnol. In index notation:
$$ A_{ij} = (b\delta_{ij} - \frac{J_0(1-\delta_{ij})}{|i-j|^\alpha}) \text{ where } \alpha \geq 1 $$
Is it possible to calculate $A^{-1}$ analytically? I hope to show that $A^{-1}$ also follows some type of power law as well as is suggested by numerics.
Any suggestions?