Hi Kristinn, do you mean something like this:
t = 0:.02:2*pi; % Graph time, interval between 0 and 2pi
f = 0*t; % Make function variable same as time
hold on
plot_n = 5;
for n = 1:1:1000
if mod(n,2)==1 % only use when n is a odd number
b_n = ((4)/(n*pi)); % calc of bn
f_n = b_n*sin(n*t); % calc of bn times summotation
f = f + f_n; % add to last
end
if ~mod(n,plot_n)
plot(t,f)
end
end
grid on;
xlabel('t');
ylabel('f(t)');
title(strcat('The building of a square wave: Gibbs'' effect'));
