;=1</span>
0017 <span class="comment">%       phierr - standard deviation for phi (note that the routine gives phierr as phierr(1,...)</span>
0018 <span class="comment">%                and phierr(2,...) in order to incorporate Jackknife (eventually).</span>
0019 <span class="comment">%                Currently phierr(1,...)=phierr(2,...). Note that phi + 2 phierr(1,...) and phi -2</span>
0020 <span class="comment">%                phierr(2,...) will give 95% confidence bands for phi - only for err(1)&g