% File: E3_045.M

clear;

Tb = 0.25;
R = 1/Tb;
f = 0:0.2:10;
x =pi*Tb/2*f;
temp1 = Sa(x);
temp2 = sin(2*x);
temp2 = temp2(:);

P = Tb*(temp1.*temp1).*(temp2.*temp2);

plot(f,P);
xlabel('f');
ylabel('Magnitude of P(f)');
title('NOTE: Compare with Figure 3-16d');
grid;




