% Emacs mode header: -*-Octave-*- source. %

%-- lo SPL (40 dB)
% ----- left edge x, low edge y, width x, height y
%axh1=axes('position',[.1  .65  .8  .3]);
%axh2=axes('position',[.1  .35  .8  .3]);
%axh3=axes('position',[.1  .05  .8  .3]);

%figure

cd ../actNL/

slid=0.11; xoff=+35;
xL=linspace(0,35,1000);
%'{\fontsize{14}string}'
hlen=1; % q3y=logspace(-49,-45,20); q3y(10)
[x,P,Q,rh,eh,data_ph] = tw_nonlin_4M3(286.5,7.8476e-48,0.1);
%axh1=subplot(3,1,1);
axh1=axes('position',[0.05  .7  .65  .25]);
p1=plot(xoff-xL,eh(hlen,:),xoff-xL,rh(hlen,:),xoff-xL,-eh(hlen,:));
set(p1(3),'Color','b');
%axis([-10 -4 -0.75 +0.75])
axis([xoff-35 xoff-14 -0.75 +0.75])
set(axh1,'XMinorGrid','on')
%grid
%pause
text(xoff-22,.5,'\fontsize{14}4.6 kHz')
%title
text(xoff-35,.75+0.05,'\fontsize{14}Traveling wave with envelope, low SPL')
%clear all

qx=5.4556e-26; hlen=1; %q2=logspace(-40,-20,20); qx=q2(15)
[x,P,Q,rh,eh,data_ph] = tw_nonlin_4M3(82,qx,0.1);
%axh2=subplot(3,1,2);
axh2=axes('position',[.05+slid  .4  .65  .25]);
p2=plot(xoff-xL,eh(hlen,:),xoff-xL,rh(hlen,:),xoff-x*3.5,-eh(hlen,:));
set(p2(3),'Color','b');
%axis([-8 -2 -0.75 +0.75])
axis([xoff-28 xoff-7 -0.75 +0.75])
set(axh2,'XMinorGrid','on')
text(xoff-15,0.5,'\fontsize{14}1.3 kHz')
%grid

%clear all
qx=1e-09; hlen=1; % q0=logspace(-9,-5,9); qx=q0(1)
[x,P,Q,rh,eh,data_ph] = tw_nonlin_4M3(23.51,qx,0.1);
%axh3=subplot(3,1,3);
axh3=axes('position',[.05+slid+slid  .1  .65  .25]);
p3=plot(xoff-xL,eh(hlen,:),xoff-xL,rh(hlen,:),xoff-xL,-eh(hlen,:));
set(p3(3),'Color','b');
%axis([-6 0 -7.5 +7.5])
axis([xoff-21 xoff-0 -7.5 +7.5])
set(axh3,'XMinorGrid','on')
%grid
text(xoff-5,5,'\fontsize{14}370 Hz')
xlabel('\fontsize{14}Distance from stapes, mm')
%clear all

cd ../mkfig/

%-- hi SPL (80 dB)

