ortho_basis

Routines that facilitate integration by gaussian quadrature through various orthonormal basis

Functions

void gauleg(double x1, double x2, double x[], double w[], int n)

Routine from Numerical Recipes that computes the abscissas and weights for Gauss-Legendre polynomials.

Modified version of Numerical Recipes in C: changed from fortran style indexing to c style starting at index 0 and going till index N-1 for length N