computeprob function

Computes the probability that a model parameter is <=k from inla marginals