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 vote164 views1 answer
Pandas read_csv drop column
Python
2 votes147 views2 answers
How to add time with time delta in python
Python
1 vote145 views1 answer
Array of dic
Python
3 votes194 views3 answers
How to load images from folder in python
Python
32 votes165 views4 answers
Copy to clipboard python
Python
0 vote150 views1 answer
How to print a matrix in python
Python
91 votes149 views6 answers
Time.sleep
Python
1 vote166 views4 answers
How to input in a list in python
Python
1 vote164 views1 answer
Python regex get string before character
Python
1 vote162 views1 answer
Olst = [] a = int(input()) b = int(input()) for ele in range(a,b+1): if ele%2 != 0: olst.append(ele) print(olst[::-1])
Python