1-4-3-GS迭代法-程序框图.pdf

课程教案 代数方程及代数方程组的求解在化学中的应用 Gauss—Seidel迭代法 ——程序框图 Gauss—Seidel迭代法——程序框图 GSD(A,B,N,X,EPS,MAXIT) n no  aij  aii j 1,j i yes 初值X(I)=0 I=1,2,…,N IT=0 1 J1 N X (B  A *X(K1)  A *X(k) ) (I 1,2,, N) I I  IJ J  IJ J AII J 1 J I1 (K) (K1) no max X(K1) X(K) EPS X X I I I I 1IN A：数组A(N,N)，存放线性方 程组矩阵的系数； yes B：数组B（N），存放线性方 IT=IT+1 程组右端的系数 N：线性方程组的维数； X：数组X（N），存放线性方 no 程组的根； IT>MAXIT EPS: 计算精度。 MAXIT：最大迭代次数 yes RETURN 有缘学习更多+谓ygd3076考证资料或关注桃报：奉献教育（店铺） Gauss—Seidel迭代法——子程序 SUBROUTINE GSD(A,B,N,X,EPS,MAXIT) G-S子程序开始 DOUBLE PRECISION A,B,X,EPS DIMENSION A(N,N),B(N),X(N) 定义实型变量，双精度 DO 20 I=1,N S=0.0 定义数组，存放方程组系数 DO 10 J=1,N IF(J.EQ.I)GOTO 10 判断是否主对角占优势 S=S+ABS(A(I,J)) 10 CONTINUE IF (S.GE.ABS(A