Comment creer une fonction integrale dependant dun parametre

Bonjour
Je voudrais creer une fonction mathematique qui est une integrale dependant d'un parametre.
Par exemple la fonction gamma:
f:=x->int(t^(x-1)*exp(-t)dt,t=0..infini)
est un parfait exemple de ce que je veux faire

quel est selon vous meilleur moyen pour coder une telle fonction dans en C++ (classes, template.....)
Je ne veux pas recourir a des solutions detournee, comme par exemple l'utilisation de series

Je compte calculer l'integrale en utilisant une methode du type quadpacK (GNU GSL).

merci