clc
syms a;
syms b;
syms c;
A=[a a^2 a^3-1;b b^2 b^3-1;c c^2 c^3-1]
det(A)