clc
A=[5 -2 3 0;-3 9 1 -2;2 -1 -7 1; 4 3 -5 7];
b=[-1 2 3 0.5]';
x0=[0 0 0 0]';
tol=1e-5; itr=0;
[x]=jacobi_dejen(A,b,x0,tol,itr);