lt;- diag(1, 1 + 2*N) --&gt;</span>
<span class="co"># &lt;!-- L2[cbind(1+(1:N), 1:N)] &lt;- -alpha --&gt;</span>
<span class="co"># &lt;!-- L2[cbind(1 + N + (1:N), 1 + 1:N)] &lt;- -beta --&gt;</span>
<span class="co"># &lt;!-- K2 &lt;- L2 %*% t(L2) --&gt;</span>
<span class="co"># &lt;!-- V2 &lt;- solve(K2) --&gt;</span>
<span class="co"># &lt;!-- sparsify(K2) --&gt;</span>
<span class="co"># &lt;!-- sparsify(V2) --&gt;</span>
<span class="co"># &lt;!-- ``` --&gt;</span>
<span class="co"># </span>
<span cl