I've been messing around with my code for a while, but the x value it outputs never comes out correct.
function [x] = jacobi(A,b,x0,n_iter)%A=[10 -1 2 0;-1 11 -1 3;0 3 -1 8; 2 -1 10 -1];
%x0=[1;1;1;1];
%b=[6 25 -11 15];
%n_iter=10;
n=length(x0);x=zeros(size(x0));for k=1:1:n_iter for i=1:1:n temp=0; for j=1:1:n temp=temp+A(i,j)*x(j); end temp=temp-A(i,i)*x0(i); x(i)=(1/A(i,i))*(b(i)-temp); end x0=x; end disp(x)
I'm just unsure which line I have messed up.
Best Answer