clear?all;
clc;

d?=?4.5;
b?=?1.5;
N?=?5;
lamga?=?0.8;
m?=?1;

for?i=-0.5:0.001:0.5
????seta?=?i*pi;???
????temp_1?=?（sinc（b*sin（seta）/lamga））^2;
????temp_2?=?（N*sinc（N*d*sin（seta）/lamga）/sinc（d*sin（seta）/lamga））^2;
????I（m）?=?temp_1*temp_2;
????
????m?=?m+