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 으로
컴파일하면 된다.
