Help:=proc(): print(` Eq(n,a,x) , Sol(n,a), MTT(n,a), Arthur(n),Cayleynx(n,x) `): print(`Cayley(n) `): end: eq3:= {(a[2,1]+a[3,1])*x[1]-a[1,2]*x[2]-a[1,3]*x[3], (a[1,2]+a[3,2])*x[2]-a[2,1]*x[1]-a[2,3]*x[3], (a[1,3]+a[2,3])*x[3]-a[3,1]*x[1]-a[3,2]*x[2]}: Eq:=proc(n,a,x) local i,j: {seq((add(a[j,i],j=1..i-1)+add(a[j,i],j=i+1..n) )*x[i]- add(a[i,j]*x[j],j=1..n)+a[i,i]*x[i],i=1..n)} : end: Var:=proc(n,x) local i: {seq(x[i],i=1..n)}:end: Sol:=proc(n,a) local x,sol: sol:=solve(Eq(n,a,x),Var(n,x)): [seq(normal(subs(sol,x[i])),i=1..n)]: end: with(linalg): #MTT: Matrix Tree Theorem MTT:=proc(n,a) local i,j,M: M:=[seq([seq(-a[i,j],j=2..i-1), add(a[j,i],j=1..n)-a[i,i], seq(-a[i,j],j=i+1..n)], i=2..n)]: det(convert(M,matrix)): end: Arthur:=proc(n) local a,i: a:= [seq([1$(i-1),0,1$(n-i)],i=1..n)]: MTT(n,a): end: Cayley:=proc(n) local a,i: a:= [seq([-1$(i-1),n,-1$(n-i)],i=1..n)]: det(convert(a,matrix)): end: Cayleynx:=proc(n,x) local a,i: a:= [seq([-1$(i-1),x,-1$(n-i)],i=1..n)]: factor(det(convert(a,matrix))): end: