cosmo pca hdr

function [pca_samples, params] = cosmo_pca(samples, retain)