Code to call C function, paretodens.c, in R