function [G_vec] = servant(num_term, t_vec, period)
num_sample = length(t_vec);
coeff = [1:2:(2*num_term+1)]';
sine_matrix = sin(coeff * t_vec * pi / period);
recip = 1/coeff;
G_vec = 4.0 * recip * sine_matrix / pi;
