合同式法で乱数を10万個生成し,分散を表示する問題

Question

###前提・実現したいこと
乱数を10万個生成し,その分散を求めたい.

###発生している問題・エラーメッセージ
今発生している問題は,分散値を表示させると,"inf"と表示され,答えが出ない状況です.

1.501441
inf

###該当のソースコード

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

int main(void)
{

  double average;
  double bunsan1=0;
  double bunsan2=0;
  double bunsan_new;
  int a=1229;
  int c=351750;
  int m=1664501;

  double x=0;
  double x1;
  int r=1;
  int i=0;

  for(i=0;i<100000;i++){
    r=(a*r+c)%m;
    x1=(double)r/m;
    x=x+x1;
    average=x/i+1;
    bunsan_new=(x1-average)*(x1-average);
    bunsan1=bunsan_new+bunsan1;
    bunsan2=bunsan1/i+1;
      }

  printf("%f\n",average);
  printf("%f\n",bunsan2);

  return 0;
}

###補足情報(言語/FW/ツール等のバージョンなど)
言語はcを使いました.

Answer 1

+の優先度は/より下です。
次のように置き換えてください。

average=x/i+1; → average=x/(i+1);
bunsan2=bunsan1/i+1; → bunsan2=bunsan1/(i+1);

質問者様のプログラムでは、i=0による割り算が発生しています。

Comments

[RYUi]

初歩的なミスですみませんでした.ご指摘ありがとうございます.