Using ADI method to solve 2D Laplace equation

Could anyone help me to code this problem, i solve this qeustion by , Jacobi, gauss seidel and SOR. but i am stucking in ADI coding

Consider Laplace’s equation in the rectangular domain with Dirichlet

boundary conditions specified in Figure 1. Use square mesh (h = k), and start with 100 subdivisions of the vertical interval (0,1). Calculate the solution using ADI methods, comparing the convergence rate (in terms of the number of iterations).