F..2 配列の代入 b=a はコピーでなく参照を作るだけ

これは事前に読んで知っていたつもりなのだが、つい

  au = al = - theta*lambda*ones(N-1)

というコードを書いて、バグになかなか気づかなかったという話である。

(al を書き換えると、au も変わってしまって、 結果が無茶苦茶になる。)

  al = - theta*lambda*ones(N-1)
  au = - theta*lambda*ones(N-1)

とするか

  al = - theta*lambda*ones(N-1)
  au = copy(al)

とすべきであった。