I am working on a program to calculate steady-state probability in a given Markov chain. My code is pasted as bellow:

CalculateSteadyStateProbability := proc (m::integer, r::integer, CV, ROU) 

local P, solutions, pai, l, eqns; 

P := CalculateMatrix(m, r, CV, ROU);    #here we calculate a (m+r)x(m+r) transition martix

eqns := {`$`('pai[i+1] = add(pai[j+1]*P[j+1, i+1], j = 0 .. m+r)', i = 0 .. m+r)

 ...