RCCOUT3.mw

> ############ START HERE ######################

> Delta:=evalm((R-B)/2):A:=evalm((R+B)/2): `R`=matrans(R);`B`=matrans(B);

R = vector([4, 3, 1, 2])

B = vector([3, 4, 4, 3])

        AA is (numerical) Asup2. Use a2 for numerical A_2, j2a for numerical J2-AA

> A2:=evalm((1/2)*(R2+B2)):a2:=evalm(AA+Del2):iszero(A2-a2),"det",det(evalm(J2-AA));

true,

> print(sympow(R,2));print(sympow(B,2));

matrix([[0, 0, 0, 0, 0, 1], [0, 0, 1, 0, 0, 0], [0, 0, 0, 0, 1, 0], [0, 1, 0, 0, 0, 0], [0, 0, 0, 1, 0, 0], [1, 0, 0, 0, 0, 0]])

matrix([[0, 0, 0, 0, 0, 1], [0, 0, 0, 0, 0, 1], [0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1], [0, 0, 0, 0, 0, 1]])

> j2a:=evalm(J2-AA):evalm(j2a),det(j2a),evalf(det(j2a));

matrix([[1, 0, 0, 0, 0, (-1)/2], [0, 3/4, (-1)/4, 0, 0, (-1)/4], [0, 0, 1, (-1)/4, (-1)/4, (-1)/4], [0, (-1)/4, (-1)/4, 1, 0, (-1)/4], [0, 0, 0, (-1)/4, 3/4, (-1)/4], [(-1)/4, (-1)/4, 0, 0, (-1)/4, 3/...

> print("LEFT NULL",Q2,"DELTA2",Del2,"RIGHT NULL",P2);"DET",det(Del2);

> rdd:=read2(Del2):for i to nops(rdd[1]) do print(convert(choose(n,2)[i],set),i,rdd[1][i],"               ",rdd[2][i]) od;

{1, 2}, 1, [[6, 6], [

{1, 3}, 2, [[3, 6], [2,

{1, 4}, 3, [[5,

{2, 3}, 4, [[2,

{2, 4}, 5, [[4, 6], [5,

{3, 4}, 6, [[1, 6], [2, 5]],

        j2a2 is the numerical I-A_2

> j2a2:=evalm(j2a-Del2):"A_2",evalm(a2),"I-A_2",evalm(j2a2)," DET ",det(j2a2);

> multiply((a2),uu);

vector([1, 1, 1/2, 1/2, 1, 1])

> unassign('x'):pi2:=evalm(1/x[1]*(linsolve(transpose(j2a2),vector(NN,0),'r',x)));u2:=evalm(1/x[1]*(linsolve(j2a2,vector(NN,0),'r',x)));

pi2 := vector([1, 0, 0, 0, 0, 2])

u2 := vector([3, 2, 1, 1, 2, 3])

                                OMEGA2 and E's coming up

> Omega2:=abel(a2);readcycles(Omega2);abel(xtend(a2));#THIS LAST INCLUDES THE ABSORBING STATE

Omega2 := matrix([[1/3, 0, 0, 0, 0, 2/3], [2/9, 0, 0, 0, 0, 4/9], [1/9, 0, 0, 0, 0, 2/9], [1/9, 0, 0, 0, 0, 2/9], [2/9, 0, 0, 0, 0, 4/9], [1/3, 0, 0, 0, 0, 2/3]])

{{1, 6}}

matrix([[1/3, 0, 0, 0, 0, 2/3, 0], [2/9, 0, 0, 0, 0, 4/9, 1/3], [1/9, 0, 0, 0, 0, 2/9, 2/3], [1/9, 0, 0, 0, 0, 2/9, 2/3], [2/9, 0, 0, 0, 0, 4/9, 1/3], [1/3, 0, 0, 0, 0, 2/3, 0], [0, 0, 0, 0, 0, 0, 1]]...