Matlab code for solving 2D Laplace pde

I already discretized pde by central difference method, but not sure how to create a code for iterative methods such as Jacobi and Gauss Seidel.

Answer 1

[Precise Simulation]

Sparse linear solvers are already available in Matlab, the backslash operator per default uses the Umfpack direct solver

   y = A\x;
   help \

for iterative solvers you can for example use gmres or bicg

   help gmres
   help bicg

Answer 2

[Vishwanatha S]

Gauss Seidel iteration