Metode Gauss-Seidel

function [X1,g,H] = seidel(A,b,X0,T,N)