#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;

#define SEED RAND_MAX-1 //乱数の種
#define R 10 //実験の繰り返し回数

//関数のプロトタイプ
double frand(void) ; //実数乱数
int main(int argc, char *argv[])
{
    int integral ; //積分値
    int n; //試行回数n
    int i, r ; //繰り返し制御変数
    double x, y; //点の座標

    //試行回数
    if(argc<2){//試行係数nの指定がない
    cout<<"使い方 hello (試行回数) \n";
    exit(1);
    }

    if((n=atoi(argv[1]))<=0){
        cout <<"試行回数nが不正です" << "\n";
        exit(1);
    }

    //乱数の初期化
    srand(SEED);

    //積分実験の繰り返し
    for(r=0; r<R; r++){
        integral=0;
        //積分値の計算
        for(i=0; i<n; i++){
            x=frand(); y=frand();
            if((x*x+y*y)<=1)//円の内部の点を数える
            ++integral ;
        }
    //結果の計算
    cout<< "積分値I: "<<(double)integral/n<<"\n4I: "<<(double)integral/n*4<<"\n";
    }
    return 0;
}


/***********************/
/*    frand()関数; 乱数発生*/
/**************************/
double frand(void)
{
    return (double)rand()/RAND_MAX ;
}