Index: trunk/psLib/src/math/psPolynomial.c
===================================================================
--- trunk/psLib/src/math/psPolynomial.c	(revision 793)
+++ trunk/psLib/src/math/psPolynomial.c	(revision 822)
@@ -1,2 +1,17 @@
+/** @file  psFunctions.c
+ *
+ *  @brief Contains basic function allocation, deallocation, and evaluation
+ *         routines.
+ *
+ *  @author George Gusciora, MHPCC
+ *
+ *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-06-01 22:24:55 $
+ *
+ *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ */
+/*****************************************************************************/
+/*  INCLUDE FILES                                                            */
+/*****************************************************************************/
 #include <stdlib.h>
 #include <stdio.h>
@@ -18,4 +33,8 @@
 #include <gsl/gsl_rng.h>
 #include <gsl/gsl_randist.h>
+/*****************************************************************************/
+/*  FUNCTION IMPLEMENTATION - PUBLIC                                         */
+/*****************************************************************************/
+
 /*****************************************************************************
     Evaluate a non-normalized Gaussian with the given mean and sigma at the
@@ -41,8 +60,7 @@
 /*****************************************************************************
     psGaussianDev()
- 
  This routine creates a psVector of the specified size and type F32
- and fils it with a random Gaussiang distribution of numbers with
- the specified meand and sigma.  This routine makes use of the GSL
+ and fills it with a random Gaussian distribution of numbers with
+ the specified mean and sigma.  This routine makes use of the GSL
  routines for generating both uniformly distributed numbers and the
  Gaussian distribution as well.
@@ -65,4 +83,6 @@
         gauss->vec.f[i] = mean + gsl_ran_gaussian(r, sigma);
     }
+
+    // GUS: Should I free r as well?
     return(gauss);
 }