[nt1]=length(pAB(:,1));
subplot(111);a=diff(pAB(:,2));amx=max(abs(a));plot([0:nt1-1]*dt,[0; a(1:nt1-1)/amx],'-');
hold on;nt1=length(seismo1(:,sn));plot([0:nt1-1]*dt,seismo1(1:nt1,sn)/max(abs(seismo1(:,sn))),'--');hold off
xlabel('Time (s)');ylabel('Normalized Amplitude');title('G(B,t|A,0)');
text(.05,.6,'Solid  = Analytic');text(.05,.7,'Dashed = Interferometric');