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]); }
}
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。