Two algorithms, and corresponding FORTRAN computer programs, for the computation of posterior moments and densities using the principle of importance sampling are described in detail. The first algorithm makes use of a multivariate Student t importance function as approximation of the posterior. It can be applied when the integrand is moderately skew. The second algorithm makes use of a decomposition: a multivariate normal importance function is used to generate directions (lines) and one dimensional classical quadrature is used to evaluate the integrals defined on the generated lines. The second algorithm can be used in cases where the integrand is possibly very skew in any direction.