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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

8.4 Test of the package<br />

Out[24]= 0<br />

Example 8.4.11 Outside the barrier we must comb<strong>in</strong>e the <strong>in</strong>cident and <strong>re</strong> ected waves to obta<strong>in</strong><br />

the desi<strong>re</strong>d <strong>re</strong>sult. He<strong>re</strong> the potential vanishes of course.<br />

In[25]:= Schroed<strong>in</strong>ger[0][TriaIncNeg[-t,b,x,d,w,k][xi] +<br />

TriaRefNeg[-t,b,-x,d,w,k][xi]] //Together<br />

Out[25]= 0<br />

Perform<strong>in</strong>g these tests on all di e<strong>re</strong>nt functions, we see that they satisfy the <strong>re</strong>qui<strong>re</strong>ments.<br />

However, with this formal veri cation we can only ensu<strong>re</strong> the cor<strong>re</strong>ct implementation of the<br />

<strong>in</strong>tegrand functions. We cannot check if those parts of the <strong>in</strong>tegrals that comprise several<br />

dist<strong>in</strong>ct terms (i. e. all transmitted portions of the waves) a<strong>re</strong> put together without error. Nor<br />

can we prove the cor<strong>re</strong>ctness of the quadratu<strong>re</strong> modules di<strong>re</strong>ctly. This must be done with<br />

numerical experiments like the ones above. But tak<strong>in</strong>g together the <strong>re</strong>sults of both formal<br />

and numerical veri cation, we conclude that the implementation is right.<br />

213

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

Saved successfully!

Ooh no, something went wrong!