How to find the Dirichlet parameters from Liouville copula with R?

Question

I want to fit data (d=3) to a Liouville copula using this function from lcopula package:

liouv.maxim(data, family, interval, boundary = NU