Solution
= Solution
The proposed stability property is false. Set $u^*=0$, $f^*=0$, and take $u_n=e_n$. Then
$$
f_n=Du_n=\frac1n e_n,
\qquad
\|f_n-f^*\|_2=\frac1n\longrightarrow0,
$$
while $\|u_n-u^*\|_2=1$ for every $n$. Hence the inverse is not continuous on its range and the recovery problem is not <Hadamard well-posedness>[stable with respect to perturbations]. This is the standard <unbounded inverse on a nonclosed operator range>.