I came across this problem.
For $X=(x_1,x_2,x_3,x_4)$, define $A(X)=(x_1^2-x_2^2-x_3^2-x_4^2,2x_1x_2,2x_1x_3,2x_1x_4)$. For any $Z\in\mathbb{R}^4$ with $|Z|\leq 1/16$ prove that the equation $X-A(X)=Z$ has a unique solution such that $|X|\leq 1/8$. $|X|$ is the Euclidean norm of $X$.
My thought is to use the inverse function theorem. Let $I$ be the identity on $\mathbb{R}^4$. the derivative of $I-A$ is $\begin{bmatrix}1-2x_1 & -2x_2 & -2x_3 & -2x_4\\ -2x_2 & 1-2x_1 & 0 & 0\\ -2x_3 & 0 & 1-2x_1 & 0\\ -2x_4 & 0 & 0 & 1-2x_1 \end{bmatrix}$. The determinant is
$(1-2x_1)^2((1-2x_1)^2-4x_2^2-4x_3^2-4x_4)$. So this quantity needs to be nonzero for $|X|\leq 1/8$? And if this is true, $I-A$ is bijective from the ball of radius $1/8$ onto its image. So I'd then have to show that the image of this ball is contains the ball of radius $1/16$? Is this the idea?
Write $X=(t,{\bf x})$ and $Z=(s,{\bf z})$ with ${\bf x}, {\bf z}\in{\mathbb R}^3$. Then $$A(X)=(t^2-|{\bf x}|^2, 2t{\bf x})\ ,$$ and the equation $X-A(X)=Z$ expands to $$\bigl(t-t^2+|{\bf x}|^2, (1-2t){\bf x}\bigr)=(s,{\bf z})\ .$$ It follows that necessarily $${\bf x}={1\over1-2t}{\bf z}\ ,$$ which is unproblematic as long as $|t|<{1\over2}$. In addition we now have the following scalar equation for $t$: $$t-t^2+{z^2\over(1-2t)^2}=s\ ,\tag{1}$$ which we have to solve for given $Z$ with $|s|\leq{1\over16}$, $\>z:=|{\bf z}|\leq{1\over16}$. Put $$f(t,s,z):= t-t^2+{z^2\over(1-2t)^2}-s\ .$$ Then $f(0,0,0)=0$ and $$f_t(0,0,0)=1+4z^2\ne0\ .$$ The implicit function theorem now allows the following conclusion: There is a neighborhood $U$ of $(s,z)=(0,0)$ and a neighborhood $V$ of $t=0$ such that for all $(s,z)\in U$ there is exactly one solution $t=\psi(s,z)\in V$ of $(1)$, whereby $\psi$ is a real analytic function of $s$ and $z$.
Of course, for exact estimates concerning the sizes of $U$ and $V$ more work is necessary. The following figure shows the graph of the LHS of $(1)$ for various values of $z$. It shows that the claim made in your question holds true.