Check this code
ao=1/2;
an=0;
%% Problem 1_Practice_Problem_17.1
disp('Problem 1_Practice_Problem_17.1')
N = 10000;
Wo = pi;
t = -1:.001:11;
f = zeros(1,length(t));
for n = 1:2:N
bn = (1-(-1)^n)./(n*pi);
f = f + an*cos(Wo*n*t) + bn*sin(Wo*n*t);
end
% Plotting and plot formatting is done here
plot(t,f)
axis([-1 11 -1.1 1.1])
title(['Fourier Series Representation w/ ' num2str(N) ' Terms'],'FontSize',16)
xlabel('Time')
ylabel('Function Value')