HP50g Using library functions in programs

[erazor]

So I have the formula n=log(g/(1-w))/log(g/(g-1)) which estimates the number n of needed clones to pick from an amplified genetic library to have all g different clones with 1 > w >= 0,9 probability.
I need to use this with g values above 1e14. That doesn't work with HPs 12 digits and produces wrong results with, for example, a 15 digits casio.
I transferred this to the following user rpl program:

<< -> G W <<
G R<-->F 'G' STO
W R<-->F 'W' STO
G 1 W FSUB FDIV FLN
G G 1 FSUB FDIV FLN
FDIV
R<-->F
>>
>>

which now works well. For calculating w I would transpose the formula, but just for learing purpose as you mentioned it:
How do I need to formulate the program so I can use it with the numeric solver like the formula? The Num.solv lets load me the program but doesn't work otherwise and I found nothing in the manual about this.