cosmo corr hdr

function c = cosmo_corr(x, y, corr_type)