mdf_ig function

Marginal Density for Given Scale Parameter and Inverse Gamma Prior for $\tau^2$