Wave Propagation in Linear Media | re-examined

8.4 Test of the package
-200 -150 -100 -50
2
1
-1
-2
The examples show that the initial conditions are met. In order to test if the boundary
conditions are satis ed, we sum up the parts of the wave function on either side of the barrier
edge and compute their di erence for a number of points along the time axis. We then plot
the absolute value j (,0;T), (+0;T)j of this di erence. Let us start with the triangular
wave.
Example 8.4.4
In[4]:= x = 0;
w = 0.2;
k = 5;
boundary = Table[ft,PhiInc[x,t,w,k,Shape->Tria]+
PhiRef[x,t,w,k,Shape->Tria]-
(PhiTrans[x,t,w,k,Shape->Tria]+
PhiEvan[x,t,w,k,Shape->Tria])g,ft,0,20g];
ListPlot[boundary/.fx_,y_g:>fx,Abs[y]g,PlotJoined->True,
PlotRange->All];
-7
6. 10
-7
5. 10
-7
4. 10
-7
3. 10
-7
2. 10
-7
1. 10
5 10 15 20
207