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