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. OK, I Understand
 
Top