High-order relaxation schemes for nonlinear degenerate diffusion problems