GSL special functions

Hi all,

I've made Octave bindings to a large part of the special functions in the
GNU Scientific library (http://www.gnu.org/software/gsl/). The code can be
found from the Octave-forge CVS (http://octave.sf.net/) under main/gsl.

The bindings are made semi-automatically from a few template files, which
causes the docstrings to be a little too generic sometimes (i.e. the integer
indices of spherical harmonics are named n,m instead of traditional l,m),
otherwise the calling convention is Octave-like: Two output arguments result
in the error estimate being returned in the second value and input arguments
of differing sizes in the multi-argument functions are treated somewhat
intelligently.

The binding code is mostly trivial, but there's still a possibility of bugs
present, perhaps even in the underlying GSL code. Therefore, I'd eventually
like to have some tests to check for correct values in a few well selected
points in the functions. Does anyone know of an accurate source of such
values for a good number of special functions?

Teemu