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