figure(101); figNo = 101;
hold on;
displot = disAlphaBeta(1:stopind);
x = linspace(min(regplot),max(regplot),100);
lw = 2; % line width
h1 = plot(regplot,displot,'x','MarkerSize',15);
h2 = plot(x,tau*delta+0.*x);
set(h1,'LineWidth',lw);
set(h2,'LineWidth',lw);
hold off;
axis square;
set(gca,'FontSize',20) % seti.pubFontSize = 20;
filePng = sprintf('output/%s_autoparam3.png',dirDatetime);
fileEpsc = sprintf('output/%s_autoparam3.eps',dirDatetime);
fileFig = strrep(filePng,'.png','.fig');
figNostr = sprintf('-f%d',figNo);
print(figNostr,'-dpng',filePng);
print(figNostr,'-depsc',fileEpsc);
savefig(figNo,fileFig);