main.m

gauss-seidel算法实现对方程组的求解

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)