2d acoustics with a radially-symmetric solution.
The data is a smooth radial hump specified in qinit.f
The 1d equation with source term is solved in subdirectory 1drad
to create an accurate reference solution.

Use compare.m to plot a scatter plot of the 2d solution along with the
reference solution.

error2rad.m also computes and plots the error.