Ax = b, A = C++의 하위 삼각 행렬

Question

선형 시스템을 해결하는이 함수를 수행하려고하면 A = x 삼각형 행렬, 선형 독립 행렬 및 단 하나의 해가됩니다. 그러나 결과는 항상 0 0 0 0 ... ... 나는 합계의 인쇄했고, 항상뿐만 아니라 0 보여줍니다

#include <iostream> 
using namespace std; 

void solve(int n, float a[][MAX], float b[], float x[]){ 
int i,j; 
float s; 

for(i = 0; i < n; i++){ 
     s = 0; 
     for(j = 0; j < n; j++){ 

      s = s + a[i][j]*x[j]; 
      cout<<s<<endl; 
     } 
    x[i] = (b[i] - s)/a[i][i]; 
} 
} 

Answer 1

void solve(int n, float a[][MAX], float b[], float x[]){ 
    int i,j; 
    float s; 

    for(i = 0; i < n; i++) { 
     s = 0; 
     for(j = 0; j < i; j++) { 
        ^
      s = s + a[ i][ j] * x[ j]; 
     } 
     x[ i] = (b[ i] - s)/a[ i][ i]; 
    } 
} 

BackSubstitution.pdf

compiled example

Answer 2

이 선을 보여

for(j = 0; j < n; j++){ 

이해야 be :

for(j = 0; j < i; j++){ 

그러면 피벗이 항상 0이 아닌 것으로 가정하면 잘 작동합니다.