function v = integral_Monte_Carlo ( a, n )
% モンテカルロ法により積分を計算
% 積分 \int_{0}^{a} 1/(1+x^2) を計算
% (x,y)をn点発生させる
    m = 0;
    for k=1:n
        x = a*rand;
        y = rand;
        if (y <= 1/(1+x*x))
            m = m+1;
        end
    end
    v = m/n*a;
end
