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;aaa=aa+b1;plot([0:nt1-1]*dt,aaa(1:nt1)/max(abs(aaa(:))),'--');hold off
xlabel('Time (s)');ylabel('Normalized Amplitude');title('G(B,t |A,0)','fontsize',14);
text(.3,.6,'Solid  = Analytic');text(.3,.7,'Dashed = Interferometric');axis([0 .7 -1 .88])
text(-.03,.95,'b).','fontsize',14)
print -dpsc GAB.ps
