27.12.2012 Views

Wave Propagation in Linear Media | re-examined

Wave Propagation in Linear Media | re-examined

Wave Propagation in Linear Media | re-examined

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

7.5 Test of the quadratu<strong>re</strong> rout<strong>in</strong>e<br />

Example 7.5.3<br />

Z 1<br />

In[20]:= f[x_] := S<strong>in</strong>[x]/Sqrt[1+x];<br />

SetP<strong>re</strong>cision[OscInt[f,#&,0],16]<br />

Out[20]= 0.8095254817472858<br />

Example 7.5.4<br />

0<br />

Z 1<br />

1<br />

s<strong>in</strong> x<br />

p 1+x dx =0:809525481747<br />

s<strong>in</strong>(x + 1<br />

x )<br />

p dx =0:232948197094<br />

x<br />

In[21]:= f[x_] := S<strong>in</strong>[x+1/x]/Sqrt[x];<br />

SetP<strong>re</strong>cision[OscInt[f,(#+1/#)&,1],16]<br />

Out[21]= 0.23294819709403225<br />

Example 7.5.5<br />

Z 1<br />

1<br />

s<strong>in</strong>(x + 1 p ) x<br />

p dx =0:0416328516893<br />

x<br />

In[22]:= f[x_] := S<strong>in</strong>[x+1/Sqrt[x]]/Sqrt[x];<br />

SetP<strong>re</strong>cision[OscInt[f,(#+1/Sqrt[#])&,1],16]<br />

Out[22]= 0.04163285168937769<br />

In the last two examples it would have also been possible to use only the l<strong>in</strong>ear part of the<br />

argument function for the determ<strong>in</strong>ation of the partition po<strong>in</strong>ts by substitut<strong>in</strong>g the pu<strong>re</strong><br />

function #& for the <strong>re</strong>spective arguments <strong>in</strong> the function call. The <strong>re</strong>sults would <strong>in</strong> fact have<br />

been the same.<br />

Some other functions with l<strong>in</strong>early grow<strong>in</strong>g numbers of oscillations have been exam<strong>in</strong>ed by<br />

Hasegawa and Torii [127].<br />

Example 7.5.6<br />

Z 1<br />

0<br />

e ,x cos xdx = 1<br />

2<br />

183

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!