du-arti. (2)进行国代。眼 dfa→x d-a441→4话丝用-1，，2，l 在本子程序中，为了便于三对角矩陈的输入，用…个实型一维数组B们：3跟一2)以行 为亦对角矩阵A中三条对角线上的元素。其中为矩阵A的阶数。三对角线矩阵 A与一…维敢组B之间有知下关系： B[2-1)+],i-16i+1 Ai.p-o ,其它 根据以上关系，可以得到如下计算过程 (1)对于=1,2，，n一1作的一化及消元处理.即 B3-t)/B(3-2)→B3-1) d/B(3k-2D→d B(3+1)-B(3)B(3k-1)→(3最+1) d-B(3k)d,→d (2)进行回代即 d/83x-2→d, d4-B(3t-1d41→d--1m,2】 录后的解向量在D中 特别要指出，由于途赶法本质上是没有选士元的高斯清去法，因此其计算过程是不移 定的。包当三对免矩年情是下列条件： leal≥lai+leld=23m-1 au>a 时，追赶法的计算过程中不会出现中间结果数量级的巨大增长和舍入误差的严置积。 三、子程序诺句 SUBROUTINE ATRDECB.N.M.D.L) 四、形参说明 议精度实型一维数组，长度为M=3:N-2,输人参教。以行为主存放三对角 知阵A中三条对角线上的元素，即在B中依次传敬元素， 时滨数组将被酸坏】 N…藤型变螫输入参教，方程组阶数 M 鉴型变地，输入参数。M一3N一2为三对角矩昨三条对角线上的元素个 双格度实型一维数组，长度为N,输入兼输出多数。闾用时存放方程组右端的 常数向世：返回方程组的解响址。 整型变量，输出参数。若惑1.<0，说明M的慎不正确（应为M=3N-2 若L-0,说明程序工作失效：者L>0,表示正常道回。 16