Tính tích phân bằng phương pháp monte-carlo

[duyquang6]

1- Tạo ngẫu nhiên n tọa độ (x,y) với x thuộc (0.1, 0.2), y thuộc (-200, -50).
2 - Xác định và thống kê các điểm trên nằm trên hay nằm dưới đường cong hàm f(x) = (1 / (log(1 / sqrt(x*x + 1))))
trong khoảng (0.1, 0.2).
3 - Tính tỷ số số điểm nằm trên và n điểm tạo được, nhân với diện tích xy xác định bởi khoảng giá trị trên, ta có giá trị xấp xỉ của tích phân hàm f(x) trong khoảng (0.1, 0.2) (phương pháp Monte Carlo).
Mình thì không biết tí gì về monte carlo cả, đọc trên google nhưng không tìm thấy cái cần tìm. vì thế mong ai đó giúp đỡ mục số 3
Code của mình đây :http://codepad.org/OEtszjW8

[INTP]

tính kết quả ko ra gần đúng à?

cái code kia ko đúng ở mấy điểm:
- y phải là float.
- để sinh số thực ngẫu nhiên thì xài
(double) rand() / RAND_MAX
để lấy ngẫu nhiên trong đoạn [0, 1], rồi chuyển thành đoạn [a, b] dễ dàng:
a + (b-a) * (double) rand() / RAND_MAX



tuy nhiên lấy "tỷ số số điểm nằm trên và n điểm tạo được, nhân với diện tích xy xác định bởi khoảng giá trị trên" thì chưa phải ra tích phân của f(x) trong khoảng (0.1, 0.2) đâu @_@ Còn thiếu 1 mảng hình chữ nhật nằm phía trên nữa...

[prog10]

Chắc ý nói hcn có cạnh là |f(a) - f(b)| và |a-b| đó muh.
Mà cái này dùng h` thang hội tụ nhanh