Matlab 近轴近似下的菲涅耳-基尔霍夫积分

我想用数值积分这个方程

b = 62.5*10^-6;
a = 4*10^-6;
n_co = 1.473;
n_cl = n_co - 0.016;
n_s  = 1.37868*10^-5;
L = 457.9*10^-9;
k = 2*%pi/L;
z = 0.4 ;
x =0.3;

int(exp(-1i*(2*k*(n_cl-n_s)*(sqrt(b^2 - x_p.^2))+2*k*(n_co-n_cl)*(sqrt(a^2-x_p.^2))))*exp(1i*((x-x_p).^2)/2*z),-a,a);

那个么我应该怎么做才能得到数值结果呢。因为它包含指数函数和复杂函数，很难计算

---

提前感谢

您可以将区间

（from，to）

中的函数

f

与

积分（f，from，to）

集成。例如，可以使用以下公式计算exp（-x^2）从x=1到5的积分值：

f = @(x) exp(-x^2);
integral(f, 1, 5)