KTUG마당은 KTUG를 방문하는 모든 이용자가 대화를 나누고 소식을 전하는 곳입니다.
- 로그인 없이 자유롭게 글을 읽고 쓸 수 있는 철학은 처음과 같이 계속됩니다.
- Team Blog의 글을 이곳 게시판의 "정보글"로 모았습니다. Team blog는 기고자가 올린 글에 질문을 받는 부담을 줄이기 위하여 댓글을 허용하지 않았습니다. 그러나 이곳 게시판으로 모으면서 댓글을 달 수 있습니다. 게시물을 작성하실 때 댓글을 원하지 않으시면 댓글을 허용하시지 않으시기를 바랍니다. 또한 불필요한 소모성 댓글을 달지 않도록 주의하여 주시기를 바랍니다.
- TeX과 관련된 질문이나 답변은 QnA 마당을 이용하십시오. TeX과 관련된 질문은 지웁니다
- MathJax를 이용한 수식조판을 사용하실 수 있습니다. 여기를 참조하세요.
- 스팸 글을 막기 위하여 짧은 시간 내에 다시 글이 등록되는 IP를 막거나, 광고 글을 막기 위하여 금지어로 .com, .net 등을 설정하고 있습니다. 다소간의 불편함이 있으시더라도 양해 바랍니다.
- 금지어에서 stackexchange, stackoverflow, ctan, overleaf, , github, google.com, gmail.com, .org, .io, sil.org, wiki.com, tistory.com등은 해제하였습니다.
- 사용하는 편집기는 CKeditor입니다. 편집기에서 [enter]를 누르면 <p> 태그가 들어가고, 문단으로 생각하고 한줄을 비웁니다. 글줄만 바꾸려면 shift-enter 를 누르시면 <BR>가 들어가므로 용도에 맞게 나누어 쓸 수 있습니다.
정보글 몬테카를로 기법으로 파이 구하기: 정수 연산 버전
2011.06.15 00:34
\newcount\N \N=1000000 % 총시행횟수. 10^n \newcount\n \n=0 % 현재시행횟수 \newcount\hit \hit=0 % Hit \newcount\rr \rr\numexpr10000*10000\relax % 상수 r^2 \newcount\p \newcount\q % 난수 담을 변수 \def\rand#1{#1=\pdfuniformdeviate10001\relax} % 난수 발생기 \def\dmark#1{#1.} % 소수점 찍기 \loop \rand\p \rand\q \ifnum\numexpr\p*\p+\q*\q\relax>\rr\else\advance\hit by1\fi \ifnum\n<\N\advance\n by1 \repeat $$ \pi \approx \expandafter\dmark\the\numexpr4*\hit\relax $$ \bye
아래 도은이아빠 글을 보고 자원 소모가 많은 소수점 연산을 안 써도 얼추 가능할 거 같아서 간단히 작성해봤다.
첫줄의 \N
은 10의 제곱수라야 한다. 이 값을 증가시키면 파이 값이
더 정확해진다.
위의 예처럼 1백만 번까지는 그럭저럭 즐겁게 컴파일할 수 있지만 0을 하나 더 붙여 1천만 번을 돌리면 기다리기가 사뭇 지겹다.
굳이 레이텍으로 만들 이유가 없기에 플레인텍으로 작성했다. pdftex 으로
컴파일하면 된다.
댓글 0
번호 | 제목 | 글쓴이 | 날짜 | 조회 수 |
---|---|---|---|---|
59 | texindy-ko.bat | karnes | 2011.07.30 | 4581 |
58 | 짧게 알아보는 windows에서 TeX Live 2011 + ko.TeX 설치하기 | 샘처럼 | 2011.07.30 | 13978 |
57 | xindy 한글 처리 실험 | DohyunKim | 2011.07.29 | 4534 |
56 | 상용로그표 : 루아텍 버전 | DohyunKim | 2011.07.14 | 5234 |
55 | 상용로그표 | karnes | 2011.07.13 | 5149 |
54 | Right to Left | DohyunKim | 2011.07.12 | 4525 |
53 | windows에서 jbig2을 이용한 pdf만들기 | 샘처럼 | 2011.07.09 | 50854 |
52 | XeLaTeX 한글 문서에서 한자와의 베이스라인 맞추기 | karnes | 2011.07.06 | 5409 |
51 | \vsplit : 수직 박스 쪼개기 [6] | DohyunKim | 2011.07.06 | 6563 |
50 | \unexpanded 와 etoolbox 패키지 | DohyunKim | 2011.07.05 | 4525 |
49 | hcrnumbers를 oblivoir enumerate에서 쓰기 | karnes | 2011.06.29 | 4526 |
48 | 함초롬의 nalt 속성을 이용해보자 | DohyunKim | 2011.06.29 | 4440 |
47 | \XeTeXinputnormalization | DohyunKim | 2011.06.27 | 4413 |
46 | jbig2 사용법 | DohyunKim | 2011.06.22 | 4777 |
45 | graphicx 패키지 안 쓰고 이미지 삽입하기 | DohyunKim | 2011.06.22 | 4723 |
44 | $$ 를 쓸까? \[ 를 쓸까? | DohyunKim | 2011.06.21 | 4454 |
43 | 파이 구하기: 시각화 버전 | karnes | 2011.06.18 | 4658 |
42 | 난수로 괘를 뽑아보았다 | karnes | 2011.06.15 | 4403 |
» | 몬테카를로 기법으로 파이 구하기: 정수 연산 버전 | DohyunKim | 2011.06.15 | 5006 |
40 | "a unit of 65536" | karnes | 2011.06.14 | 4415 |