程序名稱
　　　 ldlt(a,n,iwb,p) ――對(duì)稱帶型方程組的解法
　 功　能
　　　本程序用以求解帶型方程組
　　　其中K（N，IWB）為等帶寬存儲(chǔ)的帶型矩陣，IWB為半帶寬。
　 使用說明
　　　 子程序語句 subroutine ldlt(a,n,iwb,p)
　　　 參數(shù)說明
　　　　◎ a 輸入?yún)?shù)，n×iwb個(gè)元素的二維實(shí)數(shù)組，存放系數(shù)矩陣K的上半帶區(qū)內(nèi)的元 素。
　　　　◎ n 整變量，輸入?yún)?shù)，方程組的階數(shù)。
　　　　◎ iwb 整變量，輸入?yún)?shù)，半帶寬。
　　　　◎ p 輸入、輸出參數(shù)，載荷陣。開始時(shí)存放向量P；結(jié)束時(shí)存放解向量X。
　 方法簡介
　　　二維等帶寬存儲(chǔ)的三角分解遞推公式為（半帶寬為D）：
　　　　　　　　 　　　
　　　 行列的循環(huán)周界為：
　　　　　　　　　　　　　　
　　　回代的求解步驟為：
　　　　◎ 由LV＝P向前回代求解V；
　　　　◎ 求對(duì)角陣D的逆陣D－1；
　　　　◎ 由V'＝ D－1V求修正分解后的載荷陣V'；
　　　　◎ 由LTX＝V'向后回代求解基本未知量X。