以下のプログラムを作ったんですが実行したときに答えがnanになります。どなたか間違いがわかる方よろしくおねがします。
const double e =1e-10;
enum {N=100};
const double x0 =1;
const double y0 =1;
double f(const double x,double y) {return x*x-y*y;}
double g(const double x,double y) {return 6*x*y-5;}
double dfdx(const double x,double y){return 2*x;}
double dfdy(const double x,double y){return -2*y;}
double dgdx(const double x,double y){return 6*y;}
double dgdy(const double x,double y){return 6*x;}
double det(double x,double y);
main()
{
double x[N+1],y[N+1];
int i;
double det2;
x[0]=x0;
y[0]=y0;
for(i=0;i<=20;i++)
{
det2=det(x[i],y[i]);
x[i+1]=x[i]-(dgdy(x[i],y[i])*f(x[i],y[i]))-(dfdy(x[i],y[i])*g(x[i],y[i]))/det2;
y[i+1]=y[i]-(dgdx(x[i],y[i])*f(x[i],y[i]))-(dfdx(x[i],y[i])*g(x[i],y[i]))/det2;
printf("x = %lf\n",x[i+1]);
printf("y = %lf\n",y[i+1]);
printf("\n");
}
exit(0);
}
double det(double a,double b)
{
double ans;
ans=1/(dfdx(a,b)*dgdy(a,b)-dfdy(a,b)*dgdx(a,b));
return ans;
}
No.3ベストアンサー
- 回答日時:
ANo.1 = 2 です。
バラバラ書いてごめんなさい。
> x[i+1]=x[i]-(dgdy(x[i],y[i])*f(x[i],y[i]))-(dfdy(x[i],y[i])*g(x[i],y[i]))/det2;
これはいいけど、
> y[i+1]=y[i]-(dgdx(x[i],y[i])*f(x[i],y[i]))-(dfdx(x[i],y[i])*g(x[i],y[i]))/det2;
これ、プラスマイナスの付け方間違っていませんか?
というのは、ニュートン法を使うためにヤコビ行列の逆行列を使いますよね? まずは、逆行列が存在しないケース(すなわち行列式=0の場合)が抜けていませんか、というのがANo.1,2 の指摘。今回の指摘は、逆行列に f(x,y), g(x,y) を書けるのにプラスマイナスが間違っちゃってませんか、という指摘です。
プログラム云々の前に手で数式を展開して、最初の1、2回目程度手計算した結果とプログラムが吐き出す答えをつき合わしてみてはいかがでしょうか。
No.2
- 回答日時:
というか、その前に、
> double det(double a,double b)
という関数はきっと行列式を計算したいのだと思うけど、
> ans=1/(dfdx(a,b)*dgdy(a,b)-dfdy(a,b)*dgdx(a,b));
って、
det = dfdx(a,b)*dgdy(a,b)-dfdy(a,b)*dgdx(a,b);
じゃないですか?
No.1
- 回答日時:
> 実行したときに答えがnanになります。
というのは、ループの何回目からですか?
そのときのx,yの値は?
ans=1/(dfdx(a,b)*dgdy(a,b)-dfdy(a,b)*dgdx(a,b));
で、分母が0になる可能性は?
まずはここにエラー処理を入れてみませんか?
それから、
> const double e =1e-10;
というのはどこかで使ってますか?
> enum {N=100};
としながら、forループが20回なのはなぜですか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# ある線が円の範囲に入っているかの計算 1 2022/12/07 16:14
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- Excel(エクセル) 【マクロ】スクショ印刷がうまく動かない件 5 2022/12/06 17:37
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- FX・外国為替取引 mql4のコンパイルエラー箇所の修正お願いします。 1 2023/03/15 16:14
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
c言語で、繰り返し文の中で、0....
-
プログラムでの数字につく”f”の...
-
C言語のプログラムで#include<m...
-
C言語で直角三角形の斜辺を求め...
-
sin(x)の近似について
-
doubleの変数にintとintの割り...
-
Cで3乗根を求める方法
-
C言語初心者 構造体 課題について
-
至急です! マクロ定義で #defi...
-
float型とdouble型の変数の違い...
-
C言語を実行すると-infが出てき...
-
2分法で方程式の複数の解を自...
-
MATLABで画像のヒストグラムを...
-
浮動小数点数が表示されないん...
-
関数におけるif文とreturn文に...
-
2次方程式の解を求めるプログ...
-
double型とint型で三分の一乗の...
-
floating point not loadedとは?
-
C言語 関数プロトタイプ宣言の...
-
C言語で表記についの質問です
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
C 開放してるのにエラー(doubl...
-
Cで3乗根を求める方法
-
float型とdouble型の変数の違い...
-
至急です! マクロ定義で #defi...
-
C言語の型による処理速度の違い
-
int とdoubleの比較
-
関数におけるif文とreturn文に...
-
C言語初心者 構造体 課題について
-
c言語のコンパイルエラー canno...
-
C言語 関数プロトタイプ宣言の...
-
C言語を実行すると-infが出てき...
-
float?数字の後にLがつくもの
-
数値を指数部と仮数部に分離したい
-
difftime()について
-
浮動小数点数が表示されないん...
-
たくさんの数の平均を求める方...
-
DWORDの警告
-
-1.#IND00と出てしまうのですが...
おすすめ情報