> | ############ START HERE ###################### |
> | Delta:=evalm((R-B)/2):A:=evalm((R+B)/2): `R`=matrans(R);`B`=matrans(B); |
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)); |
> | print(sympow(R,2));print(sympow(B,2)); |
> | j2a:=evalm(J2-AA):evalm(j2a),det(j2a),evalf(det(j2a)); |
> | 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; |
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); |
> | 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))); |
OMEGA2 and E's coming up
> | Omega2:=abel(a2);readcycles(Omega2);abel(xtend(a2));#THIS LAST INCLUDES THE ABSORBING STATE |