Vectorize:
I=((sin(N*alpha))./(N*sin(alpha))).^2;
↑ ↑
HERE AND HERE
Use element-wise operations, .* instead of *, and .^ instead of ^.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!