gauss

    連立１次方程式：
    [aij] [xj] = [ci]
    を Gauss の消去法を用いて解く倍精度サブルーチン．

  • サブルーチン名： gauss

  • 呼び出し方：call gauss(adim,cdim,ndim,ndata,epsln,xsol,wrk1,wrk2)

  • 引　数：
    引数

    意　味

    入出力の区分

    備　考
    adim(ndim,*)係数行列 [aij]を格納しておく配列の配列名入力
    cdim(*)非同次ベクトル [ci] を格納しておく配列の配列名入力
    ndim配列 adim の整合寸法入力
    ndata未知数の個数入力ndata≦ndim でなければならない．
    epsln特異性判定子入力1.0×10-14 程度．マシン・イプシロンでも良い．
    xsol(*)連立１次方程式の解 [xi] を格納する配列の配列名出力寸法は ndata 以上
    wrk1(ndim,*)作業用２次元配列．整合寸法は ndim．作業用配列 adim と同じ寸法に宣言しなければならない．
    wrk2(*)作業用１次元配列．作業用配列 cdim と同じ寸法に宣言しなければならない．