非常に面白いですね．私も早速使ってみました．
GNUの数値計算ライブラリgslにある特殊関数（n次ベッセル関数）
を以下のように内部で呼び出すような関数が機能したので，
こりゃ楽じゃと感激しました．なお，MakefileのLIBSに，-lgsl -lgslcblas
を追加する細工が必要でした．

#include <gsl/gsl_sf_bessel.h>

DLLEXPORT struct value gsl_bessel_Jn(int nargs, struct value *arg, void *p)
{
struct value r;
double x;
int n;

RETURN_ERROR_IF_WRONG_NARGS(r, nargs, 2);
RETURN_ERROR_IF_NONNUMERIC(r, arg[0]);
RETURN_ERROR_IF_NONNUMERIC(r, arg[1]);
r.type = CMPLX;
n = IVAL(arg[0]);
x = RVAL(arg[1]);
r.v.cmplx_val.imag = 0.0;
r.v.cmplx_val.real = gsl_sf_bessel_Jn(n, x);

return r;
}

〔ツリー構成〕

【2697】 import 2014/3/4(火)13:03 竹の (852) ≪

┣【2701】 re(1):import 2014/3/6(木)23:15 松岡 (312)

┣【2702】 re(2):import 2014/3/8(土)14:59 竹の (389)

┣【2703】 re(2):import 2014/3/8(土)15:07 松岡 (361)

┣【2704】 re(3):import 2014/3/8(土)19:42 松岡 (321)

┣【2707】 re(4):import 2014/3/11(火)16:21 松岡 (378)

┣【2708】 re(1):import 2014/3/12(水)23:58 松田七美男 (684)

┣【2709】 re(2):import 2014/3/14(金)08:55 松岡 (220)