>>72
境界条件つけないというのはどうすればいいのでしょうか?
境界上の点をpとした時に、点pでの時間発展を計算するには
d^2 A/dt^2=c^2 d^2 A/dx^2
というように、点pでの二回微分(ラプラシアン)を計算する必要があります。
右辺を離散化すると
c^2 (A(p-1)+A(p+1)-2*A(p))
と成りますが、A(p+1) は計算領域の外側になり、値がありません。
もし、強制的に A(p+1)=0 とすれば、それは Dirichlet 条件と同じになり、
右側に進行していた筈の正弦波が反射して、左側に進行してしまいます。
境界上で具体的に、どう計算すればいいでしょうか?