High-Order Relaxation Schemes for Nonlinear Degenerate Diffusion Problems