。まったくわかりません。

Question

log n

発生している問題・エラーメッセージ

n = 1000のとき
r(n) = -5.907755e+000
n = 10000のとき
r(n) = -8.210340e+000
n = 100000のとき
r(n

for(j = 1; j <= M; j++){
b += 1/j;
}
printf("n = 10000のとき\n");
printf("r(n) = %e\n",b - log(M));

for(k = 1; k <= O; k++){
c += 1/k;
}
printf("n = 100000のとき\n");
printf("r(n) = %e\n",c - log(O));

}

### 試したこと
これはまだ途中ですが、出てきた計算結果が近似値と全く違います。



### 補足情報（FW/ツールのバージョンなど）
私のコードの間違いと解答例を教えてください。

Comments

[yumetodo]

teratailにはソースコードをハイライトする機能があるのでそれを使ってください。さっぱり読めません。あとインデント。

[tiitoi]

コードはインデックスが崩れてしまいますので、 ```コード``` と囲ってください。質問はあとからでも編集できます。

[nishimuu]

すいません。初心者でテラテイルの使い方がまだよくわかっていなくて・・・

Answer 1

for(i = 1; i<= N; i++){
a += 1/i;
}

C言語では、int / int は int になります。
つまり、2以上のiについて 1/i は常に0です。

演算対象の片方をdouble型にして計算してください。

for(i = 1; i<= N; i++){
    a += 1.0/i;
}

質問の評価について

瞬間風速では -5 まで評価が落ちてしまっていました。
内容と言うより、単に質問の書き方が良くなかったように思います。

次の点に注意されると良いかと。

• タイトルを適切に

タイトルは質問内容を的確に表すものにしてください。
今回のタイトルには必要な情報が一切含まれていないだけでなく、丸投げな感も与えます。

• コードブロックを使う

見づらい質問は、回答者に配慮しない身勝手な印象を与えます。
(指摘後に修正されているので、単に知らなかっただけかと思いますが)

• 地雷ワード

『私のコードの間違いと解答例を教えてください。』
宿題なのに解答例を求めたらダメだろ、と多くの人が感じるかと。

必要以上に低評価を食らっても損しかないので、次回以降気を付けてみてください。
ヘルプに目を通しておくとなお良いでしょう。teratailヘルプ - 質問するときのヒント

コードの書き方について

teratailには、コードを見やすく表示する機能があります。
質問編集画面を開き、コードを選択した状態で<code>ボタンを押してください。

Comments

[nishimuu]

なぜint/intだとおかしな数字になるのでしょうか？

[LouiS0616]

int/intはintを返すと言語仕様で決まっているからです。
正数同士の除算の場合は、小数点以下切り捨てになるはずです。

[nishimuu]

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

[LouiS0616]

いや、『宿題であることを伏せろ』という意味では無かったのですが...
解答(狭義にはそのまま動くコード) を求めるのがまずいんであって、宿題のアドバイスを常識的な範囲で求めるのは問題ないんじゃないでしょうか。

Answer 2

簡潔に書くなら関数化したほうがよいですね。

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

/**
   @brief オイラー定数を num_terms 項まで計算する。
   @param n 計算する項数
   @return オイラー定数
 */
double calcEulerConstant(int num_terms)
{
    assert(n >= 1);

    double sum = 0;
    for (int i = 1; i <= n; ++i) {
        sum += 1. / i;
    }

    return sum - log(n);
}

int main()
{
    int num_terms[] = {1000, 10000, 100000, 1000000, 10000000, 100000000};
    for (int i = 0; i < sizeof(num_terms) / sizeof(num_terms[0]); ++i)
        printf("r(%d) = %.15f\n", max[i], calcEulerConstant(max[i]));
}

r(1000) = 0.577715581568206 // 3桁まで一致
r(10000) = 0.577265664068165 // 4桁まで一致
r(100000) = 0.577220664893106 // 4桁まで一致
r(1000000) = 0.577216164900715 // 5桁まで一致
r(10000000) = 0.577215714898951 // 6桁まで一致
r(100000000) = 0.577215669900188  // 8桁まで一致

Wolframe Alpha の答え

0.577215664901532860606512090082402431042159335939923598805