% Example 4.14 
%
w = -50:.3:50;
tau = 1;
V = tau/2*sin(tau*w/4).^2./(tau*w/4).^2;
% alternately, use V1 = tau/2*sinc(tau*w/4/pi).^2;
clg
subplot(211);plot(w,V)
xlabel('Frequency (rad/sec)')
ylabel('V') 
title('Figure 4.30')
subplot(111)