C#のfor文

質問

質問内容に合う丁度いいカテゴリが見つからなかったのでここに投稿
しました。相応しくなければご勘弁ください。

C#を勉強中です。

http://ufcpp.net/study/csharp/sp_namespace.html
を読んでいたのですが、理解できないところがあったので質問します。

「エイリアス」の「サンプル」のところのコードについてです。
-----------------------------------------------------
public static double Sin(double x)
{
double xx = -x*x;
double fact = 1;
double sin = x;

for(int i=0; i<100; ++i)
{
fact *= i; ++i; fact *= i; ++i;
x *= xx;
sin += x / fact;
}
return sin;
}//Sin
-----------------------------------------------------

for文の1回めのループの最終行ではfactの値が0になってランタイム
エラーになると思うのですが、実際は6になるのでしょうか？
（プログラムを動かして確認する方法を学ぼうとしている途中なので
こういう聞き方になっています）

0になると考えた理由は以下。

-----------------------------------------------------
fact *= i; ++i; fact *= i; ++i;
-----------------------------------------------------

の最初の「;」の前ではiの値が0なので1*0=0がfactの値になり、
2番めの「;」の前ではiの値が1になり、
3番目の「;」の前では0*1＝0がfactの値になり、
4番目の「;」の前ではiが2になるのではないのでしょうか？

for文の１回目のループが終わった直後に

-----------------------------------------------------
for(int i=0; i<100; ++i)
-----------------------------------------------------

での「++i」によってiは3になるのですよね？

factが6になっているには、

-----------------------------------------------------
fact *= i; ++i; fact *= i; ++i;
-----------------------------------------------------

の最初のfactの代入文でiの値が2、次のfactの代入文でiの値が3
になっていなければいけないように思うのですが、上の考えと合わ
ないためこれ以上読み進めることができなくなりました。

たぶん、初心者ゆえの頓珍漢な質問かと思いますが、ご指導よろしく
お願いします。

回答No.1

masatsan · ベストアンサー

ですね。
数学的には合っているけど、プログラム的には間違っているんです。
sinを展開して計算してるんですが、最初に０の階乗は１に定義されるのでそれをプログラミングしなければなりません。（テーラー展開）
安直に
if (fact==0){fact=1;} を入れるんですね。

お礼

サンプルが変だったんですね。
半分安心しました。
疑ってかかりながら読むのは初心者にはキツイですが納得です。

ありがとうございました。

回答No.4

masatsan

ごめんなさい。何度も。 x*=xx; で良いのですね。　最初に　xx=-x*x やってました。お恥ずかしい。あまり見ずに回答して申し訳ない。

お礼

いえいえ、とんでもありません。
xxという名前があまりよくないということでしょう。
自分でコードを書くときには気をつけたいと思いました。

ありがとうございました。

回答No.3

masatsan

そうですね。
それから　x *= -xx; ですね。

回答No.2

masatsan

No1です。ごめんなさい。私も間違っている。単にif文いれただけではだめです。
初期化しているのに（０！で割る最初の項）iのまわし方が正しくない。
というか色々やり方はあるので考えてください。
少なくともこのままではsinが正しく出てこないばかりか、0割エラーになることは正しいです。
http://akita-nct.jp/yamamoto/lecture/2006/3E/2nd/html/node2.html
ここにsinのテーラー展開が記述されていますので、それのようにプログラミングする。

間違ってごめんなさい。

お礼

ご丁寧にありがとうございます。

多分、ループの部分を
------------------------------------------------
for (int i = 1; i < 101; ++i)
{
++i; fact *= i; ++i; fact *= i;
x *= xx;
sin += x / fact;
}
------------------------------------------------
とでもするんでしょうかね。