Matlab/interacni metoda – Matematika – Fórum – Programujte.com

Ahoj, mam program v matlabu, ktery pocita interacni jaccobiho metodou.

pro vstup A=[10 -2 -2; -1 10 -2; -1 -2 10]

               b=[6 7 8]

               x=[0 0 0]

               k= 3

vypocita, ze x1= 0,6; 0,7; 0,8 ....ale nevim jak udelat, abych spocital i x2 a x3 .....vim, ze x2 se ma rovnat 0,90; 0,92; 0,93 a x3=0,970; 0,976; 0,982 .... algoritmus spocita pouze x1, prosim poradte, kde mam chybu

function x=JACOBI(A,b,x,m)
n=size(A,1)
D=diag(diag(A));
U=triu(A,1);
L=tril(A,-1);
H=-inv(D)*(L+U);
B=inv(D)*b;

for l=1:m
    for i=1:n
	for j=1:n
      x=B+H(i,j)*x(j);  
     
    end
    end
end