/N,h);
    <span class="keyword">end</span>
    close(h);

    st = [<span class="string">'Counts per photon = '</span> num2str(mean(counts_per_photon),<span class="string">'%.2f'</span>) <span class="string">' +/- '</span> num2str(std(counts_per_photon),2)];
    disp(st);
    disp([<span class="string">'           Offset = '</span> num2str(mean(offset),<span class="string">'%.2f'</span>) <span class="string">' +/- '</span> num2str(std(offset),2)])

    hold(ax,<span class="string">'off'</span>)
    plot(ax