杨氏干涉----MATLAB

lamda=500e-9;
d=2e-3;
D=1;
ym=5lamdaD/d;xm=ym;
n=101;
ys=linspace(-ym,ym,n);
for i=1:n
r1=sqrt((ys(i)-d/2).^2+D2);
r2=sqrt((ys(i)+d/2).^2+D2);
phi=2pi(r2-r1)./lamda;
B(i,:)=sum(4*cos(phi/2).^2);
end
N=255;
Br=(B/4.0)*N;
subplot(1,2,1)
image(xm,ys,Br)
colormap(gray(N));
subplot(1,2,2)
plot(B,ys)