compare.m.html
Source file:   compare.m
Directory:   /home/rjl/git/rjleveque/clawpack-4.x/apps/acoustics/2d/example1
Converted:   Tue Jul 26 2011 at 12:58:47   using clawcode2html
This documentation file will not reflect any later changes in the source file.

 

% Scatter plot of the 2d solution compared with 1d radially symmetric solution
%
% First run code in this directory and in 1drad with same set of output times.  
%
% Use error2rad to also compute and plot the error for a particular Frame

for Frame=0:30
   PlotType=4;
   plotframe2
   qscat = q;
   rscat = r;
   cd 1drad
   plotstyle = 'r-';
   plotframe1
   x1d = x;
   q1d = q;
   axis([0 1.5 -.5 5])
   cd ..
   hold on
   plot(rscat,qscat','or')
   plot(x1d,q1d,'k')
   hold off
   query
   end