function y = rcos(x,alpha,T);

% function y = rcos(x,alpha,T)
% raised cosine pulse genaration
% T - period or symbol interval 
% alpha - rolloff 0 < alpha < 1
% x - time in analogue form

y = ones(size(x));
i = find(x);
y(i) = sinc(x(i)/T).*cos(alpha*pi*x(i)/T)./(1 - (2*alpha*x(i)).^2);