function out = c1(sd1)

global sn;
global alphaP;

s1 = sn .* (1 ./ (2.*sqrt(2.*pi).*sd1.^2));
out = s1 .* sqrt(-2.*log(alphaP));

%out = 5 * out;

