Nsfdtd

##概要
Ns-FDTD のPML版の実装. 

ただし, 完全なNS版はまだうまくいっていないので
PML領域はS-FDTD, 計算領域はNS-FDTD法で計算するようにしている.

実装時の注意としては
波数k は 屈折率nに依存するため, 岡田さんの修論 Uns を計算するときは, k*nを使うようにする必要がある.
(角周波数ωは必要無い)