program hirarchichal_detreministic implicit none integer, parameter:: N=4 !initial core integer, parameter:: ns=1 !number of steps integer, parameter:: nn = n**(ns+1) integer, dimension (1:nn,1:nn):: a !adjacency matrix integer:: i, j, k, ii, m, l real::p(nn) open(1,file='out') a = 0 print*,sum(a), a(0,0) pause !****************************************** !initial fully connected network; ns=0 a(1:n,1:n)=1 do i=1,n a(i,i)=0 enddo !****************************************** print*,sum(a), a(0,0), a(17,170) do i=1,ns do k= 2, n-1 do ii= k*(n**i)+1, (k+1)*(n**i) do j=k*(n**i)+1, (k+1)*(n**i) !print*,i,k,ii,j l = mod(ii,n**i+1) m = mod(j,n**i+1) a(ii,j) = a(l,m) print '(20(I2,2x))',sum(a), i, k, ii, j, l, m enddo ! pause enddo enddo enddo print*,"fd" Do i=1, nn !k=sum(a(1:nn,i)) !p(k)=p(k)+1. do j=1,nn write(1,*)i,j, a(i,j) enddo enddo p=p/sum(a) print*,"fd/" do i=1,nn !write(1,*)i,p(i) enddo end