dy/dx=HankelH1[0, x], y(0)=pi/2