MATLAB双缝干涉问题求解

根据波的叠加原理写一个双缝干涉实验的模型，我写了一小部分，有一部分不会写

clear all
lambda = 500e-9; %%波长500nm
d = 2e-3; %%双缝间距2mm
D=1; %%双缝距离观察屏之间的距离1m
ym=5lambdaD/d;
xs=ym;
n=101;
ys=linspace(-ym,ym,n); %%观察屏面
for i=1:n
r1=sqrt((ys(i)-d/2).^2+D^2); %%光程1
r2=sqrt((ys(i)+d/2).^2+D^2); %%光程2
phi=2pi(r2-r1)./lambda; %%相位差
B(i,:)=sum(4*cos(phi/2).^2); %%计算光强
end
N=255;
Br=(B/4.0)*N;
subplot(1,2,1) %%一行两列中的第一幅图
image(xs,ys,Br);
colormap(gray(N)); %%使图像以灰度图像显示

有个功能不会写，一个是以结构数组形式定义变量：屏幕和双缝。还有就是屏的位置在空间的角度可调。
想问一下这个功能实现的思路或编码