N=10;
wavelength=1;
d=0.5;
theta=-pi/2:0.01:pi/2;
theta0=30;
b=2*pi*d*sin（theta）./wavelength;
b1=2*pi*d*sin（pi*theta0/180）./wavelength;
y=abs（sin（N/2*（b-b1）