行列式

25-181823.cpp

#include<iostream>
#include<iomanip>

int main(int argc, char *argv[]){
  int n;
  double input;
  std::cout << "n次の正方行列? n:";
  std::cin >> n;
  double a[n][n];
  for(int i=0;i<n;i++){
    for(int j=0;j<n;j++){
      std::cout << i+1 << "行" << j+1 << "列目:";
      std::cin >> input;
      a[i][j] = input;
    }
  }

  double buf;
  for(int i=0;i<n;i++){
    for(int j=0;j<n;j++){
      buf = a[j][i]/a[i][i];
      if(i<j){
	for(int k=0;k<n;k++){
	  a[j][k] = a[j][k] - a[i][k]*buf;
	}
      }
    }
  }

  std::cout << std::fixed << std::setprecision(4);

  double det=1.0;
  for(int i=0;i<n;i++){
    det = det*a[i][i];
  }

  std::cout << "行列式:" << det << std::endl;

  return 0;
}