fourierD.m

clf reset
clear all;
clc
T = 5;
t = 0:.1:T; 
x = zeros(size(t));

for k=1:2:19, 
   x = x + sin(k*t)/k; 
   plot(x)%, 'linewidth', 4)
   pause

end;