Computational Physics, Python, Interpolation, Differentiation, Integration.
(-0.4177906356943205+9.86028723231441j) The behavior of functions for complex arguments is a very important subject x=y=np.linspace(-1,1,102) X,Y=np.meshgrid(x,y) plt.pcolor(X,Y,np.real(sp.hankel1(0.2,X+(0+1j)*Y)),cmap='RdBu') Scipy does not define all special functions for complex values: Summary