 Solution

The range of the Volterra integration operator is $\{g\in H^1(0,1):g(0)=0\}$, using the representative of a Sobolev space element that is an absolutely continuous function. Indeed $Ku$ has weak derivative $u$ and zero initial trace; conversely the fundamental theorem of calculus reconstructs such a $g$ from its derivative. This range contains smooth compactly supported functions and is dense in $L^2(0,1)$.

The given step is in $L^2$, and its value at the single midpoint is immaterial. It cannot be the image of an $L^2$ function: such an image is continuous, whereas no continuous representative agrees almost everywhere with zero on the left half and one on the right half. Its distributional derivative is a Dirac delta distribution, not an $L^2$ function. ThusFor a direct approximation, replace the jump by a linear ramp of width $1/n$ centered at $1/2$. Each ramp starts at zero and has an $L^2$ derivative, so lies in the range, while its squared $L^2$ error is $1/(12n)$. Its derivative norm is $\sqrt{n}$, illustrating unstable differentiation.

The supplied SVD gives $K({\sigma }_j{u}_j)={\sigma }_j^2{v}_j$, or more simply $K^{†}({\sigma }_j{v}_j)=u_j$. Since ${\sigma }_j{v}_j\to 0$ but $\Vert u_j\Vert =1$, the Moore–Penrose inverse is discontinuous. Its domain is the dense, nonclosed range above, and there it is the weak derivative.