DOUBLE PRECISION B.D.T f-l I(LNE.+L+1)THEN TT=_1 WRITEC*,20) RETURN 15=L+1 D0100K=1,N-1 P=0.0 DO 101=K.LS IF (ABSOBO.1)).GT.P)THEN P=ABS(B(1,1)) T5- END IF 10 CONTINUE· FC+1.QEo.1.0yTEN WRITE(.20) RETURN END IF 20 FORMAT(X,',表FAL,) D030J=1,M T-D(K,D DXK.J)- DIS,D=T 30 T-B(K.J) B(K,J)-B(,J) BOS.J)-T CONTINUE D050J=1,M DK.J)-D(K.J)/B(K.1) 60 B(K,J B(K.D /B(K.1) 090I= T=B1,1) D070J=1,M DO=D(LD-TD(K.D DD90J=2,1 BOJ-1)B(D-T*BCK,D 20
100 CONTINUE (ABS(B(N.1))+1.0.EQ.1.0)THEN T=0 WRT亚(，20) RETURN END IE DO 1dJ=1.M 10 D(N.D- D0150=N-1l,-1 D0120K=1, D0120J=2,Js 120 Dd.K)-Da,K)-BJ)+D(+J-1K IF UIS.NE.IL)I5-=JS+1 150 CONTINUE RETURN END 六.例 设五对角线方程组为 AX-D 其中 -4 1 -2-5 -1 -3 5 4= 3 1 1 -3 13 -31 D- - 31 22 - 21
与带型矩阵A对应的炬形数红招为 0 2 5 6 0 1 -12-3 5 6 -1 -3 -12 -5 6 1 =32 ...B 4 -12 1 -700 在本向题中.《-8，半#究1.一2，希宽【.-5，M一3. 主程序（文件名：ABANDO0.FR)为 ATAB/3.0.2.0,1.0.3.0.-3.0,5.b.-1.0,5.0. -40,a00,505e1.0,8.0,-1.0, -50,1,.-30,1.0.-7.0,30l.020,6.03*20 3￥0.0，-30，】.0，-5.0,90.3￥0.0 DATA DA13.o,-6..31.0,51.0.20d, 220 20.0.17.n.-6.0. .0 41.0.1a0.1.0 21.0l06.0.-05..-21.00时 CALLABAND(B.D.8.2.3.3.) IF IIT.GT.O)THEN WRIJEC.lo)f(Ddl]L8) END n FURMATK1X.6 1) 话钻果为 .674911D15 .100)+0 …,30)-l .3HI I ,7987115 ,2iC014G1 ,制0I:01 ,30000i月··.232川5力14·，13227D11 70u01÷i1 .9G8KI3D-15 .52342D-15 1.7对称方程组的分辉法 一，功能 用分解法求解系效矩养为衬称、月右端具有m钳常数向址的线性代数程组AX一C
二、方法说明 设系数阵A为对，则可以分解为下列形式 ALDI 其中L为单位下三角年库，D为对所钱矩醉。即 d 1 tw…1 d 矩车L和D中的齐元素血下列计算公式癞定： dn=en 4-ad lndadi lu=0.f>i 对干方程组 AX-B 当利D确定后，令 DITX-Y IY-R 解出Y.呼面方秘组 DIX.Y 解出X,其计算公式下： y一修 影=-6-n>时 3.-ylda 三、子程序语句 SUBROUTINE ALDLE(A.N.M.C'.L) 、形参说明 A一双梢度实型维数组，体和为N×N,辅入参数。作嫩对欣阵A,叫时将被 破坏 N一 整型李址，入参数。方程组阶数， M一整理变址，输人参敷。方礼组增常数州址的组数， C一双椭度实孔二推数组，体积为NyM输入兼轴出参数。鸡州时存放M凯 23
向量返回M组解向量。 L 整型变量，箱出意数.若返回L=0,说男子程序工作失敷，若L≠0，表示正常返 五、子超序（文件名：ALDLE.FOR) SUBROUTINE ALDLECA.N.M.C.L) DMMENSION ACN.N).CON.M) L= F (ABS(A(1.1))+1..EQ.L0)THEN L=0 WRITE(+,10) RETURN END IF 10 FORMAT(IX.FAIL) DO201-.N 29A)=A,12/A1,1 D0801±2，N-1 D00J=2,I 30 AGD-A(L.D-A(LJ-1)+A(LJ-D.A(-1J-1) D050K=I+1.N 0401=2,1 40 =AK,D-AKJ-1A,-A0-1-1) F (ABSCA(L.D)+1.0EQ.J.0)THEN L=0 WRTE(",10) RETURN END IE AK.1}=ACK.)fA红，1D 50 70 AN,N)=AN.N门-AONJ-1)sAN,J-1+A0-1,小-1D D001=1.M 001m2.N 00K-2. IJ-CG.D-A(I,K-D+C(K-1J D090I-3,M DO 90 J-1.N 90A(1-1,J2=A-1,1-1)=A0,I-1) IF (ABSCA(N,N))+1.Q 1.0)THEN L=0 24