Numerical methods for fourth order nonlinear degenerate diffusion problems