三对角矩阵

#include<cstdio>
#include<iostream>
#include<cstring>
const int N=1000+100;
void tdma(int n,double aa[N],double bb[N],double cc[N],double xx[N]){
    cc[0]=cc[0]/bb[0];
    xx[0]=xx[0]/bb[0];
    for(int i=1;i<n;i++){
        double m=1.0/(bb[i]-aa[i]*cc[i-1]);
        cc[i]=cc[i]*m;
        xx[i]=(xx[i]-aa[i]*xx[i-1])*m;
    }
    for(int i=n-2;i>=0;i--){
        xx[i]=xx[i]-cc[i]*xx[i+1];
    }
}
int main(){
    return 0;
}