Many thanks Margot! Other people had also the same suggestion. So nice that you get your answer so fast in this Forum!
--------------------------
Ehsan Motazedi
Wageningen UR
--------------------------
------Original Message------
The inverse gamma probability functions can be found in the package actuar. You can get the pdf by letting x = y^2 where y is a sequence of evenly spaced values over the domain you are interested in and using dinvgamma( x, alpha, beta ) to find the densities associated with the y's. Multiply the density values by 2y for the change of variable from x to y. Here alpha is the shape parameter and beta is the scale parameter.
------------------------------
Margot Tollefson
Consultant
Vanward Statistics
------------------------------