![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
指数関数のマクローリン展開のプログラムです。
x=1の場合はしっかりと値が収束するのですが、それ以外ではうまくいきません。
以下のように組んだのですが、どこがまちがっているのでしょうか?
1.n!を求める。
2.x^nを求める。
3.これらを割る。
4.exに足していくという手順でやっているつもりです。
よろしくおねがいします。
______________________
#include<stdio.h>
main()
{
double a,i,j,k,l,m,n,x,ex;
printf("x=");
scanf("%lf",&x);
n=1;
ex=1;
m=x;
for(i=1;i<100;i++){
m=x;
n=n*i;
for(k=1;k<i;k++){
m=m*m;
}
a=m/n;
if(a<0.000001){
return 0;
}else{
m=x;
ex=ex+a;
}
printf("%lf ",ex);
}
return()
}
No.2ベストアンサー
- 回答日時:
分子・分母を個別に計算しているので、それぞれオーバーフロー(doubleで表現できる値の範囲を越える)してますね。
途中の m とか n の値を表示させてみるとわかると思います。それ以外にもおかしな点はありますが。
もとのプログラムをあまり変えない範囲で書き直すと、
#include <stdio.h>
main()
{
double x,a,ex,i;
printf("x=");
scanf("%lf",&x);
a=1;
ex=1;
for(i=1;i<100;i++){
a=a*x/i;
if(a<0.000001) break;
ex=ex+a;
printf("%lf\n",ex);
}
return 0;
}
数値計算は奥深いので、もっといいやり方があるかもしれません。
検証してみました。
これほど簡単にできるものなのですね。
考え方として、分子分母別で計算した方が考えやすいかと思ってやっていました。
これを参考にまた作りなおしてみようと思います。
ありがとうございました。
No.1
- 回答日時:
ご提示のソースは
>return()
ここでコンパイルエラーが出るため、論評できかねます。
お手元にある、x=1以外の場合でうまく動かないというソースを、
そっくりそのまま貼り付けていただけますか?
すいませんでした
その部分は書き間違いです・・・
こちらでコンパイルは正常にできています
コピペができないのでアレなんですけども・・・
x=5の場合
6.00000000
18.5000000
122.666666
16398.70833・・・
1271582153.・・・・・・
これ以降はものすごく大きな数字になります。
そして Inf Inf Inf・・・ と表示されています。
おそらくxのn乗のあたりがおかしいからこうなるのかと思っているのですが、朝からずっといろいろやってみていてこのありさまです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・「I love you」 をかっこよく翻訳してみてください
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・昔のあなたへのアドバイス
- ・かっこよく答えてください!!
- ・あなたが好きな本屋さんを教えてください
- ・スタッフと宿泊客が全員斜め上を行くホテルのレビュー
- ・【大喜利】【投稿~8/27】 こんなガソリンスタンド二度と来るか!なぜそう思った?
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・【お題】動物のキャッチフレーズ
- ・【お題】甲子園での思い出の残し方
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・「それ、メッセージ花火でわざわざ伝えること?」
- ・自分用のお土産
- ・人生で一番お金がなかったとき
- ・一番好きなみそ汁の具材は?
- ・泣きながら食べたご飯の思い出
- ・ちょっと先の未来クイズ第1問
- ・ゴリラ向け動画サイト「ウホウホ動画」にありがちなこと
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
C++で等差数列の和に関する問題...
-
-1.#IND00と出てしまうのですが...
-
プログラムでの数字につく”f”の...
-
2分法で方程式の複数の解を自...
-
doubleの変数にintとintの割り...
-
プログラミングについての質問
-
c言語で、繰り返し文の中で、0....
-
C言語 正弦関数の級数展開のプ...
-
int とdoubleの比較
-
浮動小数点の定数
-
fgets()による繰り返し入力(c言語)
-
C言語を実行すると-infが出てき...
-
C++の質問です
-
インデックスが配列の境界外です.
-
2次方程式の解
-
C言語で内積、、、わかりません。
-
指数形式で入力するには
-
ガウスの消去法を用いたC言語...
-
C言語の乱数生成について。 0〜...
-
インデックスが配列の境界外で...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
c言語で、繰り返し文の中で、0....
-
C 開放してるのにエラー(doubl...
-
C言語を実行すると-infが出てき...
-
float型とdouble型の変数の違い...
-
至急です! マクロ定義で #defi...
-
C言語の型による処理速度の違い
-
関数におけるif文とreturn文に...
-
C言語 関数プロトタイプ宣言の...
-
浮動小数点の定数
-
int とdoubleの比較
-
doubleは常に%lfとするべきなのか
-
difftime()について
-
C言語でdouble型の小数点の引き...
-
たくさんの数の平均を求める方...
-
斜辺と角度て底辺と高さを出したい
-
C言語で台形公式を使った二重積...
-
float?数字の後にLがつくもの
-
c言語のコンパイルエラー canno...
おすすめ情報