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