このディレクトリの索引
http://hibari.2ch.net/test/read.cgi/tech/1291471791/798
#  [1] 授業単元:数値数理解析 
#  [2] 問題文(含コード&リンク): 
#  リンク先の画像に問題文を載せています。 
#  ttp://www.dotup.org/uploda/www.dotup.org1315077.jpg 
#  
#  三次元空間で
#  x > 0, y > 0, 0 < z < 1, x^2 + y^2 < (1-z)^2
#  を満たす領域は半径1,高さ1の直円錐の四分の一と同じになる.
#  

'x > 0, y > 0, 0 < z < 1, x^2 + y^2 < (1-z)^2を満たす領域は半径1,高さ1の直円錐の四分の一と同じになることをモンテカルロ法で示す'(_試行回数,_満足回数) :-
        count((
                    for(1,N,_試行回数),
                    X is (random mod 10001) / 10000,
                    Y is (random mod 10001) / 10000,
                    Z is (random mod 10001) / 10000,
                    'x > 0, y > 0, 0 < z < 1, x^2 + y^2 < (1-z)^2を満たす'(X,Y,Z,U)),
                _満足回数).

'x > 0, y > 0, 0 < z < 1, x^2 + y^2 < (1-z)^2を満たす'(X,Y,Z,1) :-
        (X^2 + Y^2) < ((1 - Z) ^ 2),!.