Solution (source code)

= Solution

First regress $X$ on $Z$, then regress $Y$ on the fitted values. With one centred instrument this equals $S_{ZY}/S_{ZX}$. Its <limit> is $\operatorname{Cov}(Z,Y)/\operatorname{Cov}(Z,X)=1/1=1$, the causal coefficient of $X$ in $Y=X+U+\epsilon_Y$.