function ll = invg_loglikelihood(x, p)
    mu = p(1);
    lambda = p(2);
    ll = .5*(log(lambda) - log(2*pi) - 3*log(x)) - ...
        lambda*(x-mu).^2./(2*mu^2*x);