c++でのテイラー展開

Question

c++を使ってe^x=1+x/1!+x^2/2!+…=Σ[n=0→∞]x^n/n!を10^(-10)の精度でx=0.2から0.2きざみでx=10.0まで計算したくてプログラムを以下のように書きました。
#include<stdio.h>
#include<math.h>
int main(void){
int n; (項数)
int k; (項数カウント)
double x; (e^xのx)
double xk; (x^n/n!)
double e; (テイラー展開の部分和)
double ex; (exp(x))
double er; (exp(x)との差)

e = xk = 1,n=0,ex=1.0,er=1.0;
for(x=0.2;x<=10.0;x+=0.2)｛
for(k=1;k<=n;k++){
xk *= x/k;
e += xk;
}
｝
ex = exp(x);
er = fabs(e-ex);
}
このように書いてコンパイル自体は通ったのですが、思うような結果が出なくて困っている状態です。
(項数)などは、実際にプログラムには書いてないです。
どこをどのように直せば思うような結果が得られますか？

Comments

[y_waiwai]

あなたの思うような結果と、実際に出ている結果を提示しましょう
また、
このままではコードが見づらいので、質門を編集し、<code>ボタンで、出てくる’’’の枠の中にコードを貼り付けてください

Answer 1

#include <stdio.h>
#include <math.h>

#define F "%-11.10g"

double f(double x)
{
	double y = 1, t = 1;
	for (int n = 1; t > y * 1e-10; n++) y += t *= x / n; 
	return y;
}

int main(void)
{
	for (int i = 2; i <= 100; i += 2) {
		double x = i * 0.1;
		printf("x =%4.1f: f(x) = " F ", exp(x) = " F "\n", x, f(x), exp(x));
	}
}

・#define F "%-13.12g" にすると、精度が 10桁だと分かるでしょう。

Answer 2

e = xk = 1,n=0,ex=1.0,er=1.0;

各xの値毎にeを計算したいのでしょうから
e = xk = 1
これを，xに関するforの中に入れる必要があるのでは．

あと，n=0だとkに関するforが回りませんよ．

誤差erを求める部分もx毎に行う事柄だと思います．

「あるxの値に対する e^xの近似値を求める」という関数でも作れば
見通しが良くなってミスが減るんじゃないでしょうか．