今日の人気記事

  1. 配列の要素を出力(表示)する

乱数を使って円周率を求めてみよう(モンテカルロ法)

当ページのリンクには広告が含まれています。

スポンサーリンク

ここでは乱数を使って(モンテカルロ法を使って)、円周率を求めてみましょう。

ここではrand関数を使って乱数を生成します。rand関数を使ったことがない方は、まずはこちらをご覧ください。

サンプルプログラム

それではサンプルプログラムをみてみましょう。

/*
 * C言語のサンプルプログラム - Webkaru
 * - 乱数を使って円周率を求めてみよう(モンテカルロ法) -
 */
#include <stdio.h>
#include <stdlib.h>
#include <time.h>

int main(void){

  int i;
  int count = 0;
  int max = 100000;
  double x,y,z,pi;

  /* 乱数の種を初期化 */
  srand(time(NULL));

  /* 0 〜 1 の乱数を生成・出力 */
  for(i=0;i<max;i++){
    x = (double)rand()/RAND_MAX;
    y = (double)rand()/RAND_MAX;
    z = x*x + y*y;
    if(z<1)
      count++;
  }

  /* 円周率を計算・出力 */
  pi = (double)count / max * 4;
  printf("%f\n", pi);

  return 0;
}

実行結果

いくつかの実行結果です。

3.144000
3.147520

3.14までは正確に計算できているようです。試行回数(変数「max」)を増やしたり、より周期の長い乱数を使うとより精度が高くなるでしょう。

サンプルプログラム の人気記事

  1. サンプルプログラム集

  2. 入力した2つの自然数の最大公約数を求める

  3. アルファベットの連続(A to Z)と数字(0 to 9)の連番を出力する

  4. 2進数を10進数に変換する

  5. 入力した自然数の各桁の数字の和を計算する

  6. 配列に格納した数値を昇順ソートする

  7. 入力した摂氏温度を華氏温度に変換する

  8. 自分の名前を出力する

  9. 入力した2つの整数を足し算・引き算する

  10. 級数の和を計算する - 1^2+2^2+3^2+...+n^2

関連記事(一部広告含む)