Transformation formulae are given in most theoretical statistics books. A simple case is a continuous one-to-one transformation (like ln) of a continuous distribution. We just divide by the absolute derivative, |dy/dx|.
In your case, dy/dx = 1/x = 1/exp(y)
fX(x) = xN-1exp(-x/K)/(KNGamma(N))
fY(y) = exp(Ny - exp(y)/K)/(KNGamma(N))
I think -Y has a Gumbel distribution which is one of the limiting extreme value distributions.
I do not have an answer for you, Huda. I do know that the Gumbel distribution is used to approximate p-values for the likelihood ratio test in the satscan cluster analysis software.
Why are you considering ln(X)?
I have found this, though: http://www.math.wm.edu/~leemis/chart/UDR/PDFs/GammaLoggamma.pdf
Many thanks for all. According to one-one transformation, Y is distributed as Log Gamma, the chart in this URL clarifies Univariate Distribution Relationships
which is the standardised form of the Gumbel distribution. So Alexander and Eliardo are correct. The log of a gamma is a slight generalisation of the Gumbel. As it is the log of a gamma, the term log gamma seems appropriate.
(Terminology can be confusing--a lognormal isn't the log of a normal, instead its log is normal--in fact you can't take the log of a normal because the log of a negative value is a complex number).
As there are two variables, this is a little more advanced. Depending on how you solve it, it can be a related problem. You can transform X and Y to U = -ln X, V = -ln Y so that XY < z becomes U + V > -ln z. Having obtained the distributions of U and V we have a two variable problem.
We can also solve it directly. We want 1 - the area between the lines y = 1, x = 1 and the hyperbola, xy = z. Alternatively it is the complemetary area which can be broken up as a rectangle with area z and the area between z and 1 under the hyperbola.
Huda, I'll leave you to do this both ways. It will help contribute to your expertise in transfromations.
Sorry for a misprint in the left hand side of the equality. The derivative operation was missed and correct statement is P'( X Y < z ) = - ln z. Demo and the Monte-Carlo test are in the attached file.