probleme de jacobi programme en matlab

voila le progrmme que j'ai fais,je vx juste votre aide la solution doit etre s=[1 1 1];
je vx l'erreur et lacorrection de cette erreur etje vx aussi l'algorithme ou le programme de crank nicolson resolution des equations aux dérivé partielle.

clc;

clear;

disp(

' %%%%%%%%%%%%% programme resolution AX=B %%%%%%%%%%%%%');

disp(

' %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%');

a=[8 1 -1;2 1 9;1 -7 2];

b=[8 12 -4];

a=[a b'];

n=3;

eps=1E-6;

x=[0 0 0];

s=zeros(1,n);

tol=1;

iter=1;

while

tol>eps

for

i=1:n

s(i)=0;

for j=1:i-1

s(i)=s(i)+a(i,j)*x(j);

end

for j=i+1:n

s(i)=s(i)+a(i,j)*x(j);

end

s(i)=(a(i,n+1)-s(i))/a(i,i);

end

tol=norm(s-x);

iter=iter+1;

end

iter

s'

T'es pas dans le bon forum je crois...

Cordialement,
Sébastien.