Wave Propagation in Linear Media | re-examined

8.4 Test of the package
0.4
0.2
-0.2
-0.4
5 10 15 20
As for the di erent implementations of the evaluation of the inde nite integrals, we compare
the results they yield as well as their computing time. Not surprisingly, we nd that it takes
longer to use the zeros as partition points and that for Gaussian waves the truncation of the
integrals is much faster.
Example 8.4.7
In[7]:= x = 1;
t = 4;
w = 0.1;
k = 2;
Timing[N[Phi[x,t,w,k,Shape->Tria]]]
Out[7]= f21.696 Second, 0.0717104 - 0.0469158 Ig
In[8]:= Timing[N[Phi[x,t,w,k,Shape->Tria,PPoints->Zeros]]]
Out[8]= f29.056 Second, 0.0717104 - 0.0469158 Ig
In[9]:= Timing[N[Phi[x,t,w,k,Shape->Rect]]]
Out[9]= f18.785 Second, 0.0188638 - 0.264976 Ig
In[10]:= Timing[N[Phi[x,t,w,k,Shape->Rect,PPoints->Zeros]]]
Out[10]= f22.903 Second, 0.0188638 - 0.264976 Ig
In[11]:= n = 6;
Timing[N[Phi[x,t,w,k,n,Shape->Gauss]]]
209