[Help-gsl] Complex integration in real plane

[Ivan Ramirez]

Hi,
I am trying to calculate the propagation of a electromagnetic wave, by
de Rayleigh Sommerfeld function.
This is a complex function (the imaginary part is related with the
phase of the wave) of 2 variables (the area of the exit pupil), so is
a double integral and the limits are the exit pupil.

Is it posible to do this with GSL? Until now I do it with Mathematica,
but I don't like black box,

Any susgestion about the way to do this is appreciate
Thanks in advance
Ivan