• API
• FAQ
• Tools
• Archive
SHARE
TWEET # c6 tema 2 a guest May 26th, 2019 86 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. cerinta6 = function(N, pE, pV, pS, pN)
2. {
3.  atleast25 = 0
4.  for(i in 1:1000)
5.  {
6.    mat = matrix(0, nrow = N, ncol = N, byrow = T)
7.    mat[40, 40] = 1
8.
9.    ard = 1
10.    auars = 0
11.
12.    while(ard != 0)
13.    {
14.      for(i in 1:N)
15.      {
16.        for(j in 1:N)
17.        {
18.          if(mat[i, j] == 1)
19.          {
20.            prob = runif(1,0,1)
21.
22.            if(prob <= pE && j > 1 && mat[i, j - 1] == 0)
23.            {
24.              mat[i, j - 1] = 1
25.              ++ard
26.            }
27.
28.            prob = runif(1, 0, 1)
29.
30.            if(prob <= pV && j < N && mat[i, j + 1] == 0)
31.            {
32.              mat[i, j + 1] = 1
33.              ++ard
34.            }
35.
36.            prob = runif(1, 0, 1)
37.
38.            if(prob <= pS && i > 1 && mat[i - 1, j] == 0)
39.            {
40.              mat[i - 1, j] = 1
41.              ++ard
42.            }
43.
44.            prob = runif(1, 0, 1)
45.
46.            if(prob <= pN && i < N && mat[i + 1, j] == 0)
47.            {
48.              mat[i + 1, j] = 1
49.              ++ard
50.            }
51.
52.            mat[i, j] = 2
53.            auars++
54.            --ard
55.          }
56.        }
57.      }
58.    }
59.
60.    if(auars * 4 >= N * N)
61.    {
62.      atleast25 = atleast25 + 1
63.    }
64.  }
65.
66.  print(atleast / 1000)
67. }
68.
69. cerinta6(40, 0.7, 0.4, 0.4, 0.2)
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top