with(linalg):

lambdaA := 0.05;
muA := 0.02;
lambdaB := 0.1;
muB := 0.05;

size := 5;

Q := array(sparse,1..size,1..size);
e := array(sparse,1..size);

Q[1,2] := lambdaA;

Q[1,3] := lambdaB;

Q[2,1] := muA;

Q[2,4] := lambdaB;

Q[3,1] := muB;

Q[3,5] := lambdaA;

Q[4,3] := muA;

Q[5,2] := muB;

for i to size do 
    s := 0:
    for j to size do 
        s := s + Q[i,j]
    od;
    Q[i,i] := -s
od;

Qt := transpose(Q);

for i to size do Qt[size,i] := 1 od;

e[size] := 1;

p := linsolve(Qt,e);