Analytical Solution of 2d acoustic pulse propagation
 

Hi all,

for an 2d aeroacoustic benchmark test case I had a numerical solution already. But now I want to compare the results with the analytical/exact solution.
The analytical solution is presented by Tam, Christopher KW, and Jay C. Webb. "Dispersion-relation-preserving finite difference schemes for computational acoustics." Journal of computational physics 107.2 (1993): 262-281.
My questions is, how do I manage it to plot the analytical values in matlab/octave/excel/ for comparision with my numerical values? I want to plot the results for a specific time step along the x-Axis. I have tried to solve the integral of the analytical solution in mahtematica but I was not successfull.
Hope anyone might have a helping hint for me. Thanks a lot in advance,
Jan