%%% Example 3.5 in O/W/N

Kmax=50;
k=-Kmax:Kmax;

ak=sin(pi*k/3)./(k*pi);
ak(Kmax+1)=1/3;

figure(1)
orient tall
subplot(311)
stem(k,abs(ak),'filled')
axis([-Kmax Kmax 0 .5])
xlabel('k')
ylabel('|a_k|')

subplot(312)
stem(k,angle(ak),'filled')
axis([-Kmax Kmax -pi pi])
xlabel('k')
ylabel('\angle a_k')

figure(2)
orient tall
subplot(311)
[x,t]=ctfs_synthesis(3,ak);
plot(t,x)
axis([-3 3 -.5 1.5])
grid on
xlabel('t')
ylabel('x(t)')