High order relaxation schemes for non linear diffusion problems