A relaxation approximation for degenerate parabolic equations