% File: Fig3_26b.M

clear;

r = input('Enter value of r between 0 and 1:');
fo = 1;
B = fo*(1+r);
fd = B*r/(1+r);
td = 3/(100*fo);
t = -100*td:td:100*td;

xo = 2*pi*fo*t;
xd = 2*pi*fd*t;

temp1 = Sa(xo);
he = zeros(length(t),1);

for (i = 1:1:length(t))
  he(i) = 2*fo*temp1(i)*cos(xd(i))/(1-(4*fd*t(i))^2);
end;

plot(t,he);
xlabel('t');
ylabel('he(t)');
title(['r=',num2str(r),'    NOTE: See Fig. 3-26b']);
grid;




