pan>
p_xy_rs = interp2(kgrid_recon.y, kgrid_recon.x - min(kgrid_recon.x(:)), p_xy, kgrid.y, kgrid.x - min(kgrid.x(:)));
</pre>

<p>The interpolated reconstructed initial pressure distribution (using a positivity condition) with the same plot scaling as the plot of <code>source.p0</code> above is shown below. A x = const slice through the center of the larger disc is also shown for comparison. The reconstructed pressure magnitude is decreased due to the limited view problem; the reconstruction is only exact 