function p = pdf_nor(x,mu,la,id)

if nargin<4
    id=1;
end
if id~=1 && id~=2
    error('specify either var (1) or prec (2)')
end

if id==1
    pre=1/la;
else
    pre=la;
end

p = 0.5*log(pre) - 0.5*log(2*pi) - pre/2*(x-mu).^2;