% Function to compute cosine bell
%
% Syntax: prof = cbell(ns)
%
% Where ns is the number of samples.
%
% By D. Holmgren holmgren@phobos.astro.uwo.ca

function prof = cbell(ns)
prof = ones(ns,1);
for k = 1:10
 indx = ns + 1 - k;
 w = 0.5 * ( 1 - cos(pi*k/10));
 prof(k) = prof(k)*w;
 prof(indx) = prof(indx)*w;
end