function prob = invgausspdf(x, p)

    mu = p(1);
    lambda = p(2);
    
    prob = sqrt(lambda/2/pi./x.^3) .* exp(-lambda *(x-mu).^2/2/mu^2./x);