First hitting time python

from numpy import *

hit_idx = (0,4)

# Define a graph by edge list
edges = [[0,1],[1,2],[2,3],[2,4]]

# Create adj. matrix
A = zeros((5,5))
A[zip(*edges)] = 1
# Undirected condition
A += A.T

# Make the final state an absorbing condition
A[hit_idx[1],:] = 0
A[hit_idx[1],hit_idx[1]] = 1

# Make a proper Markov matrix by row normalizing
A = (A.T/A.sum(axis=1)).T

B = A.copy()
Z = []
for n in xrange(100):
    Z.append( B[hit_idx] )
    B = dot(B,A)

from pylab import *
plot(Z)
xlabel("steps")
ylabel("hit probability")
show()

Source: stackoverflow.com

Similar codes

0 vote194 views1 answer
Python function pointer with multiple args
Python
1 vote156 views1 answer
Get urls from any website use python
Python
6 votes203 views4 answers
Datetime python timezone
Python
0 vote195 views1 answer
Print less than specific number in one row python
Python
0 vote255 views1 answer
Pydrive list folders
Python
0 vote202 views1 answer
Num = [7,8, 120, 25, 44, 20, 27] newnum = [] def remove_even(num): for i in num: if i%2 != 0: newnum.append(i) return newnum print("get_unevens") test(remove_even(num), [7,25,27])
Python
0 vote222 views1 answer
Python get factors of a number
Python
22 votes195 views3 answers
How to sort list in descending order in python
Python
0 vote238 views1 answer
Print A to z vy using loop in python
Python
2 votes191 views1 answer
Django models distinct
Python