function [w0]= cp0201_waveform(fc,Tm,tau);
dt = 1 / fc;           
OVER = floor(Tm/dt);    
e = mod(OVER,2);
kbk = floor(OVER/2);
tmp = linspace(dt,Tm/2,kbk);
s = (1-4.*pi.*((tmp./tau).^2)).* ...
     exp(-2.*pi.*((tmp./tau).^2));

if e           
    for k=1:length(s)
        y(kbk+1)=1;
        y(kbk+1+k)=s(k);
        y(kbk+1-k)=s(k);
    end
else              
    for k=1:length(s)
        y(kbk+k)=s(k);
        y(kbk+1-k)=s(k);
    end
end

E = sum((y.^2).*dt);    
w0 = y ./ (E^0.5);