APM1513- CODES WORTH MEMORIZING
THE JACOBI CODE: please see study guide for complete discussion
Function xnew=Jacobi(A,b,xold)
n=size(A1); extract size of the matrix A from input
At=A;
Xnew=xold; x1 is entered into its place when determining x2.
For k=1:n Work through all the equations
At(k,k)=0 when calc. x1, we made x2 and x3=0
Xnew(k)=(b(k)-At(k,:)*xold)/A(k,k)
b(k)-> vector of right-hand side
At(k,:)-> extract row
A(k,k)-> coefficient of the X we are looking at
End
Endfunction
GAUSS-SEIDEL CODE: please see study guide for complete discussion.
Function xnew=Jacobi(A,b,xold)
n=size(A1);
At=A;
Xnew=xold;
For k=1:n
At(k,k)=0
Xnew(k)=(b(k)-At(k,:)*xnew)/A(k,k) xnew is the only difference.
End
Endfunction
ITERATIVE LINEAR SOLVING METHOD: Please see study guide for complete discussion