In this paper, the numerical solution for two-dimensional nonlinear parabolic equations is studied using an alternating-direction implicit (ADI) Crank–Nicolson (CN) difference scheme. Firstly, we use the CN format in the time direction, and then use the CN format in the space direction before discretizing the second-order center difference quotient. In addition, we strictly prove that the proposed ADI difference scheme has unique solvability and is unconditionally stable and convergent. The extrapolation method is further applied to improve the numerical solution accuracy. Finally, two numerical examples are given to verify our theoretical results.