Vamos
I \subset U \subset \Bbb R^{n \times n}, \tag 1
donde U está abierto. Entonces podemos definir la matriz de la función
F(B) = B^2:U \to \Bbb R^{n \times n}, \tag 2
y nos cuenta que
F(I) = I^2 = I, \tag 3
que es, I está en el rango de F(B).
F(B) = B^2 es de hecho diferenciable; tenemos, por \Delta \in \Bbb R^{n \times n},
F(B + \Delta) = (B + \Delta)^2 = (B + \Delta)(B + \Delta) = B^2 + B\Delta + \Delta B + \Delta^2, \tag 4
F(B + \Delta) - F(B) - (B \Delta + \Delta B) = (B + \Delta)^2 - B^2 - (B\Delta + \Delta B) = \Delta^2; \tag 5
\Vert F(B + \Delta) - F(B) - (B \Delta + \Delta B) \Vert = \Vert (B + \Delta)^2 - B^2 - (B\Delta + \Delta B) \Vert = \Vert \Delta^2 \Vert \le \Vert \Delta \Vert^2; \tag 6
desde
\dfrac{\Vert \Delta \Vert^2}{\Vert \Delta \Vert} = \Vert \Delta \Vert \to 0 \; \text{as} \; \Delta \to 0, \tag 7
vemos que F(B) es diferenciable y su derivada es la lineal mapa
DF(B)(\Delta) = B\Delta + \Delta B; \tag 8
así
DF(I)(\Delta) = I \Delta + \Delta I = 2 \Delta, \tag 9
que claramente es no singular. De ello se sigue del teorema de la función inversa que hay algunos vecindario V de I y una función
R:V \to U \subset R^{n \times n}, \: R(I) = I, \tag{10}
que, por A \in V
F(R(A)) = (R(A))^2 = A; \tag{11}
si ahora nos elija \delta > 0 suficientemente pequeño podemos asegurar que el conjunto de
B(I, \delta) = \{A \in \Bbb R^{n \times n}, \; \Vert I - A \Vert < \delta \} \subset V, \tag{12}
y así, por A \in B(I, \delta) podemos establecer
B = R(A), \tag{13}
y tenemos
B^2 = (R(A))^2 = F(R(A)) = A, \tag{14}
como se desee. OE\Delta