function [F] = square_wave_series(t, n, w0)
% Function to evaluate a square wave Fourier Series to a given number of
% terms
F = 0;
for i = 1 : n
    F = F + 4/((2*i-1)*pi)*sin((2*i-1)*w0*t);
end
end