蒙特卡罗积分(Scheme语言实现)

蒙特卡罗积分(Scheme语言实现)

      具体见<计算机程序的构造和解释>第三章习题3.5。

      程序如下：

      

1. (define (estimate-pi trails)  
2.   (sqrt (/ 6 (monte-carlo trails experiment))))  
3.   
4. (define (random-in-range low high)  
5.   (let ((range (- high low)))  
6.     (+ low (random (+ 1 range)))));
7. (define (monte-carlo trails experiment)  
8.   (define (iter remain res)  
9.     (cond ((= 0  remain) (/ res trails))  
10.           ((experiment)  
11.            (iter (- remain 1) (+ 1 res)))  
12.           (else (iter (- remain 1) res))))  
13.   (iter trails 0))  
14.   
15. (define (fractional trails x1 x2 y1 y2 p)  
16.   (define (test-func)   
17.      (p x1 x2 y1 y2))  
18.   (* (* (- x2 x1) (- y2 y1)) (monte-carlo trails test-func)))  
19.   
20. (define (test-func x1 x2 y1 y2)  
21.   (< (+ (square (random-in-range x1 x2))  
22.         (square (random-in-range y1 y2)))  
23.      1))  
24.     
25. (define (p x1 x2 y1 y2)  
26.   (< (+ (square (random-in-range x1 x2))  
27.         (square (random-in-range y1 y2)))  
28.      1))  
29.   
30. (define the-pi   
31.   (lambda (trails x1 x2 y1 y2 p)  
32.     (fractional trails x1 x2 y1 y2 p)  
33.        ))  
34.   
35. (define (square x) (* x x))