d3tensor() and free_d3tensor()

The output of getallPsis() will be a three-dimensional array, one full wave function for each value of (n,l): Psi[l][n][0..N]. For some reason, Numerical Recipes provides routines for three-dimensional arrays in single precision, but not double precision. Accordingly, 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 below. You should place these routines at the top of ``p480.c'' (and their prototypes in ``p480.h''). Your new routines will have the following prototypes:

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)