          plot(ax,mn,ones(size(s)),<span class="string">'k'</span>);
            <span class="comment">%ylim([0.8 1.2]);</span>
            hold(ax,<span class="string">'on'</span>);
            plot(ax,mn,s.*s,<span class="string">'o'</span>);
            ylabel(ax,<span class="string">'Variance'</span>);
            xlabel(ax,<span class="string">'Corrected Photon Count'</span>);
            title(ax,[<span class="string">'N = '</span> num2str(m(1)) <span class="string">', Z = '</span> num2str(offset)]);
