function [H] = square_wave_harmonic(i, t, w0)
% Function to evaluate a Fourier Series term (harmonic) for the square wave
B = square_wave_amp(i);
H = B * sin(i*w0*t);
end