How to solve two-dimensional parabolic equations with variable coefficients by ADI(Alternating direction implicit) method?

Question

I tried to solve the two-dimensional parabolic equation with variable coefficients using the ADI method. The equation is as follows: [1]: https://i.stack.imgur.com/bYkfo.png