程序名稱
　　　 gs1(n,n1,a,ep,kwji)――Gauss消去法
　 功　能
　　　本程序用高斯消去法求線性方程組
　　　的解，其中A為n×n的系數矩陣，x為解向量，b為方程組右端n為列向量。
　 使用說明
　　　 子程序語句 subroutine gs1(n,n1,a,ep,kwji)
　　　 參數說明
　　　　◎ a 輸入參數，n×n1個元素的二維實數組，存放由系數矩陣和右端列向量所組成的增廣矩陣，存儲方式是按
　　　　　 列存放，亦是輸出參數，其n1列存放解向量。
　　　　◎ n 整變量，輸入參數，方程組的階數。
　　　　◎ ep 實變量，輸入參數，作控制常數用，通常為比較小的正實數。
　　　　◎ kwij 整變量，輸出參數，標志。當矩陣A（k－1）（其意義見方法簡介）的第k列的元素的絕對值均小于
　　　　　　ep時，令kwji＝1，這時認為方程組無解；否則令kwji＝0。
　 方法簡介
　　　高斯消去法由兩部分組成，即向前過程和回代過程。
　　　向前過程由n－1步組成。假設向前過程以進行k－1步，得到矩陣A（k－1），它有
　　　　　　　　 　　　
　　　 式： 右端項為
　　　　　　　　　　　　　　　
　　　 第k步得到的矩陣A（k）及右端項b（k）的元素由下列公式確定
　　　　　　　　　　　　　
　　回代過程的目的是求出方程組的解。其計算公式如下：