n=10;

B=[-1*ones(n,1), 3*ones(n,1), -1*ones(n,1)];
A=spdiags(B,[-1,0,1],n,n);
b=[2;ones(n-2,1);2];

correctx=ones(n,1);
x1=jacobi(A,b,5);
x2=jacobi(A,b,10);
x3=jacobi(A,b,50);

keyboard
clear

n=100000;
B=[-1*ones(n,1), 3*ones(n,1), -1*ones(n,1)];
A=spdiags(B,[-1,0,1],n,n);
b=[2;ones(n-2,1);2];
tic
x3=jacobi(A,b,50);
toc