%This function takes draws from the inverted gamma density.
%The inputs are the two parameters alpha, and beta, and 
%the function returns a n times m matrix. 
%The inverted gamma density function is
%f(x) propto x^-(alpha + 1) exp( -1/(beta*x) );

function [igdraws] = invgamrnd(alpha,beta,n,m);

tempp = gamrnd(alpha,beta,n,m);
igdraws = 1./tempp;
