３重対角行列に対するGaussの消去法のプログラム

Question

gaussの消去法に関するプログラムは完成したのですが、
三重対角行列に対するプログラムの作り方が分かりません。
以下のプログラムで、ここを直したら三角対角行列の
プログラムになるということを教えていただけると幸いです。

#include<stdio.h>
#define N
void gauss_elm(double a[N][N],double b[N],int n);
int main(void){
double a[N][N]={{},{},};
double b[N]={};
gauss_elm(a,b,N);
return 0;
}

void gauss_elm(double a[N][N],double b[N],int n){
int i;
int j;
int k;
double pivot,m;
double x[N];
for(k=0;k<n-1;k++)
{
pivot=a[k][k];
for(i=k+1;i<=n-1;i++)
{
m[i][k]=a[i][k]/pivot;
for(j=k+1;j<=n-1;j++){
a[i][j]=a[i][j]-ma[k][j];
}
b[i]=b[i]-mb[k];
}
}

for(i=n-1;i>=0;i--){
  if(i!=n-1){
        for (k=n-1;k>i;k--){
             b[i]=b[i] - a[i][k]*b[k];
        }
    }
    b[i]=b[i]/a[i][i];
}

for(i=0;i<N;i++)
{
    printf("%f\n",b[i]);
}

}

Answer 1

三重対角行列とか、三角対角行列の定義がよく分かりませんが、下記でどうでしょうか？
http://www.oyama-ct.ac.jp/M/yamashita/cgaus.htm