% Emacs mode header: -*-Octave-*- source. %
% Greenwood tuning function
% Gw 


omega=286; %:

plot(-x(ilist(21)),latcum(21),'mo')

omega=23.51; %:

plot(-x(ilist)+2.5,latcum,'go')

-x(ilist(5))

yi = interp1(x,Y,xi)

yi1 = interp1(-x(ilist(7:8))+2.5,latcum(7:8),-x(ilist(5)))

%===
% x/L:  0.25, 0.5, 0.75


% Fig 3A: 40 dB
% --- p0 = 0.02, p20 = 0.0004
% x/L:  0.25, 0.5, 0.75
% omega: 286.5 % q3y(10): p0=0.18, p20=0.03
% omega: 82
% omega: 23.51

% Fig 3B: 80 dB
% --- p0 = 2, p20 = 4
% x/L:  0.25, 0.5, 0.75
% omega: 286.5, gamma = -10
% hlen=show_bc1(286.5,q3y);

%===
% omega: 20,
% q2=logspace(-40,-20,20);

%===
% omega: 82
% q2=logspace(-40,-20,20);


%===
%   omega: 286.5,
%   p0=[1e-9,5e-9,1e-8,5e-8,1e-7,5e-7,1e-6,5e-6,1e-5];
%   q0=logspace(-9,-5,9);
%   q1=logspace(-21,-9,27); 
%   q2=logspace(-40,-20,30);
%   q3=logspace(-400,-20,100);

%   q3y=logspace(-49,-45,20); critic: q(15):q(16)
%>> q3y(15),data_ph(:,15)
% ans =/ 8.8587e-47/ ans =/ 1.0e+03 */ 1.7672/ 0.0060


%===

%plot(x,velo,'r@-o',x, late, 'g--s',x,0.1*pasiv,'b@-+',x,0.1*activ,'k--x')

%===
% Gw10([0,2.5, 5,7.5,10])/(2*pi)
% GinvGw10([16000 4600 1306 374 100]*2*pi)
%===