clear
close

load om_dat;
omega=om_dat(1);

load scrc_dat;
size(scrc_dat)
load scrr_dat;
load scri_dat;
sca=sqrt(scrr_dat.^2 + scri_dat.^2);
load phsr_dat;
m=size(phsr_dat,1);
for i=1:m,
  a(m)=0;
end

cmin=min(scrc_dat)*0.95;
cmax=max(scrc_dat)*1.05;

xmin=min(scrr_dat);
xmin1=min(scri_dat);
if xmin1 < xmin 
  xmin=xmin1;
end
xmin=xmin*1.02;
xmax=max(sca)*1.01

%cmin=600
%cmax=700

xmax=2
xmin=-xmax

plot(scrc_dat,sca,'-',scrc_dat,scrr_dat,'--',scrc_dat,scri_dat,'-.')
hold
plot(phsr_dat,a,'o')
grid
axis([cmin cmax xmin xmax])
xlabel('Velocity (m/sec)')
ylabel(['Secular Function of Rayleigh Wave at ',num2str(omega),' (rad/sec)'])
title('Dashed: Real Parts;  Dash-Dotted: Imaginary Parts;  Solid: Absolute Values')
explanation='Circles in the figure indicate the locations of phase velocities'

print secr.ps




