function [x,err,it]=GJ(A,b,x0,maxit,tol)
%
% En enkel implementasjon av Gauss-Jacobi
% [x,err,it]=GJ(A,b,x0,maxit,tol)
%
% Output: Løsning, feilestimat og antall iterasjoner brukt
% Input:  Matrise, høyreside, startverdi, maksantall iterasjoner og
% toleranse

D=diag(diag(A));
N=A-D;
it=0;err=2*tol;x=x0;

while it<maxit & err>tol
    it=it+1;
    x_ny=D\(b-N*x);
    err=norm(x_ny-x);
    x=x_ny;
end



