compute the pole density function for a given unimodal component
pdf = calcPDF(component,h,r,varargin)