Once Hk_comp.m has been run, we can solve the boundary integral equation for the traces of the
complex geometric optics solutions. This is done by the routine psi_BIE_comp.m, and the result is saved
to a file in the subdirectory 'data'. The next step is to evaluate the scattering transform by integration over
the boundary using the routine tBIE_comp.m.
We are now ready to plot the result and compare it to the scattering transform computed using the
Lippmann-Schwinger equation approach. Run the file t_plot.m, and you should see something like this: