view general/numerical/gauss_logpdf.m @ 6:0ce3c2070089

Removed duplicate code and fixed doc in timed_action.

function w=gauss_logpdf(x,sigma)
% gausspdf - Gaussian probability density
%
% gausspdf :: [Size->real] -> [Size->nonneg].
% gausspdf :: [Size->real], nonneg ~'stddev of Guassian' -> [Size->nonneg].
% 	returns normal (gaussian) pdf of 1st argument. 

if nargin<2, sigma=1; end;
w=-0.5*(x./sigma).^2 - (log(sigma)+2*log(pi));