Since the object returned by the paretotails is of Piecewise Distribution type, you can directly use the pdf function to query the pdf values at any point. Have a look at an example usage below:
figure, plot(Xgrid_data, pdf(pd_fit,Xgrid_data))
Here’s the documentation of pdf for your reference: