2、 x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)'-A(i,i+1:n)*x0(i+1:n)').../A(i,i); end ifmax(abs(x-x0))<=tol fid=fopen('G_S_iter_result.txt','wt'); fprintf(fid,'********用Gauss_Seidel迭代求解线性方程组的输出结果********'); fprintf(fid,'迭代次数:%d次',k); f
3、printf(fid,'x的值'); fprintf(fid,'%12.8f',x); break; end k=k+1; x0=x; end ifk==N+1 fid=fopen('G_S_iter_result.txt','wt'); fprintf(fid,'********用Gauss_Seidel迭代求解线性方程组的输出结果********'); fprintf(fid,'迭代次数:%d次',k); fp