det(Matrix([[(b+c)**2, c**2, b**2], [c**2, (a+c)**2, a**2], [b**2, a**2, (a+b)**2]]))