solve(c=a*(a^(2)-3*b^(2)) and d=(3*a^(2)-b^(2))*b, {a,b})