Index: trunk/doc/pslib/psLibADD.tex
===================================================================
--- trunk/doc/pslib/psLibADD.tex	(revision 756)
+++ trunk/doc/pslib/psLibADD.tex	(revision 758)
@@ -1,3 +1,3 @@
-%%% $Id: psLibADD.tex,v 1.16 2004-05-24 19:27:53 eugene Exp $
+%%% $Id: psLibADD.tex,v 1.17 2004-05-24 20:14:06 eugene Exp $
 \documentclass[panstarrs]{panstarrs}
 
@@ -42,6 +42,6 @@
 SLALIB Positional Astronomy Library & http://star-www.rl.ac.uk/star/docs/sun67.htx/sun67.html \\ \hline
 Numerical Recipes (NR)              & \\ \hline
-Knuth                               & \\ \hline
-Sedgewick                           & \\ \hline
+Knuth, D.E.                         & Sorting and Searching; The Art of Computer Programming \\ \hline
+Sedgewick, R.                       & Algorithms, Ch. 8 \\ \hline
 Sorting Summary                     & {\tt http://www.iti.fh-flensburg.de/lang/algorithmen/sortieren/algoen.htm } \\ \hline
 GSL                                 & \\ \hline
@@ -231,8 +231,7 @@
 valued, the natural bin size is an integer.  Otherwise, the bin should
 be a fraction of an estimate of the standard deviation.  Use the
-sample upper and lower quartiles to determine an estimate of the
-standard deviation: $\sigma_e = (U_{\frac{1}{4}} - L_{\frac{1}{4}}) /
-1.34$.  The bin size shall be set at $\sigma_e / 10$.  The remaining
-steps of the algorithm are as follows:
+$3\sigma$ clipped standard deviation as an estimator of the standard
+deviation.  The bin size shall be set at $\sigma_e / 10$.  The
+remaining steps of the algorithm are as follows:
 
 \begin{itemize}
@@ -242,5 +241,5 @@
 \item Find the bin with the peak value in the range $L_{\frac{1}{4}}$
   to $U_{\frac{1}{4}}$; this is the robust mode, $\mbox{mode}_r$.  
-\item Determine $dL = (U_{\frac{1}{4}} - L_{\frac{1}{4}}) / 8$.
+\item Determine $dL = (U_{\frac{1}{4}} - L_{\frac{1}{4}}) / 4$.
 \item Fit a Gaussian to the bins in the range $\mbox{mode}_r - dL$ to
   $\mbox{mode}_r + dL$.
@@ -253,5 +252,7 @@
 percentile value and its two neighbors.  Fit a quadratic to these
 three points.  The robust median value is the coordinate of the
-quadratic which returns the 50\% value.
+quadratic which returns the 50\% value.  For the upper and lower
+quartile points, the same process should be used, choosing the three
+bins in the vicinity of the upper and lower quartile points.