我想得到菲涅尔-基尔霍夫积分在近轴近似的圆柱衍射波的结果，这是我的代码的一部分，如果解决了，将把其他部分带出去。这是一个 Matlab 代码。
vpa(int(exp(-1i*2*k*(n_cl-n_s)*sqrt(b^2-x_p^2))*exp(1i*k*((x-x_p)^2)/2z),-b,-a))
关于 x_p。在哪里
b = 62.5*e-6; % cladding radius
a = 4*e-6; % core radius
n_co = 1.473; % core ref index
n_cl = n_co - 0.016; % cladding ref index
n_s = 1.37868e-5; % surrounding ref index
L = 457.9e-9; % wavelength
k = 2*pi/L; % wavenumber
z = 0.4 ;
x =-0.0002;
哪个命令应该用作积分器？我一直在尝试使用int()，而不是quad因为quad不满足获得复杂结果的要求。除此之外，我还遇到了众所周知的错误
Explicit integral could not be found
不知道该怎么办。请帮我。