PGauss:

------------------------------------------------------------------------
   Computes the Gaussian form of the Legendre functions. 
   Will also compute the first derivatives, if requested.

   Because there is no zero indexing in Matlab, the P's are offset
   in the p matrix as follows:

                 n,m
   p(n+1,m+1) = P 

------------------------------------------------------------------------
   Form:
   [p, dP] = PGauss( nMax, mMax, theta )
------------------------------------------------------------------------

   ------
   Inputs
   ------
   nMax               max value of first index
   mMax               max value of second index
   theta              input value (usually an angle in rads)

   -------
   Outputs
   -------                              n,m
   p                  Gauss functions, P
                        n,m
   dP                 ?P    / ?(theta)      

------------------------------------------------------------------------

Children: