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)
とすべきであった。



桂田 祐史