Hello, > >>> from scipy.special import erf > >>> erf(2+3j) > (-20.829461428577808+8.6873182710560588j) Yes, thanks. Unfortunately, the library used for this is Fortran, so I'm still open for other suggestions. > In C, I'd look into libgsl0 which provides complex number routines and > special functions. The special functions in GSL are mostly double-valued and take only doubles as input. Regards Thomas