Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- n=9999
- liczbaPi = function(n)
- {
- #ustawiamy wartosci poczatkowe dla zmiennych x i y
- x=0
- y=0
- #rysujemy wykres i ustawiamy jego ograniczenia (tez poto aby wyczyscic pole z wykresami)
- plot(x,y,xlim=c(-0.2,2.2), ylim=c(-0.2,2.2))
- #musimy miec zmienna liczaca ile jest puntkow spelniajacych zalozenie
- i=0
- #wykonujemy petle generujac dowolne punkty z przedzialu 0-2, sprawdzamy czy spelniaja zalozenie
- #nastepnie rysujemy wykres tych punktow
- for(m in 1:n){
- x = runif(1,0,2)
- y = runif(1,0,2)
- #warunek to zwykle rownanie kola tylko musielismy przesunac jego srodek do srodka kwadratu.
- if((x-1)^2 + (y-1)^2 <= 1){
- #inkrementujemy i(ilosc punktow spelniajacych zalozenie)
- i=i+1
- #dorysowujemy punkty
- points(x,y,col="green")
- }
- else
- {
- points(x,y,col="red")
- }
- }
- #zmienna i przybliza nam pole kola, a ilosc wszystkich punktow pole kwadratu
- #stosunek "pola" kola wpisanego do "pola" calego kwadratu da nam PI/4. Wystarczy tylko pomnozyc przez 4 i wychodzi PI
- return((i/n)*4)
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement