Please, write a program in Python and explain all steps!

we will solve a problem from electrostatics. In particular, we consider a capacitor and want to solve the Laplace equation in two dimensions in order to obtain the electric potential. To this end, we will use the finite differences method together with an iterative procedure for solving the system of equations. In particular, the setting is as follows: We are looking at the cross-section of a plate capacitor and therefore consider a domain [0,20][0,10] that contains no charges, where the electrodes of the capacitor (being modeled as two infinitely thin plates) are located at [4,16]4 and [4,16]6. The potentials at the electrodes are 3V and 3V, respectively. The potential is thus obtained from solving =0 in the domain (apart from the locations of the capacitor plates). In order to be able to solve the problem, we also need to provide conditions at the boundary. For simplicity, we employ homogeneous Dirichlet conditions, i.e. D=0. For the whole setup, see the sketch in Fig. 1. Write a function SORmethod (A, b, x0, omega, tol) that calculates the solution of the problem Ax=b, using the SOR method based on the Gauss-Seidel method. Here, A is the matrix, b is the right-hand-side vector, x0 is the initial guess, omega is the relaxation parameter, tol is the allowed tolerance. The function should return the approximate solution vector of the problem. Test your function with a small matrix problem, for which you can calculate the solution manually. Remark: Since we are going to look at rather big matrices, it might be necessary to use a sparse matrix format for A