Python 몬테카를로 알고리즘 구현

 - 간단한 몬테카를로 알고리즘 구현 ( pi=3.14..... )

 출처 : http://phya.snu.ac.kr/~kclee/lects/lect06/lect06.htm

문제 155. 0에서 1 사이의 숫자를 random으로 생성하시오.

답 : 

# 155 import random x = random.uniform(0,1) y = random.uniform(0,1) print(x) print(y)

문제 156. 위의 코드를 for loop문을 이용해서 10000번 반복시키시오.

답 :

# 156 import random for i in range(10000): x = random.uniform(0,1) y = random.uniform(0,1) print(x) print(y)

문제 157. x^2 + y^2 <= 1 의  안에 들어가는 점을 이용해서 pi를 구해보시오.

답 :

# 157 import random j=0 cnt=0 for i in range(1000000): j=j+1 x = random.uniform(0,1) y = random.uniform(0,1) if x**2 + y**2 <= 1: cnt = cnt+1 print(cnt/j*4)

문제 158. 정사각형 안에 찍을 점의 갯수를 물어보게끔 코드를 수정하시오.

답 :

# 158 a = int(input('정사각형 안에 들어갈 점의 갯수를 입력하세요 : ')) import random j=0 cnt=0 for i in range(a): j=j+1 x = random.uniform(0,1) y = random.uniform(0,1) if x**2 + y**2 <= 1: cnt = cnt+1 print('파이(pi)는 %f' %(cnt/j*4))