Skip to content

Instantly share code, notes, and snippets.

@komasaru
Last active December 18, 2017 08:16
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
Star You must be signed in to star a gist
Save komasaru/576acdb809dfa65ae218523e784fbf0e to your computer and use it in GitHub Desktop.
C++ source code to solve simultaneous equations with Gauss-Jorden method.
/*********************************************
* 連立方程式の解法 ( ガウス・ジョルダン法 )
*********************************************/
#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