This function calculates the RND from a cross-section of put prices using the PCA method as described in O. Bondarenko. ``Recovering Risk-Neutral Densities: A New Nonparametric Approach, UIC working paper, September 2000. -This function calc