Вычисление числа Пи методом монте-карло. В чем ошибка?

@goblin_l · Регистрация: 12.04.2011

Студворк — интернет-сервис помощи студентам

Требуется создать программку для вычисления числа Пи методом монте-карло.
Суть заключается в следующем: генерируются 2 случайных числа А и В, равномерно распределенных в интервале (-1,1). Если выполняется условие А^2+B^2 <= 1, то точка с координатами (А,В) попадает внутрь круга радиусом 1 с центром в начале координат, если нет, то не попадает. Если при N испытаниях получится n попаданий в круг, то приближенно число Пи можно вычислить по формуле Pi = 4*n/N.
Была сделана вот такая программка:

Visual Basic
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Dim n, m As Integer, s, pi As Double
 
Private Sub Command1_Click()
s = Text1.Text
pi = 0
Do Until m = s
x = Rnd * 2 - 1
y = Rnd * 2 - 1
If (x * x + y * y) <= 1 Then
Picture1.PSet (x, y), vbRed
n = n + 1
Else
Picture1.PSet (x, y), vbBlue
m = m + 1
End If
Loop
pi = n/m
Text2.Text = pi
 
End Sub

Из формулу была убрана 4, с ней число Пи равно 15, ...
По идее если сделать N = 10000 (s = 10000), то число Пи будет равно 3,14...
В данном случае число Пи равно 3,664. В чем ошибка?

@gaw · 13.04.2011, 17:57

Visual Basic
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
Dim n&, m&, s&, pi As Double
Picture1.Scale (-1.5, 1.5)-(1.5, -1.5)
Picture1.Cls
s = Text1.Text
m = 0: n = 0
pi = 0
Do Until m = s
x = Rnd * 2 - 1
y = Rnd * 2 - 1
If (x ^ 2 + y ^ 2) <= 1 Then
Picture1.PSet (x, y), vbRed
n = n + 1
Else
Picture1.PSet (x, y), vbBlue
End If
m = m + 1
Loop
pi = 4 * n / m
Text2.Text = pi

@goblin_l · 13.04.2011, 20:12 **[ТС]**

А вот в чем ошибка

Спасибо, так и сделаю.

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .

@goblin_l 0 / 0 / 0 Регистрация: 12.04.2011 Сообщений: 11
	13.04.2011, 20:12 [ТС]
	А вот в чем ошибка Спасибо, так и сделаю. 0