Last active
December 18, 2017 08:16
Star
You must be signed in to star a gist
C++ source code to solve simultaneous equations with Gauss-Jorden method.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
/********************************************* | |
* 連立方程式の解法 ( ガウス・ジョルダン法 ) | |
*********************************************/ | |
#include <iostream> // for cout | |
#include <stdio.h> // for printf() | |
// 元の数定義 | |
#define N 3 // 4 | |
using namespace std; | |
/* | |
* 計算クラス | |
*/ | |
class Calc | |
{ | |
double a[N][N + 1]; // 係数用配列 | |
// 各種変数 | |
double p, d; // ピボット係数、ピボット行x係数 | |
int i, j, k; // LOOP インデックス | |
public: | |
// 連立方程式を解く(ガウス・ジョルダン法) | |
void calcGaussJorden(); | |
}; | |
/* | |
* 連立方程式を解く(ガウス・ジョルダン法) | |
*/ | |
void Calc::calcGaussJorden() | |
{ | |
// 係数 | |
static double a[N][N + 1] = { | |
{ 2.0, -3.0, 1.0, 5.0}, | |
{ 1.0, 1.0, -1.0, 2.0}, | |
{ 3.0, 5.0, -7.0, 0.0} | |
//{ 1.0, -2.0, 3.0, -4.0, 5.0}, | |
//{-2.0, 5.0, 8.0, -3.0, 9.0}, | |
//{ 5.0, 4.0, 7.0, 1.0, -1.0}, | |
//{ 9.0, 7.0, 3.0, 5.0, 4.0} | |
}; | |
// 元の連立方程式をコンソール出力 | |
for (i = 0; i < N; i++) { | |
for (j = 0; j < N; j++) | |
printf("%+fx%d ", a[i][j], j + 1); | |
printf("= %+f\n", a[i][N]); | |
} | |
for (k = 0; k < N; k++) { | |
// ピボット係数 | |
p = a[k][k]; | |
// ピボット行を p で除算 | |
for (j = k; j < N + 1; j++) | |
a[k][j] /= p; | |
// ピボット列の掃き出し | |
for (i = 0; i < N; i++) { | |
if (i != k) { | |
d = a[i][k]; | |
for (j = k; j < N + 1; j++) | |
a[i][j] -= d * a[k][j]; | |
} | |
} | |
} | |
// 結果出力 | |
for (k = 0; k < N; k++) | |
printf("x%d = %f\n", k + 1, a[k][N]); | |
} | |
/* | |
* メイン処理 | |
*/ | |
int main() | |
{ | |
try | |
{ | |
// 計算クラスインスタンス化 | |
Calc objCalc; | |
// 連立方程式を解く(ガウス・ジョルダン法) | |
objCalc.calcGaussJorden(); | |
} | |
catch (...) { | |
cout << "例外発生!" << endl; | |
return -1; | |
} | |
// 正常終了 | |
return 0; | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment