@Michael Freimann
Lo que sigue no es una respuesta. Es sólo un comentario aparte, que podrían ser de interés para algunos lectores.
El tipo de matrices en cuenta en esta pregunta puede ser descrito por la siguiente fórmula para sus coeficientes
$$A_{ij}=a_{min(i,j)}$$ for a given sequence $a_1,a_2,...a_n$, esta secuencia de estar aquí {a,b,c,d}.
Sus inversos tienen de forma tridiagonal con simple de los coeficientes. Por ejemplo, aquí:
$A^{-1}=\begin{pmatrix} \frac{1}{a} + \frac{1}{b-a}&\frac{1}{a - b}&0&0\\
\frac{1}{a - b}&\frac{1}{b-a} + \frac{1}{c-b}&\frac{1}{b - c}&0\\
0&\frac{1}{b - c}&\frac{1}{c-b} + \frac{1}{d-c }&\frac{1}{c - d}\\
0&0&\frac{1}{c - d}&\frac{1}{d-c}\end{pmatrix}$
En particular, cuando se $a=b=c=d=1$,
$$\text{La inversa de} \ \ \begin{pmatrix}
1&1&1&1\\
1&2&2&2\\
1&2&3&3\\
1&2&3&4
\end{pmatrix} \ \ \text{es} \ \
\frac{1}{2}\begin{pmatrix} 2&-1&0&0\\
-1&2&-1&0\\
0&-1&2&-1\\
0&0&-1&1\end{pmatrix}
$$
Uno puede reconocer en la última matriz, una aproximación al contrario de la segunda derivada, una parte muy importante de la matriz de análisis numérico: consulte este (diapositivas 10, 29, 35, 37, 42...).
Esto está en clara conexión con la descomposición dada por @user1551, que, por $a=c=d=1$, se tiene:
$$A=\pmatrix{1&0&0&0\\ 1&1&0&0\\ 1&1&1&0\\ 1&1&1&1}
\pmatrix{1&1&1&1\\ 0&1&1&1\\ 0&0&1&1\\ 0&0&0&1}$$
que es equivalente a un doble discreta integración de operador.
Por supuesto, esta propiedad se extiende de forma natural a cualquier dimensión. Por otra parte, tiene una buena correspondencia con la "continua mundo": consulte este.