Recovering risk-neutral probability density functions from options prices using cubic splines and ensuring nonnegativity