Intensitätsverteilung nach Gitter mit N Spalten.
Sin(N x)^2/Sin(x)^2

Manipulate[Plot[Sin[N x]^2/Sin[x]^2,{x,-1.5 Pi,1.5 Pi}, PlotRange -> Full],{N,2,20,1}]