solve([a^2+b*c=a,a*b+b*d=b,a*c+c*d=c,d^2+b*c=d],[a,b,c,d])