main.m

来自「gauss-seidel算法实现对方程组的求解」· M 代码 · 共 6 行

A=[-4,1,1,1;1,-4,1,1;1,1,-4,1;1,1,1,-4];
b=[1,1,1,1,]';
x0=[0,0,0,0]';
emg=10^-5;
N=100;
[x,k]=gaussseidel(A,b,x0,N,emg)