fitting pareto to data?

hi

i have an experimental dataset which i think is coming
from a pareto/powerlaw distribution. 

i am trying to use the 2 parameter pareto distribution
from Jim Lindsey's rmutil package (i am open to using
any package however.)

i would like R to compute estimates for the 2 pareto
parameters based on my dataset; what commands would
achieve this?  

thanks
-gong