double ***d3tensor(long nrl,long nrh,long ncl,long nch,long ndl,long ndh)Return value:
- pointer to three dimensional array of doubles
Input:
- nrl, nrh: lowest and highest values of indices for first index (rows)
- ncl, nch: lowest and highest values of indices for second index (columns)
- ndl, ndh: lowest and highest values of indices for third index (depth)
void free_d3tensor(double ***t, long nrl,long nrh,long ncl,long nch,long ndl,long ndh)Input:
- nrl, nrh: lowest and highest values of indices for first index (rows)
- ncl, nch: lowest and highest values of indices for second index (columns)
- ndl, ndh: lowest and highest values of indices for third index (depth)
To store the values for each wave function 
, we
need a three-dimensional, double-precision array (with indices n,
, k). Numerical recipes has routines to handle such objects for
normal floating point, but not double precision. Generalize the
Numerical Recipes utility routines ``f3tensor'' and ``free_f3tensor''
to double precision by replacing all occurrences of float with
double, thereby producing routines with the two prototypes above.