A:-matrix([x−1,1,0],[0,y−1,−1],[−1,0,z−1]);
B:-matrix([x,1,0],[0,y−1,−1],[−1,0,z−1]);
C:-matrix([x,1,0],[0,y,−1],[−1,0,z]);
D:-matrix([z,1,0],[0,z,−1],[−1,0,z]);
log(determinant(A))+log(determinant(B))+log(determinant(C))+log(determinant(D));
J:jacobian([%],[x,y,z])$
factor(subst(x, x^2, subst(y, y^2, subst(z, z^2, subst(x, x^3, subst(y, y^3, subst(z, z^3,
subst(x, x^4, subst(y, y^4, subst(z, z^4, subst(x, x^5, subst(y, y^5, subst(z, z^5, expand(%))))))))))))))$
factor(subst(x, x^2, subst(y, y^2, subst(z, z^2, subst(x, x^3, subst(y, y^3, subst(z, z^3,
subst(x, x^4, subst(y, y^4, subst(z, z^4, subst(x, x^5, subst(y, y^5, subst(z, z^5, expand(%))))))))))))));
J.J=0$
factor(subst(x, x^2, subst(y, y^2, subst(z, z^2, subst(x, x^3, subst(y, y^3, subst(z, z^3,
subst(x, x^4, subst(y, y^4, subst(z, z^4, subst(x, x^5, subst(y, y^5, subst(z, z^5, (
subst(x, x^6, subst(y, y^6, subst(z, z^6, subst(x, x^7, subst(y, y^7, subst(z, z^7,
subst(x, x^8, subst(y, y^8, subst(z, z^8, subst(x, x^9, subst(y, y^9, subst(z, z^9, (

%)))))))))))))))))))))))))))$
factor(subst(x, x^2, subst(y, y^2, subst(z, z^2, subst(x, x^3, subst(y, y^3, subst(z, z^3,
subst(x, x^4, subst(y, y^4, subst(z, z^4, subst(x, x^5, subst(y, y^5, subst(z, z^5, (
subst(x, x^6, subst(y, y^6, subst(z, z^6, subst(x, x^7, subst(y, y^7, subst(z, z^7,
subst(x, x^8, subst(y, y^8, subst(z, z^8, subst(x, x^9, subst(y, y^9, subst(z, z^9, (
subst(x, x^10, subst(y, y^10, subst(z, z^10, (
%)))))))))))))))))))))))))))))));