function y = RCpulse(alpha, t, Ts) 
% Inputs: alpha



t = -6:0.01:6;
t = t + eps; % Offset slightly to prevent NaNs.

a = 0;

tau = 1; % Symbol time
tpi = pi/tau;
atpi = tpi*a;
at = 4*a^2/tau^2;

y = sin(tpi*t).*cos(atpi*t)./(tpi*t.*(1-at*t.^2));

plot(t, y)
grid on