vc++(2010)でガウスの消去法を使って連立方程式を解く、というプログラムを組みました。
正直中身自体をきちんと理解していないので、
組みながら理解しようと思って組んだのですが、
結果に-1.#IND00と出てしまいました。
これは何なのでしょうか?
下記がそのプログラムです。
#include<stdio.h>
#include<math.h>
#define N 3
int main(void){
double a[N][N+1] = {{2,3,4},
{3,-2,5},
{5,4,-7}};
double b[1][N+1] = {20,14,-8};
double x[N];
int i = 0;
int j = 0;
int k = 0;
int l = 0;
int pivot = 0;
double p = 0;
double q = 0;
double m = 0;
for(i = 0; i < N; i++){
x[i] = 0;
}
for(i = 0; i < N; i++){
pivot = i;
for(l = i; l < N; l++){
if(fabs(a[l][i]) > m){
m = fabs(a[l][i]);
pivot = l;
}
}
if(pivot != i){
for(j = 0; j < N; j++){
b[0][j] = a[i][j];
a[i][j] = a[pivot][j];
a[pivot][j] = b[0][j];
}
}
}
for(k = 0; k < N; k++){
a[k][j] = a[k][j] / p;
a[k][k] = 1;
for(j = k; j < N; j++){
a[k][j] = a[k][j] / p;
}
for(i = k+1; i < N; i++){
q = a[i][k];
for(j = k+1; j < N; j++){
a[i][j] = a[i][j] - q*a[k][j];
}
a[i][k] = 0;
}
}
for(i = N-1; i >= 0; i--){
x[i] = a[i][N];
for(j = N-1; j > i; j--){
x[i] = x[i] - a[i][j] * x[j];
}
}
//for(i = 0; i < N; i++){
//for(j = 0; j < N+1; j++){
//printf("%.1f", a[i][j]);
//}
//printf("\n");
//}
//printf("解\n");
//for(i = 0; i < N; i++){
//printf("%f\n", x[i]);
//}
return 0;
}
最後のコメントにしてある行は解を表す時と行列を表す時で使い分けているので、
実際はどちらかを外して使用しています。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
初めて見た映画を教えてください!
初めて見た映画を覚えていますか?
-
歳とったな〜〜と思ったことは?
歳とったな〜〜〜、老いたな〜〜と思った具体的な瞬間はありますか?
-
みんなの【マイ・ベスト積読2024】を教えてください。
積読、ついついしちゃいませんか?そこでみなさんの 「2024年に買ったベスト積読」を聞きたいです。
-
AIツールの活用方法を教えて
みなさんは普段どのような場面でAIツール(ChatGPTなど)を活用していますか?
-
今から楽しみな予定はありますか?
いよいよ2025年が始まりました。皆さんには、今から楽しみにしている予定はありますか?
-
C++で誤差関数erf(x)を使いたいのですが...
C言語・C++・C#
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・あなたの人生で一番ピンチに陥った瞬間は?
- ・初めて見た映画を教えてください!
- ・今の日本に期待することはなんですか?
- ・【大喜利】【投稿~1/31】『寿司』がテーマの本のタイトル
- ・集中するためにやっていること
- ・テレビやラジオに出たことがある人、いますか?
- ・【お題】斜め上を行くスキー場にありがちなこと
- ・人生でいちばんスベッた瞬間
- ・コーピングについて教えてください
- ・あなたの「プチ贅沢」はなんですか?
- ・コンビニでおにぎりを買うときのスタメンはどの具?
- ・おすすめの美術館・博物館、教えてください!
- ・【お題】大変な警告
- ・【大喜利】【投稿~1/20】 追い込まれた犯人が咄嗟に言った一言とは?
- ・洋服何着持ってますか?
- ・みんなの【マイ・ベスト積読2024】を教えてください。
- ・「これいらなくない?」という慣習、教えてください
- ・今から楽しみな予定はありますか?
- ・AIツールの活用方法を教えて
- ・最強の防寒、あったか術を教えてください!
- ・【大喜利】【投稿~1/9】 忍者がやってるYouTubeが炎上してしまった理由
- ・歳とったな〜〜と思ったことは?
- ・モテ期を経験した方いらっしゃいますか?
- ・好きな人を振り向かせるためにしたこと
- ・スマホに会話を聞かれているな!?と思ったことありますか?
- ・それもChatGPT!?と驚いた使用方法を教えてください
- ・見学に行くとしたら【天国】と【地獄】どっち?
- ・これまでで一番「情けなかったとき」はいつですか?
- ・この人頭いいなと思ったエピソード
- ・あなたの「必」の書き順を教えてください
- ・14歳の自分に衝撃の事実を告げてください
- ・人生最悪の忘れ物
- ・あなたの習慣について教えてください!!
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
ラグランジュの補間法のCプログ...
-
float型とdouble型の変数の違い...
-
c言語のプログラミングについて...
-
バイナリから実数を取得したい...
-
listに構造体を格納
-
C言語を実行すると-infが出てき...
-
2次方程式の解を求めるプログ...
-
プログラムでの数字につく”f”の...
-
C言語のpow関数の不具合
-
C#イベント中の戻り値の設定の...
-
浮動小数点の定数
-
c言語の構造体エラーについて
-
doubleは常に%lfとするべきなのか
-
C言語でsqrt(a^2+b^2)のテーブ...
-
「Aに対するBの割合」と「Aに対...
-
信頼区間の1.96や1.65ってどこ...
-
C言語での引数の省略方法
-
【gcc・cygwin】multiple defin...
-
break文でループを一気に抜ける...
-
数学のデータの相関の散布図で...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
C言語を実行すると-infが出てき...
-
C 開放してるのにエラー(doubl...
-
float型とdouble型の変数の違い...
-
c言語で、繰り返し文の中で、0....
-
至急です! マクロ定義で #defi...
-
C言語 関数プロトタイプ宣言の...
-
C言語の型による処理速度の違い
-
2次方程式の解を求めるプログ...
-
関数におけるif文とreturn文に...
-
int とdoubleの比較
-
C++で外積
-
C言語で台形公式を使った二重積...
-
指数の表示
-
C言語のpow関数の不具合
-
(C,C++言語)関数の引数は自動キ...
-
C言語のプログラムで#include<m...
-
数値を指数部と仮数部に分離したい
-
プログラミングでのテイラー展開
おすすめ情報