    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);