Ето пълно аналитично решение. Малко съм фризирал функциите за хубост.

f(1, Y) = 2 при 1/2<=Y<=1

        = 0 навсякъде другаде



За n>1:

f(n, Y) = Интеграл в граници (1/n, y/(1-y)) от f(n-1, x)*(1+1/x)*dx при 1/(n+1)<=Y<=1/n

        = Интеграл в граници (1/n, y/(1-y)) от f(n-1, x)*dx +

         +Интеграл в граници (y, y/(1-y)) от f(n-1, x)/x*dx при 1/n<=Y<=1/2

        = 1 + Интеграл в граници (y, 1) от f(n-1, x)/x*dx при 1/2<=Y<=1

        = 0 навсякъде другаде