function f = g(epsilon)
    global mean sd
    f = lognpdf(epsilon,mean,sd);
end