为了计算两个复数的 e+if-(c+id)/(a+b) 本子程序采用如下算法： P-ac.q--bd.-(a-b)(c+d)ww=d+B e=(p-q)fwf=(s-p-q)/w 三，子程序语句 SUBROUTLNE ACJDN(AR.AI.N.BR.Bl.M.L.JS> 四、形参说明 AR,AI一均为双精度实型二维教组，体积为NXN,输入参数。分别存放方程组系 数矩阵的实部与虚部，返回时将被坏。 N一整型变量，输入参数。方程组阶数。 BR,B一均为双精度实型二堆数组，体积为NXM,输入兼行出参数，润月时分别 存放M组右装数向的实部与虚都：返回时分则存放M组解向量的实部与崖部。 M一整型变量，静入参数。方程组右给常数向量的组数。 L一整型变量，输出鑫数。若返回L=0,说明方程组的系数矩阵奇异，求解失败，若 L≠0，表示正常返向。 S一型一维数组，长度为N。本子程序的工作数组。 五、芋程序（文件名，ACDN,FOR) SUBROUTINE ACIDN(AR,AI,N,RR,图，M,L.IS) DIMENSTON AR(N.N).A1(N.N).BR(N.M),BIN.M).JS(N) DOUBLE FRECISION AR.A .Bl.D.P.Q.S D-0.0 DO 101-K.N DO 10J-K.N P-AR(L.J)ARO0J》+M0D￥AIJ) IF OP.GT.D)THEN KK)-I END F 10 CONTINUE W=D IF (W+1.0.EQ.1.0)THEN WRITE(20) 1=0 RETURN FORMAT() 12