14.9通数连分式的计算 (5271 1.10厨数曲绘的着出 52 第15章复云 (533 151复散 〔533 1乐？复数底按 (53 15.3发教梁幂-…一 m转 t55) 】5.1复救的N次方根mwm t5a71 5.5复教指数-nn (539) 15.0复整对数 -50i 15.7月致正孩 (543 附录1 FORTRAN?7库管理程序的使用 附录2关于(FORTRAN常用算法程序集K第二意配奢软造的说明一m…一一(548) 参考文献 ·(549)
1.1全选主元高斯消去法 一、能 用全选主元高斯(Gauss)消去法求解线性代数方程组AX=B。 二、方法议明 高斯带去法分两步进行。 第一步消去过程 在这一过程中，为了保证数值计算的稳定性，本子程序采用了全速主元， 对于=1,2，…一1，作以下三步 ()全港主元，即从系数矩阵的第k行，第是列开始的右下子阵中选取绝对值最大的 元素，并将它交换到主元素的位置上。 《2)日一化.即 ,/n→a,=克十1，m,酒 b/au-h (3)消去.即 ,-a6→4，J=k+1,-n 第二步时代过程 a blu-r 〔2) -2ig-121 三、子程序通句 SUBROUTINE AGAUS(A.B.N,X.L.JS) 巴、形参说明 A一—双箱发实型二维数组，体积为NXN,输入鑫数。存放方程组的系数矩阵，返回 时将被破坏。 B—双度实型一数组长度为,输入参数。存放方程组右着向量，这回时将 N一整理变金，输入鸯数。存流方程组的阶数」 &一双精度实型一维数组，长度为N,物出参数.返回方程组的解向. L一盛型变量，输出参数。若返回L=0,说明方程组的系数矩阵奇异，求解失废：若 L≠0，说明正常这回。 1
S一整型一维数组，长度为N。本子程序的工作数组 五，子程序（文件名：AGAUS..FOR) SUBROUTINE AGAUS(A,B.N,X.L.IS) DIMENSION ACN.N].X(N).B(N).JS(N) DOUBLE PRECISION A.B.X.T L1 150K=N- D=0.0 DO 210I-K,N 0210J=KN IF (ABS(A(L,J)).GT.D)THEN D·A5(A1.J) JS(K)-] 210 CONTINUE IF (D+1.0.EQ1.0)THEN 1.=0 ELSE IF (S(K).NE.K)THEN D0220I- A(.K AIK)=AH.s区)》 AUJJSIKD-T CONTINUE ENDIE IF (IS.NE.K)THEN 02301=K T-A(K) A(K)=A(SJ) A(5,J刀=T 230 CONTINUE T-B(K) (IS) BOS)-T END 1F END IF IF (LEQ.0)THEN WRITE(100) RETURN 2
010J-K+1.N A(K.D=A(K.J/A(K.K) 10 E B)-)/A(K.X D0301=K+1,N D020J=k+1,N ACJ)-A(LJ-A(L,K)+ACKJ) 0 CONTINUE B(D=B0-ACKBCK) CON IF (ABS(A(N,N))+1.0.EQ.1.0)THEN 1=0 WRTE(,1Oo】 RETURN END IF X(N)=BOND/ACN.N) D0701=N-1.1.- T=00 CONTINUE X1》=BT)-T 70 CONTINUE 100 FORMAT(IX,'FAIL) ISO)=N 0150K=N.3.-1 F (S(K).NE.K)THEN X(S(K))-T END JF 150 CONTINUE RETURN END 六，创 求解下列方程组 0.2368+02471+0.2568x1+1.2671,=1.8471 0.19681+0.2071z4+1.21684+0.2271z.=1.7471 01582+1.1675a+0.1768+0.1871z4-1.6471 1.1161+0.1254+0.1397+0.1490x-1.5471
主程序（文件名：AGAUS0,FOR)为 DIMENSION A(4).B(.x(4).(4) DOUBLE PRECISION A.B.X DATA/a.8368019680.1582,1.161,0.2471.0,201.1,165.0.1254 ,0.2563,1.8168.0.1758.0.1397,1.21.0.22.0.1R1.0.14c0/ DATA B/1.B471.1.74n1.L.6471,L5471/ CALLAGAUS(A.B.N.X/L.JS) IF (LNE.O)THEN wRTE+,1o4,X0,=1,0 ENO I 10 FORMAT(IX.X(.12)=,DIS.6) END 运行结果为 X415=.30405BD+01 X2】= .9665D+00 930s3D+00 .881297D叶0 1.2全选主元高斯约当消去法 一、功能 用全选主元高斯-纳当(Gau越-Jordar)消去法词时求解系数矩阵相同面常数响量不 同的多军线性代数方组AX=B, 二、方法说明 设妮性代数方程组为 AX=B 其中 X= 9 及，及：…九 全选主元高约当清去法的计算过程如下。 对于-1,2，“，x作知下变换： 1)全主元 4