Try this —
t = linspace(0, 2*pi, 1E3); % Time Vector
x = 1:2:31; % Terms Vector
f = 2; % Frequency
s = sum(sin(x(:)*t*f)./x(:)); % Signal
figure
plot(t,s)
grid
xlabel('t')
ylabel('s')
xlim([min(t) max(t)])
I will let you figure out how it works.
.