%--------------------------------------------------------------------------
% function f = gauss(para,x)    returns gaussian distribution
%
% x     scalar variable
% para  [s,c,z] 
%   s sigma, 
%   c amplitude, 
%   z shift along x
%
% the gaussian distribution is defined by
%     f(x) = 1/sqrt(2 pi)/sigma exp(-(x-z)^2/ (2 sigma^2) ) 
%
%  2-May-2003 S.Schreiber
%--------------------------------------------------------------------------
function f = gauss(para,x)
% ini f
f = zeros(size(x));

if isempty(para) | para(1) <= 0
    return
end

s = para(1);    % sigma
c = para(2);    % amplitude scaling
z = para(3);    % shift along x

f = c/sqrt(2.*pi)/s*exp(-(x-z).^2/2/s^2);

return