使用matlab实现Jacobi迭代法与Gauss-Seidel迭代法求解线性方程组

jacobi.m

function [y,n]=jacobi(A,b,x0,ep) %A:系数矩阵；b:常数矩阵；x0:迭代初值；ep:迭代精度；
D=diag(diag(A));	%取A的主对角线上的元素建立对角矩阵；
L=-tril(A,-1);		%取A的主对角线以下（不包括主对角线）的元素建立下三角矩阵；
U=-triu(A,1);		%取A的主对角线以上（不包括主对角线）的元素建立上三角矩阵；
B=D\(L+U);
f=D\b;
y=B*x0+f;
n=1;