SciPy

scipy.special.modfresnelp

scipy.special.modfresnelp(x[, out1, out2]) = <ufunc 'modfresnelp'>

(fp,kp)=modfresnelp(x) returns the modified Fresnel integrals F_+(x) and K_+(x) as fp=integral(exp(1j*t*t),t=x..inf) and kp=1/sqrt(pi)*exp(-1j*(x*x+pi/4))*fp