
subplot(2,2,4) %D
%fXD=logspace(1.2,4.2,1000);
%fYD = fXD + 10;
%fYD2 = 0.00001*(fXD - 1000).*(fXD - 1000) + 8;
%hold all; 

fXD=logspace(1.2,4.2,1000);
fYD2 = 0.00001*(fXD - 1000).*(fXD - 1000)+6; % + 8; =quadratic fit
  fYD3 = 10^(-8)*fYD2.*fXD + 0.05*log(fXD) + 0.2; % +0.2+0.05*log(fXD)
% fYD3 = ;

semilogx(fXD, fYD2, 'k',140, 13.333,'o');
semilogx(fXD, fYD2+7.33, 'k',[140 1000], 13.333*ones(1,2),'s');
grid on
xlabel('Sound Frequency [Hz]');
ylabel('ITD [{\mu}s]');
% axis([10^1 2*10^3 0 30])
axis([20^1 10^4 0 30])
