程序名稱
　　　subroutine bandv(n,m,np,a,b,id,ir)――大型對稱變帶寬方程組的解法 　
　 功　能
　　　本程序適用于求解大型線性代數方程組
　　 的解，其中A為n×n階大型對稱變帶寬帶型矩陣，X和B均為n×m階矩陣，為了節省存儲量，對A采取變帶寬存儲，即每行只從第一個非零元素到對角線為止。
　 使用說明
　　　 子程序語句 subroutine bandv(n,m,np,a,b,id,ir)
　　　 參數說明
　　　　◎ n 整變量，輸入參數，方程組的階數。
　　　　◎ m 整變量，輸入參數，方程組右端常向量的個數。
　　　　◎ np 整變量，輸入參數，為系數矩陣A壓縮存儲的元素存儲的元素的個數。
　　　　◎ a 輸入參數，np個元素的一維數組，存放系數矩陣的壓縮存儲元素。
　　　　◎ b 輸入、輸出參數，n×m個元素的二維實數組。開始時存放方程組右端的m個n維常向量；工作結束時，存
　　　　　放m個解向量，元素按列存放。
　　　　◎ id 輸入參數，n個元素的一維數組。存放系數矩陣的各個對角線元素在壓縮存儲的一維數組中的位置。
　　　　◎ ir 整變量，輸出參數，標志。若輸出值不為零，表示矩陣奇異，不求方程組的解；正常求解結束時值為
　　　　　 零。
　 方法簡介
　　　 其基本方法就是三角分解法：
　　　　　　　　　
　　　　只是其中的元素Kij在一維等帶寬存儲A中的位置為A（n）
　　　 分解得到L和D后，回代的求解步驟為：
　　　　◎ 由LV＝P向前回代求解V；
　　　　◎ 求對角陣D的逆陣D－1；
　　　　◎ 由V'＝ D－1V求修正分解后的載荷陣V'；
　　　　◎ 由LTX＝V'向后回代求解基本未知量X。