a=-pi/4:pi/999:pi/4;
n=2:2:1000;
for ii = 1:length(a)
y1(ii)=0.3/sum((cos(n.*(pi/9))./((n.^2-1).^2))).*sum(cos(n.*(pi/9)).*cos(n.*a(ii))./((n.^2-1).^2));
end
plot(a,y1);
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!