首先，这是一个家庭作业，我实际上正在尝试复制Young的狭缝实验，我意识到这可能属于物理论坛，但我的问题是涉及数值近似，我不能做工作
我想要做的是拥有四个变量的函数，（x，y，xp，yp），我想整合那些关于x和y的函数，从-infinity到infinity，这意味着：
F（X，Y，XP，YP） - ＆GT; Q（XP，YP）
我的f是一个看起来像我的代码中的函数f的函数。当我有q = @（xp，yp）时，我想放一个常数yp，让我们说yp = 0 - ＆gt; b = @（xp）q（xp，0），然后我想针对xp绘制b，您可以查看下面的代码，看看我想做什么，但显然它不起作用。在我的代码中，我实际上是在-50到50之间进行集成，只是为了让程序运行得更快。
clear all;
clc
l=1.45;
lambda = 650*10^-9;
f = @(x,y,xp,yp) exp((1i*2*pi/lambda)*sqrt((xp-x).^2 +(yp-y).^2 +l.^2));
xp = linspace(-10,10,20)
q = @(xp,yp) integral2(f,-50,50,-50,50)
b=@(xp) q(xp,0)
plot(xp,abs(b).^2)
我似乎无法完成这项工作，有人对我有任何建议吗？