\(\begin{vmatrix}x^2&x^2-(y-z)^2&yz\\y^2&y^2-(z-x)^2&zx\\z^2&z^2-(x-y)^2&xy\end{vmatrix}\)
R1→R1 - R2, R2→R2 - R3
= \(\begin{vmatrix}(x+y)(x-y)&(x+y)(x-y)-(x-y)(2z-x-y)&-z(x-y)\\(y+z)(y-z)&(y+z)(y-z)-(y-z)(2x-y-z)&-x(y-z)\\z^2&z^2-(x-y)^2& xy\end{vmatrix}\)
R1→ R1/(x - y),R2→R2/(y - z)
= (x - y) (y - z)\(\begin{vmatrix}x+y&2(x + y + z)&-z\\y+z& 2(x +y+z)& -x\\z^2&z^2-(x-y)^2&xy\end{vmatrix}\)
R2→ R2 - R1
= (x - y) (y - z) \(\begin{vmatrix}x+y&2(x+y+z)&-z\\z-x&0&z-x\\z^2&z^2-(x-y)^2& xy\end{vmatrix}\)
R2 → \(\frac{R_2}{z-x}\)
= (x - y) (y -z)(z - x)\(\begin{vmatrix}x+y&2(x+y+z)&-z\\1&0&1\\z^2&z^2-(x-y)^2&xy\end{vmatrix}\)
C3→ C3 - C1
= (x - y)(y - z) (z - x)\(\begin{vmatrix}x+y&2(x+y+z)&-(x+y+z)\\1&0&0\\z^2&z^2-(x-y)^2&xy-z^2\end{vmatrix}\)
Expand determinant along C1
= (x - y) (y - z) (z - x) (2(x + y + z) x (xy - z2) + (x + y + z)(z^2-(x - y)2))
= -(x - y)(y - z) (z - x) (x + y + z)(2xy - 2z2 + z2 - x2 - y2 + 2xy)
= -(x - y) (y - z) (z - x) (x + y + z)(-(x2 + y2 + z2) + 4xy)
= (x - y) (y - z)(z - x) (x + y + z) (x2 + y2 + z2 - 4xy)