a=var('a')
b=var('b')
c=var('c')
A1 = matrix([[1,a,a^2],[1,b,b^2],[1,c,c^2]])
det(A1)